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Preface 


The major purpose of this book is to introduce and investigate various topics 
related to so-called generalized trigonometric and hyperbolic functions. The 
methodology and the associated analysis are essentially mine and in many 
cases are not connected mathematically with previous work on this subject. 
In general, our obtained results are derived and discussed using a “rigorous 
heuristic” style of mathematical analysis. However, in spite of what some 
might consider a limitation, this procedure allows very interesting results to 
follow. The background knowledge needed to study with profit and to under- 
stand the contents of the book consists of the reader having taken courses in 
elementary plane geometry, trigonometry, and a year of calculus. 

My interests in generalized trigonometric and hyperbolic functions had 
their genesis nearly four decades ago when I simultaneously began the study 
of and the writing of a book on nonlinear oscillations: 


e R.E. Mickens, Introduction to Nonlinear Oscillations (Cambridge Uni- 
versity Press, New York, 1981). 


For many conservative one-degree-of-freedom systems, the energy function, 
i.e., the Hamiltonian function, takes the form 


2 


H(c,y)= +2, n=0,1,2,..., (0.1) 
where dx(t) 
Xx 
y(t) = ip a(0)=1, y(0)=0, (0.2) 


and we have used scaling to eliminate unneeded parameters, and to normalize 
the energy to H(z, y) = 1. A general extension of Equations (0.1) and (0.2) is 


A(a,y) = |yl? + lal? = 1, 
x(0l=1, y(0)=0, (0.3) 
p>0, q>0, 


where p and q are real. Since in the xz-y plane the solution curves, y = y(), 
are simple, closed curves, then all the solutions, (x(t), y(t)), are periodic. It 
is important to note that in the first of Equations (0.3), the absolute value 
appears. This requirement is critical for our work, otherwise, the curves may 
not be simple and closed. 


xill 


xiv Preface 


Our next spark of interest occurred with the reading of a paper by 
Schwalm, i.e., 


e W. Schwalm, Elliptic functions sn, cn, dn as trigonometry: 
http://www.und.edu/instruct/schwalm/MAA_ Presentation_10_02/ 
handout.pdf, 2005. 


This work clearly indicated the potential of using planar geometrical con- 
cepts to analyze two-dimensional dynamic systems in such a manner that new 
classes of periodic functions could be created. This book illustrates what can 
be done using this methodology. 

Another important point is that the results reported in this book are 
greatly constrained by the condition that |y|? +|2|? = 1 is invariant under the 
transformations 


Ti:@->-2, yoy 
To:272, yory-y (0.4) 
T3 = T To = ToT, : 2 > —2, Umea 


One consequence is that periodic and non-periodic solutions have many of 
the same general mathematical features as those possessed by the standard 
trigonometric and hyperbolic functions. 


The book consists of ten chapters, an appendix, and a short bibliography. 
Chapters 1 and 2 restate and reformulate the standard results known, re- 
spectively, on the trigonometric cosine and sine functions, and the associated 
hyperbolic functions. The techniques used for this analysis are then applied, 
in Chapters 3 and 4, to the “square” and “parabolic” periodic and hyperbolic 
functions. Our major, original contribution is to demonstrate the existence of 
periodic hyperbolic-type functions. 

Chapter 5 provides a detailed discussion of a particular class of periodic 
solutions to the functional equation 


f(0)=1, g(0) =0. (0.5) 


A rather cursory presentation is given in Chapter 6 on some of the previous 
work done on the so-called generalized trigonometric functions. Most of these 
formulations were done within the framework of integral relations. We provide 
references to publications which are of direct relevance to the general topics 
presented in this book. 

Chapters 7 and 8, while the shortest chapters, are the central features of 
the book. They define a new class of generalized trigonometric and hyperbolic 
functions based on the functional equations 


yl? + lel? =1; 2(0)=0, y(0) =0; (0.6) 


Preface Xv 


ly? — |x|? =1;  2(0) =0,y(0) =1, (0.7) 
p>0, q>0. (0.8) 


Chapter 9 gives an analysis of several cases which model important dy- 
namic systems involving, in general, nonlinear oscillations. For these situations 
p=2andq>1. 

Chapter 10 lists and briefly discusses a number of unresolved issues in the 
field of generalized trigonometric and hyperbolic functions. 

The Appendix includes several mathematical listings of certain special 
functions and their properties, along with related materials on Fourier series, 
the delta and sign functions, and a brief summary of Hamiltonian dynamics. 

Finally, the short Bibliography gives a listing of books, publications, and 
websites which were helpful to the writing of this volume. 


As always, I am truly thankful to Annette Rohrs for her outstanding tech- 
nical and editorial services, which allowed my handwritten pages to be trans- 
formed into this volume. Both she and my wife, Maria Mickens, provided valu- 
able assistance and encouragement to successfully finish this writing project 
on time. I acknowledge and thank the editor, Sarfraz Khan, for his help in the 
initiation and production of this book. Further thanks and great appreciations 
are extended to my former student and now colleague, Kale Oyedeji (More- 
house College, USA), and my colleague and friend, Ivana Kovacic (University 
of Novi Sad, Serbia). Both have provided aid and insights that have enhanced 
my search for deep understandings of a range of issues in nonlinear dynamics 
and special functions. 

Finally, it should be noted that my research was especially enhanced by 
the efforts of two persons, both reference librarians at the Robert Woodruff 
Library, Atlanta University Center: Imani Beverly and Bryan Briones. They 
obtained needed copies of articles and books, often within days of my requests, 
and in addition, provided general cheerful and professional services. 


Ronald E. Mickens 
Atlanta, Georgia, USA 
August 2018 


Taylor & Francis 


Taylor & Francis Group 
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Chapter 1 


TRIGONOMETRIC AND 
HYPERBOCLIC SINE AND 
COSINE FUNCTIONS 


1.1 INTRODUCTION 


Historically, the trigonometric sine and cosine functions had their genesis 
in the study of properties of the magnitude of line segments contained within 
the perimeter of a circle. This chapter examines these functions from several 
different perspectives and ends with the introduction of a new set of functions 
which we denote the 6-periodic functions. 

Section 1.2 starts with the standard definition of the sine/cosine functions 
based on the geometric properties of the circle. Section 1.3 is based on an 
analytic definition of sine/cosine which follows from the Euler relation 


e” = cos + isin 9. (1,1) 


Another interesting and important way of examining the properties of 
sine/cosine functions is to consider them solutions to a dynamic system mod- 
eling the simple harmonic oscillator. These results are presented in Section 1.4. 
The next three sections are devoted to the associated hyperbolic functions. In 
the simplest, geometric case, this is accomplished by replacing @ by 70, where 
i = /—1. This is done in Section 1.5. The hyperbolic functions can also be 
examined by considering them to be solutions to a particular dynamic system, 
and these considerations are given in Section 1.6. While the usual, standard 
trigonometric hyperbolic functions are not periodic, we show in Section 1.7 
how to construct periodic generalizations; we call them 6-periodic hyperbolic 
functions. Finally, we conclude the chapter with a summary of the obtained 
results, along with the inclusion of comments on what these results (might) 
mean. 


2 GENERALIZED TRIGONOMETRIC & HYPERBOLIC FUNCTIONS 


FIGURE 1.2.1: Plot of circle x? + y? = 1 and the definitions of cos@ and 
sin 0, respectively as «/r and y/r. 


1.2 SINE AND COSINE: GEOMETRIC 
DEFINITIONS 


The unit circle, in the x-y plane, is defined by the formula 
oe +y?=1, (1.2) 


Let a ray be drawn from the origin O, such that it intersects the perimeter of 
the circle at point P(x, y); see Figure 1.2.1. Denote by @, the angle that this 
ray makes with the positive z-axis. The sine and cosine functions are defined 


to be 
sin@d = 4 
e : (1.3) 
cos 0 = = 


TRIGONOMETRIC & HYPERBOLIC SINE & COSINE FUNCTIONS 3 
Since r = 1, for the unit circle, these definitions reduce to 
sind=y, cos@=«a. (1.4) 


Taking into consideration the geometry in Figure 1.2.1 and further taking 
into consideration the fact that « and y are “signed” (positive or negative) 
real numbers, the following is a listing of some of the major properties of the 
sine and cosine functions: 


(cos 0)? + (sind)? = 1 (1:5) 
cos(—@) = cos@, sin(—@) = — sind (1.6) 
—1<|cos@}<+1, -1< |sin@| <+4+1 (1.7) 
cos(? + 27) =cos@, sin(@+ 27) = sind (1.8) 
cos(0)=1, cos(¥)=0, cos(m) = —1 
- (1.9) 
cos (37) =0, cos(27) = 1 
ont) = i sin (Z) = sin(7) = 0 (1.10) 
sin (34) =—-1,  sin(27) =0 
(+), 0<0<a7 
0) = 1.11 
sored) fa T<O0< 20 = 
(+), 0<0< § 
sgn(cos0) = 4 (-), 3<0< (1.12) 
(+), %<0<29 
Max) /|cos 6| Min \ /|sin 0| 
= 1.1 
os on A te s 0 a) 


1.3 SINE AND COSINE: ANALYTIC DEFINITION 


An alternative way to introduce and investigate the sine and cosine func- 
tions is to start with the Euler formula 


e” = cos6 + isin 9. (1.14) 


Note that this is the same as Equation (1.1), but for convenience, we rewrite 
it again, with a new numbering, for ease of reference to it in this section. 
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FIGURE 1.2.2: Plots of cos@ and sin@ vs @. 
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Taking the complex conjugate of Equation (1.14) and multiplying it by 
Equation (1.14) gives 


(ec) (e”) = (cos — isin) (cos 6 + isin @) 


1 = (cos 6)? + (sin 6)’, (1.15) 
and this is just the previously obtained result stated in Equation (1.5). 
Since 
(9) = ery, (1.16) 
it follows that 
cos(—@) + isin(—0) = cos@ — isin 9, (1.17) 
and we conclude that 
cos(—0) = cos@, sin(—0@) = —sin@, (1.18) 


i.e., the cosine and sine functions are, respectively, even and odd functions of 
0. 

In the next several subsections, we derive other properties of the co- 
sine/sine functions. 


1.3.1 Derivatives 


Taking the derivatives of both sides of Equation (1.14) gives 


d , 4 dcos@  _.dsin@ 
—/(e!”) = La 
wo = ae + aa on) 
But 
(eit) = ie!” 
= —sin#@+icos6, (1.20) 


and comparing the corresponding real and imaginary parts of Equations (1.19) 
and (1.20), allows the following conclusions to be made 


dcos@ dsin 0 
=—gj — 50. 1.21 
77 sin 0, 7 cos 6 (1.21) 


If now the derivative is taken of the results expressed in Equation (1.21), we 
obtain 


d?cos@ ___—dsin# 
doz dé 


d*sin@  dcosé ; 
ae = = aOs = —sind. (1.22b) 


=—cos6, (1.22a) 
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Consequently, cos@ and sin@ both satisfy the following second-order, linear 
differential equation, with constant coefficients 


d?w(0) 
—— 0) =0. 1.23 
+ wl) (1.23) 

Observe that the general solution for Equation (1.23) is 
w(0) = Acos@ + Bsin9, (1.24) 


where A and B are arbitrary real constants. Thus, the cosine and sine functions 
are solutions to the following initial value problems to Equation (1.24): 


cos@: w(0) = 1, auto) =0, (1.25a) 
sin8 : w(0) = 0, aut) =1. (1.25b) 


1.3.2 Integrals 


It follows directly from the results presented in Equation (1.2) that 
[soa =—cos6, (1.26a) 
i cos 6 d@ = sin 0. (1.26b) 


1.3.3 Taylor Series 


The exponential function exp(u) has the Taylor series 


= 5. a (1.27) 
k=0 
Therefore, for u = 76, we have 
io _ <> (i8)* 
eS aT (1.28) 
k=0 
If we now use 
(i)° = 1, (i)* = t, (i)? = =; (i)? = =A; @) = 1, (1.29) 


' 2k ; 
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Comparing Equations (1.14) and (1.30) gives the following expressions for the 
Taylor series of the cosine and sine functions 


co g2k 
cos 9 = dL) aa 
Sas (1.31) 
sind = SO Gea 
Spe eae (1.32) 


1.3.4 Addition and Subtraction Rules 


We now derive formulas for the cosine and sine functions corresponding to 
the sum and difference of two angles. This can be easily done by starting with 
the relation 

(#1) (e%2) — e(91 +62) (1.33) 


Expanding the terms gives 
(ec?) (e%2) = (cos, + isin 61)(cos 2 + isin 62) 
= (cos 6, cos 02 — sin 6; sin 2) 
+ i(cos 61 cos 2 + sin 6; cos 02), (1.34) 
and 
e'(91+92) — cos(O, + 02) + isin(O, + 42). (1.35) 
Comparison of the real and imaginary parts of Equations (1.34) and (1.35) 
gives 
cos(6; + 62) = cos 6; cos 62 — sin 6; sin 02, (1.36a) 
sin(@; + 62) = sin 6; cos 62 + cos 6; sin 02, (1.36b) 


and these are called the addition relations for the sine and cosine functions. 
Changing the sign of 62, i.e., 02 + —02, gives the subtraction relations 


cos(6; — 62) = cos 6; cos 62 + sin 6; sin 02, (1.37a) 
sin(@; — 62) = sin 1 cos 62 — cos 6; sin 02. (1.37b) 


1.3.5 Product Rules 
Adding Equations (1.36a) and (1.37a) gives 


1 
cos 61 cos 82 = (5) [cos(0; + 62) + cos(@1 — @2)]. (1.38) 
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Likewise, adding Equations (1.36b) and (1.37b) gives 

sin 0) cos 62 = (5) [sin(O, + 62) + sin(9, — 2)]. (1.39) 
Similarly, subtracting Equation (1.36a) from Equation (1.37a) gives 

sin 0; sin 02 = (5) [cos(, — 62) — cos(; + 02)]. (1.40) 


The formulas given in Equations (1.38), (1.39), and (1.40) are the product 
relations for the cosine and sine functions. 


1.4 SINE AND COSINE: DYNAMIC SYSTEM 
APPROACH 


A third way of investigating the cosine and sine functions is to consider 
them to be solutions to a dynamic system. For this situation, the dynamic 
system turns out to be the undamped harmonic oscillator studied in every 
elementary physics course. 

It is important to note that the independent variable is now the “time,” 
rather than the angle 6. 

To begin, let Equation (1.2) be rewritten to the form 


H(z,y) = (5) (a? + y’) = > (1.41) 


where H(az,y) is the Hamiltonian of the system. This particular functional 
structure is a rescaled version of the harmonic oscillator energy function 


Ae a x DA 28 1 2 

A(z,y)=—+ (5) kin (5) kA’, (1.42) 
where Z = z(t) and y = y(t), the initial conditions are 7(0) = A and y(0) = 0, 
and m is the mass of the object subject to a linear spring characterized by a 
spring constant k. The variable Z(t) is the location of the mass m at time f, 
and y(t) = mdz(t)/dt is the momentum. 

Returning to Equation (1.41), the dynamics is determined by the following 

pair of first-order differential equations 


dx OH(x,y) dy OH(x,y) 


RS yee 1.48 
dt Oy ° dt Ox (Ee) 
which upon evaluation gives 
dx dy 
eal —~=-_y, 1.44 
ae es ape (1.44) 
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These two differential equations are equivalent to the following single, second- 
order equation 


a +a2=0, (1.45) 
for which the general solution is 
x(t) = Ci cost + Cy sint. (1.46) 
Using the initial conditions 
x(0)=A=1, y(0)=0, (1.47) 


Cy=1, C.=0 (1.48) 


and 
x(t) =cost, y(t) =—sint. (1.49) 


We thus see that the cosine and sine functions appear as solutions to the 
harmonic oscillator dynamic system. 


1.4.1 «-y Phase-Space 


Let us now examine in more detail the phase-space for the harmonic os- 
cillator. This phase-space is the set of all trajectories in the z-y plane. These 
curves satisfy a first-order differential equation and we now provide two related 
ways to derive this equation. 

First, note that the curves are given by y = ¢(x) and taking its time 
derivative gives 


dy _dédx _ dy de 


dt da dt dx dt’ ee) 
which can be solved for dy/dx to obtain the result 
dy  dy/dt 
a . 1.51 
dx  da/dt ay 


From Equation (1.44), it follows that 


at --(). (1.52) 


A second way to derive Equation (1.52) is to take the derivative of the 
Hamilton, Equation (1.41), with respect to x, where y = y(x), to obtain 


d dy 
ag (& y(@)) =E+yT = 0, (1.53) 


which upon solving for dy/dx gives Equation (1.52). 


10 GENERALIZED TRIGONOMETRIC & HYPERBOLIC FUNCTIONS 


The fixed point, (a*, y*), are constant solutions and they can be calculated 
from the conditions 


dx dy 
ale —=0. 1.54 
dt m dt (eos) 
Referring to Equation (1.44), we see there is only one fixed point, i-e., 
x*(t)=0, y*(t)=0. (1.55) 


For the harmonic oscillator, the fixed point represents the rest state of the 
oscillator. 


1.4.2 Symmetry Properties of Trajectories in Phase-Space 


Inspection of either the Hamiltonian, Equation (1.41), or the differential 
equation determining the trajectories in phase-space, Equation (1.52), shows 
that both are invariant under the following transformations 


Ty:t@7>-2@, yoy, (1.56) 
To:072, yrr-y, (1.57) 
T3 =TTn: "44-2, yo-y. (1.58) 


These transformations have the following interpretations (see Figure 1.4.1): 
(i) T; corresponds to reflection of a point P through the y-axis. 
(ii) T2 is reflection through the z-axis. 
(iii) T; = T,T> is reflection through the origin. 
Figure 1.4.2 shows the results of these transformations on segments of curves 
in phase-space. 


1.4.3. Null-Clines 
For the Hamiltonian H (a, y), the equations of motion are 


dx _ OH (x,y) 


2 = Rey) = g(x,y). (1.60) 


Thus, the differential equation for the trajectories in phase-space is 
dy _ dy/dt _ g(x,y) 


dx dx/dt  f(a,y) 


(1.61) 


The null-clines are curves in phase-space where the derivative, dy/dz, is either 
zero or unbounded (infinite). They are determined by the requirements 


w Ha deeaelayy =o. (1.62a) 
OY — 00 : f(,yoo()) =0. (1.62b) 
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Yo T:2> ayy 
et P 
O z 
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FIGURE 1.4.1: Symmetry transformations, in the x-y phase-plane. 
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14 1 
Sg a 
Ty 
1 
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1’ 


FIGURE 1.4.2: Symmetry transformations for segments of curves. 
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Applying these considerations to the harmonic oscillator, we have 


1 dx dy 

A(2,y) = (5) Ce) oe (1.63a) 
dy x 
= 1.63b 
rE (1.63b) 
d 
~ = 0:2 =0 or along the y-axis, (1.63c) 
d 
a = oo: x =0or along the z-axis. (1.63d) 


1.4.4 Geometric Proof that All Trajectories Are Closed 


We now demonstrate that just using knowledge of the major features of the 
trajectories in phase-space, it follows from the differential equation dy/dx = 
(—ax/y) that all the trajectories for the harmonic oscillator are simply closed 
curves. The major impact of this conclusion is that this implies that all the 
solutions are periodic. 

Consider the top, left diagram in Figure 1.4.4, where the point P is selected 
to be on the positive y-axis. The trajectory, passing through P, in the first- 
quadrant, must start with zero slope on the y-axis and then have a negative 
slope in the remainder of this quadrant. When it makes contact with the x- 
axis, it does so with a vertical, negative slope. Applying the transformation 
T2, reflection in the z-axis, gives the third diagram. Note that the point P’ is 
the reflection of P and thus lies below the origin the same distance as P lies 
above the origin. 

If the transform 7}, reflection in the y-axis, is applied to the third diagram, 
then we obtain the result of the fourth diagram, i.e., a simple closed curve. 
This result is very important because it implies that x(t) and y(t) are peri- 
odic solutions, with the possibility that the period may depend on the initial 
conditions. 

The discussions of this section are very general and are applicable to other 
dynamic systems. We will use them in other investigations to be given in this 
book. 
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FIGURE 1.4.3: Null-clines and signs of dy/dx in the four quadrants for the 
harmonic oscillator dynamic system. 


1.5 HYPERBOLIC SINE AND COSINE: DERIVED 
FROM SINE AND COSINE 


Using the Euler formula and its complex conjugate, see Equation (1.14), 
the cosine and sine have the representations 


i0 1 ,-i0 
cos 6 = os (1.64a) 
10-70 
sind = ——, (1.64b) 
i 


The functions obtained by replacing @ by 76, where 1 = /—1, are related to 
the so-called hyperbolic functions, defined as 


04,8 

hyperbolic cosine = — = cosh6, (1.65a) 
0 — e-9 

hyperbolic sine = = sinh 6. (1.65b) 


The exact relationship between the cosine/sine and the cosh/sinh is given by 
the relations 
cosh(#) = cos(i#), i sinh(6) = sin(#@). (1.66) 
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FIGURE 1.4.4: Geometric proof of simple, closed phase-space trajectories 
for the harmonic oscillator dynamic system. 
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10, Cosh @ 


-10 


FIGURE 1.5.1: The hyperbolic cosine and sine functions. 


Note that cosh @ and sinh @ are, respectively, even and odd. Further, unlike 
the cosine and sine, they are not periodic and become unbounded as 0 +> 
(+)oo. See Figure 1.5.1. 

The cosh/sinh satisfy many relations which are similar to those for the 
cosine and sine function; see the Appendix for a listing of the more important 
formulas. However, one of the most important of these will now be derived. 

Observe that (cosh @)? and (sinh @)? are 


0 -~0\ 2 20 | 20 
(cosh 8)? = (=) ie se eae (1.67a) 
2 4 
0 .-0\2 20, 4-20 _ 
(sinh 0)? = (=) = ote? (1.67b) 


and, if we subtract the second expression from the first, then the following 


relation results 
(cosh 0)? — (sinh 6)? = 1, (1.68) 
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which is to be contrasted with (cos)? + (sin @)? = 1 for the cosine and sine 
functions. 
Returning to Equation (1.67), we have 


1 20 —20 1 
(cosh)? = 5 + ae a (5) [1 + cosh(26)], (1.69a) 
1 20 —20 1 
(sinh)? = — (5) ca calle (5) [cosh(20) — 1], (1.69b) 
2 4 2 
and the interesting relation 
(cosh 0)? + (sinh @)? = cosh(26). (1.69c) 


All of the hyperbolic formulae, given in the Appendix, can be derived using 
similar techniques. 
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SYSTEM DERIVATION 


Consider a system having the Hamiltonian 


1 1 
He) = (5) W-2) = 5. (1.70) 
2 2 
The equations of motion are 
dx _ OH(x,y) dy OH (x,y) 
dt Oy yd Ox me el 
with the initial conditions 
x(0)=0, y(0) =1, (1.72a) 
or 
x(0)=0, y(0)=-1. (1.72b) 


(For some insight into these particular selections, see Figure 1.6.1.) Therefore 


dx dy 
Ss = i, 
dt Y, wv, (1.73) 


and both x and y satisfy the second-order differential equation 


Oe ag (1.74) 
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FIGURE 1.6.1: Plot of y? — 2? = 1. 


whose general solution is 
w(t) = Aet* + Be. (1.75) 
For the initial conditions given in Equation (1.72), we have the solutions 


x(0) =0, y(0) =1: a(t) =sinh(t), y(t) = cosh(t), (1.76a) 
x(0) =0, y(0) = —1: a(t) = sinh(t), y(t) = — cosh(t). (1.76b) 


1.7 0-PERIODIC HYPERBOLIC FUNCTIONS 


We now consider the nonstandard case of constructing periodic hyperbolic- 
type solutions for the equation 


y—x? =1, (1.77) 


Figure 1.7.1 defines the variables r and 6 associated with a point (x, y) 
located on the curve represented by Equation (1.77). 
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FIGURE 1.7.1: Definition of the variables r and 6. 


We define three periodic hyperbolic functions as follows: 


pch(0) = x, (1.78a) 
psh(@) =y, (1.78b) 
pdh(@) =r (1.78c) 


The expressions on the left side are, respectively, called the periodic hyper- 
bolic cosine, the periodic hyperbolic sine, and the periodic hyperbolic dine 
functions. 

There are two functional relations connecting x, y and r; they are 


Dian AD ag 
ee Ty (1.79) 
yo =a = 1. 
Note, we can also express x and y as 
x(0) =r(0)cosé, y(0) =r(6)sin@, (1.80) 


where r(@) can be determined by substituting these representations into the 
second of Equation (1.79), i.e., 


r(0)?[(sin 0)? — (cos 0)?] = 1. (1.81) 


Therefore, 
1 
r(0) = ———_., (1.82) 
(sin 6)? — (cos 6)? 
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but, restrictions must be placed on this result. Since r(6) is the distance from 
the origin to a point on the curve y? — x? = 1, the following angular intervals 
are not allowed 

On <4: 

Src 9 < 2, (1.83) 

a <0 < 2n. 


The reason for this is that for these angular intervals, there is no curve to be 
intercepted by a ray from the origin and r(@) is not defined. 

An alternative way to obtain the non-allowed angular regions is to examine 
the expression in the radical for Equation (1.82), ice., 


(sin 0)? — (cos 0)? = (5) (1 — cos 20) — (5) (1 + cos 20) 
= —cos 26 (1.84) 


The requirement that 
—cos26 > 0, (1.85) 


gives the results listed in Equation (1.83). 
The explicit analytical expressions can now be given for the periodic hy- 
perbolic functions; they are 


1 
r(@) = pdnh(@) = Vein)? = (con)? (1.86a) 
cos 6 
(0) = peh(0) = Tandy? = (en (1.86b) 
y() = psh(8) = ae (1.86¢) 


(sin 0)? — (cos 0)?’ 


where the angular intervals for which they are defined are listed in Equation 
(1.83). 

Finally, it can be easily seen that the three functions become unbounded at 
the boundaries of the angular intervals where they are defined; see Figure 1.7.2, 
where r(@), (0), and y(@) are shown over the interval, 0 < 0 < 2r. 


1.8 DISCUSSION 


Let us examine what we have found in the work given in this chapter. Of 
importance is the fact that the equation 


e+y=1, (1.87) 
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FIGURE 1.7.2: Plots of r(@) and «(@) from Equations (1.86). 
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can be interpreted in several different ways. First, from the perspective of 
plane geometry, this equation represents a circle. Second, this equation can 
be considered (proportional to) the Hamiltonian for the harmonic oscillator 
dynamic system, which gives a model for a broad range of phenomena in 
the natural sciences. Third, Equation (1.87) can be taken as a functional 
equation and, with this interpretation, it allows for the construction of many 
different classes of solutions, periodic and non-periodic, and, continuous and 
non-continuous. Chapter 5 gives a discussion of certain classes of periodic 
solutions. 
Similar results hold for the equation 


y’—a? =1, (1.88) 


for which unbounded periodic and non-periodic solutions exist. 
In the chapters to come, we apply the general methodology of this chapter 
to investigate solutions to 
ly|? + ||? = 1, (1.89) 


where 
p>od, q>O0. (1.90) 


NOTES AND REFERENCES 


Sections 1.2, 1.3, and 1.5: The elementary properties of the trigono- 
metric cosine and sine, and hyperbolic cosine and sine functions can be found 
in any standard textbook covering elementary geometry and trigonometry. 

Section 1.4: A brief introduction to Hamiltonian dynamics is given in the 
Appendix, Section E. 

Sections 1.6 and 1.7: As far as I am aware, the results given in these 
sections is new, and generalize the discussions in Sections 1.2, 1.3, and 1.4. 


Chapter 2 


ELLIPTIC FUNCTIONS 


2.1 INTRODUCTION 


The main purpose of this chapter is to introduce the standard Jacobi 
elliptic functions and show that other related elliptic functions also exist. 
Further, we demonstrate that some such functions are periodic, while others 
are non-periodic. 

Section 2.2 examines the construction of what we name 6-periodic elliptic 
functions. These functions are based on the geometric properties of the ellipse 


2 a 
Ytaal, a> (2.1) 


In Section 2.3, we investigate the case where the expression in Equation 
(2.1) is taken to be the Hamiltonian of a one-space-dimensional dynamic sys- 
tem. 

The standard Jacobi elliptic functions are investigated in some detail in 
Section 2.4. We also provide a listing of many of their important mathematical 
properties in Section 2.5. 

Section 2.6 presents a discussion of the Jacobi elliptic functions as the 
solutions to a three-dimensional dynamic system, for which the equations of 
motion are three coupled, nonlinear, first-order, ordinary differential equa- 
tions. 

The next two Sections, 2.7 and 2.8, consider “hyperbolic” generalizations 
of the work given in Sections 2.2 and 2.3. 

Finally, in Section 2.9, we give a general summary of what was done and 
its significance. 


2.2 6-PERIODIC ELLIPTIC FUNCTIONS 


Consider the graph of the elliptic curve given in Figure 2.2.1. Let P(x, y) 
be a point on this curve and let r denote a ray from the origin to P(2,y), 
where @ is the angle between the ray and the positive z-axis. The following 
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(0,1) 


(0, =) 


FIGURE 2.2.1: Graph of y2+ 2; = 1, and the definitions of r and 
denotes a point on the ellipse, and a > 1. 


two relations hold 


2 x 
+> ’ 

y a 

ra=atty’, 


where x(0), y(@), and d(@) are 


Substituting Equations (2.3a) and (2.3b) into Equation (2.2a) gives 


2 


1 =r? sin(0)? + (5) (cos 6)? 


=r? [1 — (cos6)?] + (5) (cos 0)? 
=r? [1 —k?(cos6)"] , 


6. P(x, y) 


(2.4) 


ELLIPTIC FUNCTIONS 


where k? is defined to be 
Pet 
ae 


Since a > 1, we have the bounds 
Ok ie oe Sa: 


Note that a = 1 reduces the ellipse to a circle. 
Solving Equation (2.4) for r(@) gives 
1 


sO 1 — k?(cos 0)?’ 


and it follows that 
1l<r(0)<a, 0<0< 2z, 


where we have suppressed, for the time being, the dependence on k?. 
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(2.5) 


(2.6) 


(2.7) 


(2.8) 


The elliptic dine, cosine, and sine periodic functions are defined by the 


expressions 


edn(0, k2) : 


r(@) = 


0 
ecn(0, k?) = r(0) cos@ = = 


esn(0,k”) = r(0)sin@ = ee 


These three functions have the properties: 
1) bounds 


1 < edn(0,k?) <a 
—a < ecn(0,k?) <a 
—1<esn(0,k?) <1 


2) even/odd-ness 


edn(—0, k”) = edn(0, k”) 
ecn(—6, k?) = ecn(6, k?) 
esn(—0, k”) = (—)esn(0, k?) 


3) periodicity 
edn(@ + 7, k?) = edn(6, k?) 
ecn(@ + 2r, k?) = ecn(6, k”) 
esn(6 + 27, k?) = esn(0, k?) 


1 — k2(cos@)?" 
1 — k?(cos 6)?’ 


1 — k?(cos 0)? 


(2.9) 
(2.10) 


(2.11) 


(2.12) 


(2.13) 


(2.14) 
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4) special values 


2\ — x 2 = 2 = 
edn(0,k2) =a, edn(%,k?)=1, edn(x,k*) =a (2.15a) 
edn (32,k?) =1, edn(27,k?) =a 
ecn(0, k?) = . ecn (3, k?) = 0, ae k?) = —a (2.15b) 
con (38,12) = 0, eon(2m, 4?) = a 
esn(0,k?)=0, esn(Z,k?)=1, esn(7,k?) =0 (2.15¢) 
esn (22,k?) =—-1, esn(2n,k?) =0. 


Figure 2.2.2 provides graphs of these three functions. 


2.3 ELLIPTIC HAMILTONIAN DYNAMICS 


Consider a dynamic system having the Hamiltonian 


H(z,y) = (5) (+5) =>. (2.16) 


(Compare this to the expression in Equation (2.1).) The equations of motion 
are 


dx  OH(zx,y) 
ee Ne 2.17 
dt Oy Y; ( a) 
dy OH (x,y) 1 
2 if 2.1 
dt Ox a} ettP 
and the initial values are selected to be 
x(0)=a, y(0)=0. (2.18) 
Equations (2.17) implies that x(t) satisfies the second-order, linear, differential 
equation 
dx 1 
pala a= = 2.19 
cs +()e 0 (2.19) 
with initial conditions iO 
«(0) =a, a io, (2.20) 


Therefore, x(t) is 


a@ eres (=) aes (=) (2.21) 
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FIGURE 2.2.2: Graphs of edn(6), ecn(9), and esn(0) for a = % or k? = 2. 
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and imposing the initial conditions gives the periodic solutions 


Mi dees (=) , A =sin (=) . (2.22) 


Physically, x(t) and y(t) correspond to the position and velocity of a linear 
harmonic oscillator having a period 


T = 2na. (2.23) 


2.4 JACOBI, CN, SN, AND DN FUNCTIONS 


We will now construct and examine the properties of the so-called Jacobi 
elliptic functions. Refer to Figure 2.2.1 for the definitions of r, 8, and a > 1. 
The eccentricity of the ellipse is defined to be 


dl 
e=k=y1——, a>. (2.24) 


Note that for a = 1 the ellipse becomes a circle. 
The three Jacobi elliptic functions are defined in the following way, in 
terms of (x, y,1, a): 


en= =, sn=y, dn=_>0. (2.25) 
a a 


Examination of Figure 2.2.1 shows that all three of these functions are periodic 
in 6, with period 27. However, this is not the independent variable selected 
for these functions. The standard independent variable, wu, is selected to be 


0 
u(e) = f r()dy, (2.26) 


where r(@) is the length of a ray from the origin to the ellipse, where @ is 
the angle that the ray makes with the positive z-axis. While not explicitly 
indicated, u depends also on k?, as defined by Equation (2.5). 

This particular definition of u(@) is based on construction of the Jacobi 
elliptic functions in terms of integral representations. See the note to this 
section for pertinent references on this issue. 

One consequence of this definition is that when e =k =0O0,ie.,a=r=1, 
u(@) is 


6 
u(0) =| (1)dw = 8, (2.27) 


and 
cn=cos#, sn=siné, dn=1. (2.28) 
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Thus, as expected, the elliptic functions reduce to trigonometric functions and 
the above definition of u(@) is consistent with the a priori known geometrical 
features of the problem. 

From the fact that 


r(—6) = r(6), (2.29) 
it is an easy calculation to show that 
u(—0) = —u(8). (2.30) 


Using the results in Equations (2.25), (2.26), and (2.30) allows the following 
results to be stated: 


(i) there are three elliptic functions and they are denoted (cn, sn, dn); 


(ii) these three elliptic functions depend on the independent variable, u, and 
the parameter k, i.e., 


cn=cn(u,k), sn=sn(u,k), dn = dn(u,k); 
(iii) the parameter k, called the modulus or the eccentricity, €, satisfies the 
bounds, 0< k <1; 


(iv) when a > 1, or k = 0, the elliptic functions reduce to the trigonometric 
functions; 


(v) the three elliptic functions are periodic in 0, with period 27; however, 
they are periodic in u with period T (which for the time being is un- 
known, but will be calculated later). 


The three elliptic functions, as we have constructed them, are called the Ja- 
cobi, “cosine,” “sine,” and “dine” functions. 


2.4.1 Elementary Properties of Jacobi Elliptic Functions 


From the definitions of the Jacobi elliptic functions, see Equations (2.25), 
(2.29), and (2.30), we conclude that cn and dn are even functions of u, while 
sn is an odd function of u, i.e., 


cn(—u) =cn(u), sn(—u) =—sn(u), dn(—u) = dn(u). (2.31) 
Using the definitions given by Equation (2.25) in the relations presented 
in Equations (2.2) gives 
sn(u)? + en(u)? = 1, (2.32a) 
dn(u)? + k?sn(u)? = 1. (2.32b) 


The Jacobi elliptic functions satisfy the bounds 


-1l<cn(u) <1, -1<sn(u) <1, 


ale 
IA 
Q 
3 
— 
~ 
wa 
x 
= 
— 
a 
Ww 
w 
wm 
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2.4.2 First Derivatives 


From Figure 2.2.1, we have 
6 = tan! (4) , (2.34) 
x 
and from this it follows that 
I 

dé = (=) (x dy — ydz), (2.35) 

and from taking the derivative of Equation (2.27) 

1 

du = rd@ = (z) (a dy — ydz). (2.36) 


Likewise, from Equation (2.1) 


d 
ydy+ = = 0, (2.37) 
a 
and from this, we obtain 
x a’y 
dy=—|-——]dx or dx=—|— } dy. (2.38) 
a*y x 


Substituting, respectively, the results in Equations (2.38) into (2.36), and sim- 
plifying gives 


=—ry, (2.39) 


dx _ ry 
du | (2)? 4 y2 


=5 (2.40) 
or F . 
a Sai sn, ost = dn en. (2.41) 
du du 
Now 
d ay _ 2,2 
$9) = Fe +) 
edge acme | (2.42) 
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i -()#+O# e 


Using the definitions from Equations (2.25), this equation can be written as 


and 


dd 
ce = ken sn. (2.44) 
du 
Observe that if we take the limit a + 1 or k? — 0, then these expressions 
for the derivatives become 
dcos@ _ d sin 6 dr(@) 


= —sij = ee 2.4 
sin 0, 7 cos 0, 76 0, (2.45) 


i.e., they reduce to the standard results for the usual trigonometric functions. 


2.4.3 Differential Equations 


All three Jacobi elliptic functions satisfy nonlinear, second-order, ordinary 
differential equations. We now derive these equations for (sn,cn,dn), where 
it is understood that sn = sn(u, k?), etc. 

Using Equations (2.32a),b, we have 


d 
ade 


du 


= (V1— sr?) (V1 —ksn®) (2.46) 


and replacing sn(u) by y(u), this equation becomes 


(zy = (1—9?)(1 — Fy): (2.47) 


Taking the derivative and simplifying the resulting expression gives 


d’y 2 2,3 
—5+(14+k*)y —2k-y’ =0, (2.48) 
du? 

which is the nonlinear, second-order, differential equation for which with the 


initial values 
dy(0) 


0) =0 
y(0)=0, — 
gives the solution y(w) = sn(u). 

Let z(u) = acn(u), then substituting this into Equation (2.41) and squar- 
ing, gives 


= dn(0)en(0) = 1, (2.49) 


(=) = [(1 — k*) + k*2”] (1-27). (2.50) 
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Taking the derivative of the expression, canceling common factors, and sim- 
plifying gives the following nonlinear, second-order differential equation for 
x(u) 
ao | (1 — 28?)n +282? = 0 2.51 
ae +(1- ya + xz’ =0, (2.51) 


for the initial conditions 


=0. (2.52) 


2 
ae — (2—k?)w + 2w* =0, (2.53) 
du 
with Anh 
“ori, ay (2.54) 
du 


In summary, if V(u) is any of the three Jacobi elliptic functions, 
(cn, sn,dn), then V(u) satisfies first- and second-order derivative relations 
which take the form 


d 2 
(=) =A+BV*?+CV*, (2.55a) 
U 
d2 
ad +aV + BV? =0, (2.55b) 


where (A, B,C) are functions of k?, and (a, ) are determined by (A, B,C), 
and therefore are functions of k?. 


2.4.4 Calculation of u(#) and the Period for cn, sn, dn 


The calculations and analysis given above makes it clear that u, as defined 
by Equation (2.26), is the appropriate independent variable for the Jacobi 
elliptic functions. We now calculate u(@). 

From Figure 2.2.1, we have « = rcos@ and y = rsin@. Therefore from 


ae 
rP=ety, yt+sa=l, (2.56) 


it follows that 
(0) = — (2.57) 


and therefore 


ay db 
u(O) = af Vn Dena (2.58) 


Using the “look-up” method, i.e., seeing if the integral appears in one of the 
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standard reference books on mathematical tables and formulas, we find this 
integral in Gradshteyn and Ryzhik (see page 173, formula 2.597.1), and obtain 
the result 


asin 6 


u(0) = F 2 sin-? | —_“" —_| 
1+ (4) (sin 0)? 


where F'(¢,k) is an elliptic integral of the first kind. Remember that a and k 
are related; see Equation (2.24). 
If 0 = 7/2, then 


(2.59) 


u (=) =F (=. k) = K(k), (2.60) 


where K(k) is the complete elliptic integral of the first kind. Both F(¢,k) 
and K(k) are tabulated for values of ¢ and k, respectively, in the intervals, 
0<¢<7/2and0<k<1. 

Since, in terms of the variable 0, the period of the Jacobi elliptic functions 
is 27, and the corresponding period in the variable u is 

T =4u (5) =4F (5. k) = 4K(k). (2.61) 

Thus, if ef(u) denotes any of the three Jacobi elliptic functions, (cn, sn, dn), 
then the following periodicity condition holds 


ef(u+4K(k)) = ef(u). (2.62) 


2.4.5 Special Values of Jacobi Elliptic Functions 


In addition to the upper and lower bounds, expressed in Equation (2.33), 
these functions are zero or take their respective maximum/minimum magni- 
tudes at the following indicated multiples of K = K(k): 


sn(0) = cn(0) =1 dn(0) =1 
sn(k) =1 cn(i) =0 dn(K) = + 
sn(2K) =0 cn(2K) =—-1 dn(2K) =1 
sn(3K)=—-1  en(3K)=0 dn(3K) = + 
sn(4k) = cn(4k) =1 dn(4k) =1 


Also, note in summary: 

i) When k = 0, then K(0) = 7/2, and the three Jacobi elliptic functions 
reduce to (cos 6, sin 6, 1). 

ii) The fundamental period of cn and sn is 4K(k), while that of dn is 
2K (k). 

iii) sn has zeros at even multiples of AK, while cn has zeroes at odd mul- 
tiples of K. 

iv) sn has maximum magnitudes at odd multiples of A’, while cn has 
its maximum magnitude at even multiples of K. For each, the maximum 
magnitudes are one. 
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(v) dn is positive for k values in the interval, 0 << k <1. 
(vi) Since 1 < r(@) <a, then 


du 
a r(8) > 0, (2.63a) 
0<u(0)<a0, a>. (2.63b) 


Consequently, u(@) is a monotonic increasing function of 6, bounded by the 
two functions ui(0) = 6 and u2(@) = ad. 
See Figure 2.4.1 for graphs of sn(u), cn(u), and dn(u). 


2.5 ADDITIONAL PROPERTIES OF JACOBI 
ELLIPTIC FUNCTIONS 


2.5.1 Fundamental Relations for Square of Functions 


cn? + sn? = 1 (2.64) 
cn? + (1 — k*)sn? = 1 (2.65) 
dn? + k?sn? =1 (2.66) 


2.5.2 Addition Theorems 


ety) = 1 — k?sn?(x)sn?(y) co, 
sn(xt+y) = ee (2.68) 
iene dn(x)dn(y) — k?sn(x)sn(y)en(x)en(y) (2.69) 


1 — k?sn?(x)sn?(y 
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FIGURE 2.4.1: Graphs of the three Jacobi elliptic functions for k = 0.8. 
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2.5.3 Product Relations 


sn?u — sn?v 
sn(u+v)-sn(u—v) = hes anti 
cn?u — sn?udn7u 


A NM OY re oir 
dn?u — k?cn?u sn2u 
saa ae a pe en re 


snucnudnv + snucnvudnu 


1— k?sn2u sn2v 
snudnucnv+ snudnvcnu 
1— k?sn2usn2v 


cnudnucnv + k?snusnv 


1— k?sn?2u sn2v 


2.5.4 cn, sn, dn for Special k Values 


sn(u,0) = sin u sn(u, 1) = tanhu 
cn(u, 0) = cos u cn(u, 1) = sechu 
dn(u,0) =1 dn(u, 1) = sechu 


2.5.5 Fourier Series 


Let 


oo 


— mu 
U= 3K? 


(2.70) 
(2.71) 
(3:78) 
(2.73) 
(2.74) 


(2.75) 


(2.76) 
(2.77) 
(2.78) 


(2.79) 


where K = K(k?) and K’ = K(1—k?), then cn(u), sn(u), and dn(u) have 


the expansions 


T = n+3 
cn(u) = (4 :) 55, () cos[(2n + 1)v] 


n=0 


dn(u) = (==) “fs (=) y (<,) cos(2nv) 


E— K(k‘)? Qn? \ = nq” 
cn? = AS + Ssoany Cos 2nv 
K - k? FORA a Lesgee nt) 


n= 


(2.80) 


(2.81) 


(2.82) 


(2.83) 
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sn? = (=) = (Sa) ys (; w.) cos(2nv) (2.84) 
dn? = (=) 2 (3) 3 (, <.) cos(2nv). (2.85) 


n=1 


In the above, K and E are complete elliptic integrals of the first and second 
kind. 


2.66 DYNAMICAL SYSTEM INTERPRETATION OF 
ELLIPTIC JACOBI FUNCTIONS 


An alternative derivation and interpretation of the three Jacobi ellip- 
tic functions, (cn,sn,dn), is to consider them as the solution to a three- 
dimensional dynamical system. This means that they are represented as the 
solutions to three coupled, nonlinear, first-order differential equations. This 
methodology was summarized in a paper by Meyer who based his work on the 
earlier results of Hille and Chicone. This section provides a concise summary 
of Meyer’s presentation, and for the most part, his notation will be used. 

Anyone desiring the full proofs and other analytical background details 
can consult the following references: 


1. E. Hille, Lectures on Ordinary Differential Equations (Addison-Wesley, 
Reading, MA, 1969). 


2. C. Chicone, Ordinary Differential Equations with Applications (Dover, 
New York, 1961). 


3. K. R. Meyer, Jacobi elliptic functions from a dynamical systems point 


of view, American Mathematical Monthly 108 (2001), 729-737. 


2.6.1 Definition of the Dynamic System 


Let k be a real number such that 


0<k<1. (2.86) 
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The dynamic system is defined by the three coupled, nonlinear differential 


equations 
dx 
= yz 
dt YZ, 
dy 
— = -22 
dt : 
dz 
— = _fk? 
at LY, 


where the initial conditions are taken to be 


x(0)=0, y(0)=1, 2(0)=1. 


2.6.2 Limits k > 0* and k > 17 


We have 
x(t) sint 
Lim (t)}] = | cost], 
k-0 (t) 1 
and 


2.6.3 First Integrals 
The functions I(x, y) and J(x, y), defined as 
T(e,y) =a? +47, 
J(a;z) = kx? + 27, 


are first integrals of Equations (2.87), i.e., let 


x(t) or(t) 
y(t) | = | o2(t) 
z(t) os(t) 


be any solution of Equation (2.87), then 


I(¢1(t), 62(t)) = constant, 
J(¢1(t), 63(¢)) = constant. 


Proof: 
d dx dy 
—TI(x(t t)) = 2a— + 2y— 
i (elt), y(t) = 2a + Qu 
= 2e(y2) + 2y(—22) 
= 0, 


(2.87a) 
(2.87b) 


(2.87c) 


(2.88) 


(2.89) 


(2.90) 


(2.93) 


(2.94a) 
(2.94b) 


(2.95) 


and 
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d g dx dz 
eet = — 427 
a (tte), (4) 2k t + an 
= 2k? x(yz) + 22(—k?2y) 
=0. (2.96) 


2.6.4 Bounds and Symmetries 


The following statements are true: 


(ii) 


(iii) 


(vi) 


The functions, (x(t), y(t), z(t)), are periodic in t. 
For fixed k, 0< k <1, then for allt > 0 


a(t) +y7(t) =1, 
fee +2(t) = 1. oe 


The following inequalities are satisfied 


If (x(t), y(t), z(t)) is a solution of Equations (2.87), then so are 


( 
(2(-t), —y(-t), 2(-t)) (2.99) 
)) 


For fixed k, 0 < k <1, we have 
x(-t)=-2(t), y(-t)=y(), (8) = 200). (2.100) 
There exists a constant, K(k), such that 


a(t +4K) = a(t), 
y(t + 4K) = y(t), (2.101) 
a(t ree 2K) = y(t), 


ie., z(t) and y(t) are periodic, with period 4K, while z(t) is periodic, 
with period 2K. 
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2.6.5 Second-Order Differential Equations 


The following nonlinear, second-order, differential equations are satisfied 
by (x(¢), y(t), 2(¢)) 


dx 2 2.3 
{-4+(1+k*)x -— 2k’ =0, 
ge INE (2.102) 
x(0) = 0, ae i 
2 
TH + (1 — 2k?)y + 2k?y? = 0, 2-103 
_ dy(O) __ (2. ) 
y(0) a 1, dt - 0, 
dz 2 3.2 
—(2-—k 22° = 0 
ae — ( Lae nee (2.104) 


2.6.6 Discussion 


Comparison of the results derived by Meyer, as summarized in this section, 
with the results given in Section 2.4, shows that the correspondence between 
(sn, cn, dn) and (a, y, z) is 


snor, ney, dnoz. (2.105) 


It is this correspondence which allows Meyer to state that the Jacobi ellip- 
tic functions can be interpreted as solutions to a three dimensional dynamic 
system. 

A major feature of Meyer’s article is that he demonstrates that the Ja- 
cobi elliptic functions can be applied to several important problems arising in 
mathematics, physics, and engineering; these include: 


e the pendulum system, 
e evaluation of elliptic integrals, and 


e physical systems having quadratic or cubic forces. 


2.7 HYPERBOLIC ELLIPTIC FUNCTIONS AS A 
DYNAMIC SYSTEM 


Consider a dynamic system modeled by the Hamiltonian 


H(2,y) = (5) ? = (=)' - 7 (2.106) 
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The equations of motion are 


< = =y, (2.107a) 
“ = = = (=) aL, (2.107b) 
with the initial conditions 
x(0)=0, y(0)=1, (2.108) 
or 
x(0)=0, y(0)=-1. (2.109) 
First, observe that both x(t) and y(t) satisfy the second-order differential 
equation : 
= (=) o (2.110) 
for which the general solution is 
w(t) = Ae/® 4 Be“(/%), (2.111) 


where A and B are arbitrary integration constants. 

Second, from Equation (2.111), the solutions corresponding to the two 
sets of initial conditions can be calculated. These solutions are given by the 
expressions 


), (419) 


x(t) = —asinh (4) ‘ 
y(t) = —cosh (4). 


Figures 2.7.1 and 2.7.2 present the graphs of these two solutions. 

Note that these solutions are the corresponding elliptic hyperbolic func- 
tions for the case where Equation (2.106) is considered the Hamiltonian of a 
one (space)-dimensional system. 


x(0)=0, y(0)=-1: (2.113) 


2.8 HYPERBOLIC 6-PERIODIC ELLIPTIC 
FUNCTIONS 


Consider the curve , 
yo (=) =1, (2.114) 
a 
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go ett) 


-2 -1 0 1 2 t 


FIGURE 2.7.1: Graphs of Equations (2.112). 
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FIGURE 2.7.2: Graphs of Equations (2.113). 
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FIGURE 2.8.1: Graphs of y? — ea = 1, and definition of r, 0, and P(z, y). 


which forms the basis for the treatment of the hyperbolic 6-periodic elliptic 
functions. 

The variable r is the segment of the ray from the origin (0,0), to a point 
P(a,y) on the curve given by Equation (2.114). Using 


Dees er 
Pe ee (2.115) 
xz=rcosé, y=rsin0, 
then r(@) can be calculated and its value is 
1 
(0) = ————.. (2.116) 


1— (4) (cos6)? 


Since r(@) is only meaningful for r(@) real and positive, this implies that the 
relevant @ values are those for which 


1+? 
1- ( - ) (cos 0)? > 0 (2.117) 


2 
|cos 6| < View (2.118) 


A(a) = cos! ( a ; (2.119) 


or 


If A(a) is taken to be 
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then the allowed angles are those for which 


(5 - A) <0< (5 + A) ; (2.120a) 
(= - A) <0< (= + a) ; (2.120b) 


We define the hyperbolic 6-periodic elliptic functions using the following 
notation 


6 — deh = r(8), 
6 — ceh = r(6) cos9, (2.121) 
6 — seh = r(0) sin 9, 


where (0) is defined. The graphs of these functions are similar to those shown 
in Figure 1.7.2. 


2.9 DISCUSSION 


This chapter, just as in Chapter 1, allowed us to explore what it means to 
have periodic solutions to a two-term quadratic equation in x and y. It turns 
out that both standard and non-standard types of both periodic and non- 
periodic solutions exist. Further, the periodic solutions may not be bounded. 
Also, of great interest and great value, is that often explicit solutions can be 
calculated. Similarly, we found that considering these “systems” as Hamilto- 
nian also provided additional periodic- and hyperbolic-type functions. 


NOTES AND REFERENCES 


Sections 2.3, 2.4, and 2.5: The discussions and results in these sections 
are based on standard presentations such as appear in these books: 


(1) E. H. Neville, Jacobian Elliptic Functions (Oxford University Press, Ox- 
ford, 1944). 


(2) P. Byrd and M. D. Friedman, Handbook of Elliptic Integrals for Engi- 
neers and Physicists (Springer-Verlag, Berlin, 1954). 


For completeness, we give below the definitions of, respectively, the complete 
elliptical integrals of the first and second kind: 


m/2 do 
T 
F k =f k,= = ———— 
(x) ( 5) | 1 — k?(sin 6)2 
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where r 
do 
4) ——*—__.. 
Fs?) [ 1— k2(sind)2 
and je 
E(k) =E (x, =) = | 1 — k2(sin 0)? do 
where 


E(k, ¢) = a V1 — F2(sin 6)? 6. 


An alternative way to construct the Jacobi elliptic functions is to start 
with the relation 


0 
d 
os 
0 1 —k(sindy 


and define (sn, cn, dn) as follows 
sn(u,k) =sin0@, cn(u,k) = cosd 


dn(u, k) = V/1— k?(sin 6)?. 
Section 2.6: As previously noted, this section is essentially a summary of 
some of the results reported in the article 


1) K. R. Meyer, Jacobi elliptic functions from a dynamical systems point 
of view, American Mathematical Monthly 108 (2001), 729-737. 


See also, 


2) K. R. Meyer and G. R. Hall, Introduction to Dynamical Systems and the 
N-body Problem (Springer-Verlag, New York, 1992). 


Chapter 3 
SQUARE FUNCTIONS 


3.1 INTRODUCTION 


The square periodic functions are derived from the geometry of the curve 
ly| + |a] = 1. (3.1) 


This curve corresponds to a square in the z-y plane; see Figure 3.1.1(a). 
Likewise, in Figure 3.1.1(b), the variables r and @ are defined. In particular, 
if P(x, y) is a point on the boundary of the square, then r is the segment of a 
ray from the origin to P(a, y) and @ is the angle that r makes with the positive 
x-axis. Thus, it follows that, in addition to Equation (3.1), we also have 


rae ty’. (3.2) 


We now define three periodic square trigonometric functions associated 
with the geometry of Figure 3.1.1 by the expressions 


sqc(@) = «(0) = r(0) cos 9, (3.3a) 
sqs(0) = y(0) = r(6) sin 9, (3.3b) 
sqd(0) = r(6), (3.3¢) 
where, from its definition, 
r(@)>0, O<@0<2n. (3.4) 


The three functions sqc, sqs, and sqd, are called respectively, the square- 
cosine, square-sine, and square-dine functions. 


Note that 
sqc(@ + 27) = sqc(9), (3.5a) 
sqs(0 + 27) = sqs(0), (3.5b) 
sqd C + =) = sqd(@). (3.5c) 


Plots of these functions are given in Figure 3.1.2. 


47 


48 GENERALIZED TRIGONOMETRIC & HYPERBOLIC FUNCTIONS 


FIGURE 3.1.1: (a) Plot of |x| + |y| = 1. (b) Definition of r and 0. 


sqd 


sqs 
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1.0 


- 1.0 


-1.0 


FIGURE 3.1.2: Plots of sqd(@), sqc(@), and sqs(@) vs 0. 
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3.2. PROPERTIES OF THE SQUARE 
TRIGONOMETRIC FUNCTIONS 


Explicit, analytic expressions can be determined for sqc(@), sqs(0), and 
sqd(@), by beginning with Equation (3.1) and substituting « = rcos@ and 
y =rsin@; doing this gives 

ly| + |a| = |r sin 6| + |r cos 6| 
= 1r(0)[|sin 8] + |cos 4]] 


=, (3.6) 
Therefore r(@) is 
1 
¢) = —————_. ; 
ue) |cos 6| + |sin 0|’ et) 
and 
cos 6 
0) = ————_ . 
sqc(8) \cos 6| + |sin 6|’ ee) 
sin 0 
9) = ——_—____. 3.8b 
2a8(9) — a8] + sind] Gee) 
1 
d(é) = —————_. : 
saa\@) |cos 6| + |sin 6| is:88) 


From the definitions of the square trigonometric functions, in the x-y plane 
or from the explicit expressions given in Equations (3.8), the following prop- 
erties are easily obtained: 


(i) The three square trigonometric functions are periodic; see Equation 


(3.5). 
(ii) These functions also satisfy the nonlinear functional relations 
|sqc(0)| + |sqs(A)| = 1. (3.9) 
[sqc(@)]* + [sqs(9)]? = [sad(6)]?. (3.10) 


(iii) sqc(@) and sqd(@) are even functions, while sqs(@) is an odd function of 
0, i.e., 
fh = sqc(),  sqs(—8) = —sqs(0), on 


sqd(—0) = sqd(6) 
(iv) These functions satisfy the bounds 
—1 < sqce(6) < +1, 
—1< sqs(@) < +1, (3.12) 
< 


a < sqd(0) < +1. 
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(v) The zeros, and min/max values for sqc(@) and sqs(0) are at the same 6 
values as those for the trigonometric cos @ and sin @ functions. 


(vi) A direct, but somewhat long calculation gives, for the first derivatives, 
the expressions 


=le Os Ges 
dsqc(9) <1, - 220 <9 


= 2 2 
dé +1, << % [sqd(9)]°, (3.13a) 
+1, 3a a 0 < 27 
+l, “Osos 
dsqs(9) | —1, ReQen : 
do |-l, ~<0<% [sqd(0)]”, (3.13b) 
7 1, is < é < 27 


sqs(0) — sqc(6), 0<0<45 

dsqd() _ | —sqs(0) — sqc(0), 5 <O<n 
dé —sqs(0) + sqc(0), T<O0< an 
sqs(9) + sqc(0), 32 <0< 27 


sqd(0). (3.13¢) 


These derivatives do not exist at 
T 37 


= (3.134) 


3.3 PERIOD OF THE SQUARE TRIGONOMETRIC 
FUNCTIONS IN THE VARIABLE u(6) 


While we have used 6 as the independent variable for the square- 
trigonometric functions, Schwalm suggests that a more fundamental indepen- 
dent variable is u(@), defined 


0 
u(a) = f r(w)dw, (3.14) 


where r(@) is the expression given in Equation (3.7). 
Note that for the circular or standard trigonometric functions, where 


2 as y? = 1, 
= cosé 
ase sae (3.15) 
y = sind, 
r=din@d=1, 


we have 


0 
u(0) =| (1)deb = 0. (3.16) 
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Since, in general, we expect dine functions to be even functions of #, then 
u(@) will be an odd function of 6, ice., 


r(—0) = r(0) > u(—6) = —u(8). (3.17) 


Since the square trigonometric function sqd(@) = r(@), satisfies the upper 
and lower bounds given in Equation (3.12), it follows that the bounds for u(@) 
are 


6 
— <ul(9) < 6. 3.18 
5 (9) (3.18) 
We now calculate u(@) for the case where r(@) is given by Equation (3.7), 
Le., 
6 
dp 
0) = ——_——_—_—.. 3.19 
ule) | |cos 6| + |sin | eee) 
A somewhat more than elementary calculation gives the result 
O.n 
if A gg |e 8) 
uy) (+) - | tan($) |}? (3.20) 
Va0s 5 


The restrictions in the range of # are a consequence of the fact that |z|+]y| = 1, 
does not have continuous derivatives, in 6, at 6 = 0, 7/2, 7, 37/2. 
The period, T, in the variable u, is 


T = 4u (5) = 4,985 801922..., (3.21) 


where the factor 4 is a reflection of the “time” to go from 0 = 0 to 7/2, and 
is one-fourth the total “time” to go from 6 = 0 to 6 = 27, based on the 
symmetries reflected in the curve presented in Figure 3.1.1. Another way of 
stating this result is to observe that the values of u(@) for the intervals (3, T), 
ar, sn). and (=, 2m) can be calculated by applying the shift and reflection 
properties of the trigonometric functions, and the result of doing this gives 
Equation (3.21). 
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3.4 FOURIER SERIES OF THE SQUARE TRIGONO- 
METRIC FUNCTIONS 


Let us denote the Fourier series of x(0), y(@), and d(@) by the expressions 


x(0) = r(0) cos@ = e ax cos(k8), (3.22) 
k=1 
y(0) = r(0) sind = 3 b;, sin(ké), (3.23) 
k=1 
1 c oS 


Inspection of Equations (3.22) and (3.23) suggests that the way to proceed 
is to first obtain the Fourier series representation for d(@), and then use this 
result to determine the Fourier series for x(@) and y(@). 

From Figure 3.1.1, we conclude that d(@) is periodic, with period 7/2. 
Therefore, its Fourier coefficients, the cz, can be calculated from the expres- 


sions 
m/4 
Ce (=) / d(0) cos (=) dé 
Tp od pe 1/4 


= (=) (2) il “" 4(0) cos(4k6)d0 


TT 
/4 

(2) (oath an, 
Now for 0 <@ < 7/4, 

Ee ral ee eee (3.26) 

V2 V2 
and, on this interval, 

|cos 6| + |sin 6] = cos@ + sind > 0. (3.27) 


Therefore, to calculate the cz, we must evaluate the integral 


— (=) ROSEY oa (3.28) 


T cos@ + sind 
Using the trigonometric relation 


cos(6; — 62) = cos 6; cos 62 + sin 6; sin 02, (3.29) 
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and setting, 0; = 6 and 02 = 7/4, gives 
T 


cos (0 -_ “) — (=) (cos @ + sin 0), (3.30) 


a result based on cos(7/4) = sin(w/4) = 1/./2. Therefore, Equation (3.28) 
can be expressed as 


m/A 
Ck = (4) | cos(4ké) do. (3.31) 
T 9 Cos (0 — =) 
Making the linear transformation 
G+ y=0--, (3.32) 
gives 
4/2\ f° 4k(w t+ % 
Ck = (*) i cos[4k(y + 2)] - dw. (3.33) 
T eye cos w 
Also, 
a k 
cos [4k (v + *)| = cos(4ky + km) = (—1)" cos(4kw), (3.34) 
and thus allows Equation (3.33) to now be written as 
4/2 °  cos(4kw) 
k 
fe eas - dw. . 
‘i OO Nw, se a (3.35) 
If we replace w by (—@), ice., 
~=-4, (3.36) 
then ii 
= 3(4 
T 0 cos } 


So, what can be done to evaluate the integral in the last equation? A 
way to proceed is to construct the Fourier series for (cos@)~'. We now give 
arguments for the following expression: 


1 
= 2[cos @ — cos 30 + cos 56 — cos 76 + --- + (—1)* cos(2k + 1)0 +--+] 
= 25 °(-1)* cos(2k + 1)0. (3.38) 
k=0 
Proof: 
1=1 


= 1 + cos 20 — cos 20 — cos 46 + cos 40 + cos 66 — cos60 + --- 
= (1+ cos 20) — (cos 20 + cos 40) — (cos 40 + cos 60) + --- 
+ (—1)"[cos(2n@) + cos2(n+ 1)6]+---. (3.39) 
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Now 
cos(2n@) + cos 2(n + 1)0 = 2cos@cos(2n + 1)6, (3.40) 


and this relationship allows Equation (3.39) to be written as 


l= yw cos #) cos(2k + 1)0 


k=0 
= 2cos@ 5 “(—1)* cos(2k + 1)6. (3.41) 
k=0 


Finally, dividing by cos@ gives the result expressed in Equation (3.38). 
Substituting Equation (3.38) into Equation (3.37) gives 


cp = (—1)* (4) 2 (—1)” < cos(4k@) cos|(2m + 1)6]dé. (3.42) 
m=0 


After a direct, but tedious calculation, we find that the integral has the value 


ie cos(4k@) cos[(2m + 1)6]d6 
0 


= (5) (ape) [et 2m +9 (F) 
4 (5) (G5) sin [(4x ~2m—1) (5) (3.43) 


Note that since (2m +1) is odd, then the denominators (4k + 2m + 1) cannot 
be zero. Therefore, we have for the Fourier coefficients, cz, the result 


ae. (22) .s = [(4k + 2m + 1) (4)] 


Te foot (4k + 2m + 1) 
sin |(4k — 2m — 1) (4)] 
r (4k — 2m — 1) \ ee) 


This means that the three square-trigonometric functions have Fourier 
series representation (see Equations (3.22), (3.23), and (3.24)), which can be 
calculated from the following relations: 


sqd(0) = oa aS. cr cos(k0), (3.45a) 
k=1 

sqc(0) = r(6) cos 8, (3.45b) 

sqs(0) = r(0) sin 0 (3.45c) 


We leave as an exercise, for the reader to complete, the explicit calculation of 
the Fourier coefficients for sqc(0) and sqs(@). 
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3.5 DYNAMIC SYSTEM INTERPRETATION OF |z|+ 
lyf =1 
Consider a dynamic system modeled by the following Hamiltonian 
A(a,y) = |e] + |yl = 1. (3.46) 
The corresponding equations of motion are given by 


dx  OH(x,y) dx OH(x,y) 


aie —= A 
dt Oy ° dt Ox” eet) 
or ; 
aan dy 
89"), Ge = ~sgn(z), (3.48) 
with the initial conditions 
x(0)=1, y(0)=0. (3.49) 
Note that in the z-y phase-space the “motion” takes the path 
(1,0) > (0,-1) > (—1,0) > (0,1) > (1,0). (3.50) 


The closed path indicates that this is a periodic solution. 

To determine the explicit solutions, as a function of t, we calculate the 
solution from (1,0) to (0,—1), ie., x(t) goes from x = 1, when t = 0, tox = 0, 
at time t = t1; and y(t) goes from y = 0, at t = 0, to y = —1, when ¢t = fy. 
Since in the x-y phase-plane, this motion is in the fourth quadrant, it follows, 
see Equation (3.48), that the equations of motion are 


dx dy 
ee Ok tap 3.51 
dt OG ) 
Integrating these equations gives 
a(t)=Ci—t, y(t)=C2-t, (3.52) 


where C; and C2 are integration constants. Imposing the conditions 


a(0)=1, y(0)=0, (3.53) 
gives 
Cy=1, C=), (3.54) 
and 
a(t)=1-t, y(t) =-t. (3.55) 


Now, if we required 
a(ti)=0, y(ti) =—1, (3.56) 


SQUARE FUNCTIONS 


then it follows that 
ty = 1. 


x(t) =1-t, 


Therefore, we have 


forO<t<l. 
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(3.57) 


(3.58) 


Continuing this procedure, the full solution can be determined and it is given 


by the following formulas: 


(1,0) — (0,-1): a % ford <t<L, 
(0,-1) > (-1,0): ee for le 9, 
(—1,0) + (0,1): faeces for 2<t<3. 
(0.1) -9 0,0): {0 OF for 3<t<4. 


See Figure 3.5.1 for plots of x(t) and y(¢) vs t. 


(3.59a) 


(3.59b) 


(3.59c) 


(3.594) 


From Equations (3.46) and (3.59), we find that x(t) and y(t) have the 


following properties: 
(i) They are bounded, i-e., 
—lsaQeth -lLey@ <1 


(ii) x(t) and y(t) are, respectively, even and odd, ie., 


(iii) z(t) and y(t) are periodic, with period T = 4. 
(iv) a(t) and y(t) satisfy the following basic integral relations 


(3.60) 


(3.61) 


(3.62a) 
(3.62b) 


(3.62c) 
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x(t) 


1.0 


0.5 


- 0.5 


- 1.0 


y(t) 
1.0 


0.5 


0.0 


- 0.5 


- 1.0 


FIGURE 3.5.1: Plots of x(t) and y(t) vs t. See Equations (3.59). 
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3.6 HYPERBOLIC SQUARE FUNCTIONS: DYNAM- 
ICS SYSTEM APPROACH 


The hyperbolic square functions are solutions to a dynamic system whose 
Hamiltonian is 


A(z,y) = |y| — |2| = 1. (3.63) 
A plot of this function appears in Figure 3.6.1. 
— de 0H dy oH 
4 y 
a ee SE 3.64 
dt Oy’ dt Ox’ eRe) 
we obtain the following equations of motion 
dx dy 
— = sgn(y), B= sgn(c). (3.65) 


Close inspection of the results given in Equations (3.65) shows that the flow 
(in terms of the evolution in time) takes place as in the directions indicated by 
the arrows in Figure 3.6.2. To summarize this situation, there are two classes 
of motions for the Hamiltonian given by Equation (3.63): 


(a) For #(0) = 0 and y(0) = —1, the dynamics evolves on the lower trajec- 
tory, with the general motion from right to left. 


(b) For «(0) = 0 and y(0) = +1, the dynamics takes place on the upper 
trajectory, with the general motion from left to right. 


The general, explicit solutions, in terms of t, for these two behaviors are, 
for —co <t< ow, 


; a(t) =t, 
upper trajectory : (3.66a) 
. =1+(¢, 
t)=-t 
lower trajectory : ay) , (3.66b) 
y(t) = —(1 + [é)). 


Figure 3.6.3 plots these functions with time. 
Finally, observe that x(t) is an odd function, while y(t) is an even function 
of t. 
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FIGURE 3.6.1: Plot of H(z,y) = |y| —|z| =1. 
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FIGURE 3.6.2: The flow along the trajectories in the z-y phase-space for 
H(a,y) = |y| — |a| = 1. See Equations (3.65). 
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0.5 


-0.5 


- 1.0 


t 
-4 -2 0 2 4 


FIGURE 3.6.3a: Plot of x(t) and y(t) for upper trajectory where (0) = 0 
and y(0) = 1. 
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0.5 


-0.5 


-10 


FIGURE 3.6.3b: Plot of x(t) and y(t) for lower trajectory where x(0) = 0 
and y(0) = —1. 
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FIGURE 3.7.1: Definition of r and @ for the periodic hyperbolic square 
functions. 


3.7 PERIODIC HYPERBOLIC SQUARE FUNCTIONS 


For this situation, we have 


ly| — |e] =1 (3.67a) 
yta=r’, (3.67b) 

with 
x(0) =r(0)cosé, y(0) =r(0)sin@. (3.68) 


See Figure 3.7.1 for the definitions of r and 6. 
If Equations (3.68) are substituted into Equation (3.67a), then the follow- 
ing result is obtained 


|r sin 8| — |r cos @| = r[|sin 6| — |cos|] = 1, (3.69) 
and solving for r gives 


1 
~ |sin 6] — |cos 6] 


r(6) 


Inspection of the expression for r(@) shows that the requirement 


(3.70) 


r(0@) > 0, (3.71) 
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permits only the following allowed angles 


T 37 57 Gis 
_— — —<0< —. (2 
rE. ri Fi (3.72) 


In other words, r(@) is not mathematically meaningful (for our purposes) for 
the angular intervals 


Tt 39 on 70 
< — —_—_, — : : 
Ne eae an ork 7 <9< 20 (3.73) 
We now define the fundamental three periodic functions, associated with 
the periodic hyperbolic square functions, as follows 


1 


hsqd(6) = r(¢) = ———_— .74 

phsqd(@) = r(0) = a (3.74) 
cos 0 

phsqc(0) — x(0) = r(0) cos 9 = Jsin | — ]cos 6)’ (3.75) 
in 0 

phsqs() = y(0) = r(0) sin@ = ——— (3.76) 


|sin 6| — |cos 6| 
Note that we have “named” these periodic functions, respectively, the periodic, 
hyperbolic, dine, cosine, and sine functions. The graphs of these functions are 
similar to those presented in Chapter 1, Section 1.7. 

Finally, where defined, we have 


r(—0)=r(0),  2(-@)=2(9),  y(—8) = —y(0), (3.77) 


and all are periodic with period 27. 
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Chapter 4 


PARABOLIC TRIGONOMETRIC 
FUNCTIONS 


4.1 INTRODUCTION 


The purpose of this chapter is to introduce the parabolic periodic- and 
hyperbolic-type trigonometric functions. A search of the literature indicates 
that little significant work has been done on such functions. In fact, much of 
what will be presented is new. The critical formula for our initial investigations 


is 
i\ a7 
= =, 4.1 
wit ($)2=5 (4.1) 


Inspection and analysis of this relation shows that its graph is a simple, closed 
curve; see Figure 4.1.1. Sections 4.2 and 4.3 are based on this equation. The 
first treats Equation (4.1) as a Hamiltonian and from it derives the equations 
of motion for x(t) and y(t), which can be explicitly solved. The second inves- 
tigates solutions from a geometric perspective in which x and y are functions 
of an angular variable @. In a similar manner, Sections 4.4 and 4.5 are based 


on the relation i i 
—~{io)7y2-2 4.2 
ly| (5) a= 5 (4.2) 


which gives rise to the corresponding hyperbolic-type solutions: (a(t), y(t)) 
and (2(4), y(9)). 


4.2 H(x,y) =|y|+(}) 2? AS A DYNAMIC SYSTEM 


Consider the Hamiltonian function 


(x,y) = |y|+ (5) w= 5. (4.3) 
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FIGURE 4.1.1: Plot of |y| + ($) 2? = $. 
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The equations of motion are 


dt OH (ay) 
or oy = sgn(y), (4.4a) 
dy OH (x,y) 
ign NE ae Sa 4.4b 
dt Ox ( ) 


Note that the point P[x(t),y(t)] travels around the perimeter of the curve 
given by Equation (4.3). If we select the initial conditions 


z(0)=1, y(0)=0, (4.5) 
then the path is 
(1,0) > (0 5) + (-1,0) 9 (0 5) , (1,0). (4.6) 


Explicit expressions can be calculated for (a, y) as functions of t. We now show 
how this is done. 
To begin, there are two general cases to consider, y > 0 and y < 0: 


dx dy _ 
Gots gee 
y>0:¢a(t)=tt+a, (4.7) 
2 
y(t) => (S) = ct + C2, 
dx d 
Book aa ao: 
y <0: 4 a(t) = -t+c3, (4.8) 


y(t) = a = c3t + C4, 


where the c’s are integration constants. 
Starting at (0) = 1 and y(0) = 0, the phase-space trajectory goes into 
the fourth quadrant. Using these initial conditions, it follows that 


ce =1, a =0, (4.9) 


and, for the fourth quadrant, we have the solutions 


a(t) =1-t, 
y(t) = 4 -t, (4.10) 
0<t<l 


The ¢ interval is gotten from the fact that there must exist a t = t,, such that 


a(ti)=0, y(ti) = 5, (4.11) 


and from the first two expressions in Equations (4.9), it follows that 


ty=1. (4.12) 


70 GENERALIZED TRIGONOMETRIC & HYPERBOLIC FUNCTIONS 
Now going to the third quadrant, we have 


a(t)=—-t+es, y(t) ="- = + 6, 
=i, 10) 


where the initial values come from inspection of Figure 4.1.1. We conclude 
that 


(4.13) 


C= 1; c= 0, (4.14) 
and 
x(t) =1-¢t 
y(t)= 5-1, (4.15) 
1<t<2. 
Using 
= ree 
a(4)=1, y(4) =0, 


we obtain the following expressions for x(t) and y(t), respectively, in the second 
and first quadrants: 


x(t) =t—3 
y(t) = — 5) SEA, (4.17) 
2<t <8, 
x(t) =t—3 
y(t) =— (5) 438-4, (4.18) 
3<t<4. 

The graphs of x(t) and y(t) are displayed in Figure 4.2.1. 


The functions x(t) and y(t) have the following properties: 


(1) a(t) and y(t) are, respectively, even and odd functions, i-e., 


a(—t) = x(t), y(—t) = —y(t). (4.19) 
(2) Both x(t) and y(t) are periodic, with period T = 4, ice., 
a(t+4)=<2(t), y(t+4)=y(t). (4.20) 


(3) a(t) and y(t) satisfy the bounds 
enh eii. = (5) < y(t) <+ (5) (4.21) 


(4) Special values 
ex(0)=1, 2(1)=0, 2(2)=-1, 2(3)=0, 2(4)=1 
ey(0)=0, y1)=—-3, y2)=0, y(3)= 5, (4) =0. 
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1.0 


-0.5 


- 1.0 


y(t) 
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- 0.2 
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FIGURE 4.2.1: Plots of x(t) and y(t) vs t. See Equations (4.9), (4.15), (4.17) 
and (4.18). 
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FIGURE 4.3.1: Definition of r and @ for the curve |y| + (5) 2? = 3. 


4.3 GEOMETRIC ANALYSIS OF |y| + (4) 2? = 4 


Figure 4.3.1 defines r and @ for this case. The fundamental algebraic rela- 
tions are 


is 4 
= == 4.22 
i+ (5)2=5 (4.22) 
rae ty’, (4.23) 
where 
x(0) =r(0)cosé, y(0) =r(0)sin@. (4.24) 


Since, by inspection of Figure 4.3.1, we have 


r(0)>0, O0<0<2r, 
r(—6) =r(9), (4.25) 
r(0+7) =r(9), 


and it follows that 7(@) and y(@) are, respectively, even and odd, i-e., 
a(—8) =2(9), y(—8) = —y(9), (4.26) 
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and are periodic, with period T = 27, i.e., 
a(6+2r)=2(6), y(O+2r) = (8). (4.27) 


The function r(@) can be determined by substituting the expressions in 
Equation (4.24) into Equation (4.22); doing this gives 


(cos 0)?r? + (2|sin 0|)r — 1 = 0. (4.28) 


This latter equation has, for a given 6, two solutions. Since r(@) > 0, we select 
the solution satisfying this condition and it is 


r(0 


wa 


1 
— |. 2 | gies —. : 
7 oe 2|sin | + /4|sin 6|? + 4(cos @) 


_ 1—|sind| 
~ (cos 0)?’ 
where the latter equation can also be written as 
1~—sin@ 
Ho) = | Gear 0<0s%, 


1+sin 0 
Teos 0)? TT < 7 < 27. 


(4.29) 


(4.30) 
Using (cos #)? = 1 — (sin @)?, Equation (4.30) can be rewritten as 


— <O0< 
r(9) = ¢ Tesin® OSOS", (4.31) 
Note that 
rT\ — 2 — 
‘i =. (3) =5, r(x) =1 (4.32) 


r() = 3, r(x) =1. 


The parabolic trigonometric functions are defined as follows: 


pbd(9) = r(9), 

pbc(?) = r(@) cos 0, (4.33) 
pbs(0) = r(0) sin @, 

where they are “called,” respectively, the parabolic dine, cosine, and sine func- 


tions. 
Graphs of these three functions are given in Figure 4.3.2. 


4.4  |y|— (})«?=} AS A DYNAMIC SYSTEM 


Consider the Hamiltonian 


A(a,y) =|y| — (5) 7 = 7 (4.34) 
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pbd(@é@) 


pbs(@) 


FIGURE 4.3.2: Graphs of the parabolic trigonometric functions. 
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which gives rise to following equations of motion 


de = SEN) = sgn(y), en 
dy _ _ 9H(xy) _ , 
dt Ox _ 
Two possible sets of initial conditions are 
1 
A:2x(0)=0, y(0)= 5 (4.36) 
1 
B:2(0)=0, y(0)=—- (5) . (4.37) 


As Figure 4.4.1 indicates, initial conditions A and B, respectively, correspond 
to motions along the upper and lower phase-space trajectories, i.e., 


1 1 
y>osy=5+(5) x, (4.38) 


y<O:y=- (5) - (5) x”. (4.39) 


The explicit solutions for initial conditions A, where y > 0, will now be 
calculated. For this case 


dx 
ae +1, «(0) = 0, 
cf 7 + (4.40) 
dt x, y(0) = 99 
and the solutions are given by the expressions 
1 1\ 5 


Likewise, for initial conditions B, where y < 0, we have 


3 =-1, x«(0)=0, (4.42) 


z(t) =-t, y(t)=—- (5) - (5) i, (4.43) 


If we extend ¢ to the interval, —co < t < oo, then z(t) and y(t) are, respec- 
tively, odd and even functions of t. These functions are plotted in Figure 4.4.2. 
Observe that both x(t) and y(t) are unbounded. 


with solutions 
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FIGURE 4.4.1: Graphs of H(z, y) = |y| — (4) x? = 4, along with the direc- 
tions of the phase-space flow. 
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1.0 a(t) 


0.5 


y(t) 


-1.0 -0.5 0.0 0.5 
FIGURE 4.4.2a: x(t) and y(t) vs t for 2(0) = 0 and y(0) 
—co <tE< oO. 
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top MO) 


0.5 


-O.5 


- 1.0 


FIGURE 4.4.2b: x(t) and y(t) vs t for 7(0) = 0 and y(0) = 
t<o. 


, and —oo < 


i 
2 


4.5 GEOMETRIC ANALYSIS OF || — (4) 2? =} 


The parabolic hyperbolic-type functions associated with 


wi- (3) =5 (4.44) 
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y 
1.0 


P(z.y) 


x 
- 1.0 
FIGURE 4.5.1: Definition of r and 6 for the curve |y| — ($) x? = $. 
will now be examined; see Figure 4.5.1. We have 
7? (0) = 0(0)" + y(0)?, 
x(0) = r(0) cos, (4.45) 
y(0) = r(6) sind 


Substituting (9) and y(0) into Equation (4.44) gives the following expression 
for r(0) 
(cos 0)?r? — (2|sin 6|)r +1 =0. (4.46) 


Inspection of Equation (4.46) indicates that it has either two positive roots 
or a pair of complex conjugate roots. Solving this quadratic equation for r(@) 
gives 


aoe Feod [sin 8| + Bin OP? — (cos 8] 


__ [sin | + \/(— cos 20) 


sos OF? (4.47) 
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where 


|sin 0|? — (cos)? = (sin)? — (cos 6)? 


. (5) (1 — cos 26) — (5) (1 + cos 26) 


= — cos 26, (4.48) 


was used to simplify the item within the square root. 
Observe that the results in Equation (4.47) allow the following conclusions 
to be reached: 


, . Sot ; 1 1 
(i) The ray from the origin can intersect the curve |y| — (5) 2? = 3 


where, once, or twice. 


(ii) In the following angular intervals, rays from the origin do not intersect 
the curve 


Tt 3m on) Tt 
<@<-, —<é0<—, — 27. 4.4 
VSO a 4 tsar Reet (4.49) 
(iii) For the angles 


4=-, —, —, — 4, 
A o] A ? 4 ? 4 ? ( 50) 
a ray from the origin is tangent to the curve. 


(iv) Given r(@), the smallest value is the one to be selected for the defini- 
tions of the parabolic hyperbolic functions, i.e., 


_ __ [sin @| — \/(— cos 26) 
r(0) =r_(0) = tees (4.51) 


With this in mind, the three parabolic hyperbolic functions are defined to be 


Pbhd(6) = r(0), (4.52a) 
Pbhc(0) = r(8) cos 9, (4.52b) 
Pbhs(0) = r(@) sin 9, (4.52c) 


where r(@) is only defined in the angular intervals 


T 37 51 as 
FOS GS fo = (4.53) 
It is important to point out that these three functions are periodic, with period 
27. But they are not defined for all values of 8; only for those given in Equation 
(4.53). The graphs of these functions are similar to those given in Chapter 1, 
Section 1.7. The three functions, defined in Equations (4.52), are respectively 
called the parabolic hyperbolic dine, cosine, and sine functions. 
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NOTES AND REFERENCES 


Section 4.1: The only relevant reference that I could locate that was 
directly concerned with parabolic-trigonometric functions was the following 
paper: 

1. G. Dattoli, M. Migliorati, M. Quattromini, and P. E. Ricci, The 

parabolic-trigonometric functions, arXiv: 1102.1563v1 [math.ph]. 


This work turned out to be of little value for my purposes. 
Sections 4.2—4.5: As far as Iam aware, my presentation has not appeared 
in the published works of others. 
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Chapter 5 


GENERALIZED PERIODIC 
SOLUTIONS OF f(t)? + 9(t)? =1 


5.1 INTRODUCTION 


The equation 
POH+gPt) =1 (5.1) 


has a range of solutions, both periodic and non-periodic. Further, solutions 
may or not be continuous, may or not have derivatives, etc. In this chapter we 
investigate a class of solutions to Equation (5.1) based on the assumption that 
both f(t) and g(t) have continuous second-order derivatives. The discussion 
presented is based largely on one of our previously published papers. 

As previously noted, in Chapters 1 and 2, two “elementary” solutions to 
Equation (5.1) are 


fi(t) = cos(Qt), gi (t) = sin(M2), (5.2) 
and 
fo(t) = cn(at), ga(t) = sn(at), (5.3) 


where 2 and a are real parameters, cosine and sine are the standard trigono- 
metric functions, and cn and sn are the Jacobi elliptic functions. Note that 
these functions are solutions to the following, respectively, linear and nonlin- 
ear, second-order differential equations 


a2 f 


f= (fio): oe +0? f =0, (5.4) 
and Pp 
f = (fasga): SE + Pf) =0, (5.5) 


where P3(f) is a cubic polynomial. 

Our major goals are to derive and analyze the fundamental properties of a 
particular class of periodic solutions to Equation (5.1) and to explicitly deter- 
mine the functional forms for these solutions. Again, as discussed previously, 
we find that associated with f(t) and g(t) is a third periodic function, d(t), 
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which is denoted the “dine function.” For the standard trigonometric situa- 
tion, the dine function is just the radius of the circle, and, as a consequence, 
has the constant value r(t) = 1, while for the Jacobi elliptic functions, it is the 
function dn(t). An important consequence of our work is how dine-type func- 
tions naturally appear in the investigation of periodic solutions to Equation 
(5.1), within the formulization of our methodology. 

The following is a summary of the central results found in this chapter: 


(i) For T real and positive, Equation (5.1) has periodic solutions of period 
T, ie., 


(5.6) 


f(t+T) = f(b), 
g(t +T) = g(t). 


(ii) The periodic functions, f(t) and g(t), are determined by a function, 
O(t), whose functional dependence on t is restricted by the requirement 
of periodicity. 


(iii) Both f(t) and g(t) are solutions to the same second-order, linear, non- 
autonomous ordinary differential equation. They differ, however, in their 
initial conditions, and this is reflected in the fact that f(t) can be con- 
structed such that it is an even function of t, while g(t) is odd in t. 


er 


The functions, f(t) and g(t), may be considered particular generaliza- 
tions of the standard trigonometric cosine and sine functions. 


(iv 


This chapter is organized as follows: Section 5.2 provides preliminary re- 
sults on the definition of functions, f(t) and g(t), in terms of a phase function 
@(t). Section 5.3 shows how 6(t) can be determined and gives its essential 
properties. An illustrative explicit example is studied in Section 5.4, including 
relevant graphs of the associated functions. In Section 5.5, we derive the linear, 
second-order, non-autonomous differential satisfied by the generalized cosine 
and sine functions, f(t) and g(t). The chapter ends with a brief discussion of 
non-periodic solutions of Equation (5.1). 


5.2 GENERALIZED COSINE AND SINE 
FUNCTIONS 


Let the real function 6(t) have a continuous second-derivative. Let c(t) and 
s(t) be defined as follows 


e) = c(t) +is(t), i= V1. (5.7) 


Now 
elt). e104) = c(t)? + s(t)? (5.8) 
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and we have 
c(t)? + s(t)? = 1. (5.9) 


Comparison of Equations (5.1) and (5.9) shows them to be the same equation, 
except for the labelng of the dependent functions. 
Since, for real 0(t), 


e°) — cos (t) + isin A(t), (5.10) 
we conclude that 
c(t) =cos@(t), s(t) = sin A(t). (5.11) 
For the case where 
20 
O(t) = (#) t, T>0O, (5.12) 


the standard trigonometric functions are obtained, i.e., 
Qat Qrt 
c(t) = cos (=) , s(t) =sin (=) (5.13) 


both functions of period T. Note that here we have c(t) and s(t), respectively, 
even and odd functions of t. To make this a specific feature, in general, for 
c(t) and s(t), we required 6(t) to be an odd function, i-e., 


6(—t) = —O(t). (5.14) 
This may be shown as follows 
c(—t) = cos 6(—t) = cos[—0(t)] = cos 6(t) = c(t), (5.15) 


and 
s(—t) = sin 0(—t) = sin[—@(t)] = — sin @(t) = —s(t). (5.16) 


If we now differentiate Equation (5.9), the following result is obtained 
c—+s— =0, (5.17) 


and this implies that 


5.18 
“= a(te(t), _ 


where, for the momentum, d(t) is not known. However, for the trigonometric 
cosine and sine functions, d(t) is a positive constant, not depending on f, i.e., 


d(t) =? > 0, (5.19) 


86 GENERALIZED TRIGONOMETRIC & HYPERBOLIC FUNCTIONS 
see Equation (5.4). This implies that d(t) should be selected such that 
d(t) >0, d(—t) = d(t). (5.20) 


In other words, d(t) is a non-negative, even function. We name d(t) the “dine” 
function associated with c(t) and s(t). 
Squaring the two expressions appearing in Equations (5.18) and adding 


gives 
de\*  (ds\” 
(=) + (=) = (ce +s”) =d’, (5.21) 


and this allows the determination of d(t) in terms of the derivative of c(t) and 
s(t); carrying out this calculation provides the expression 


a= f(#) + (2). 20 


Note that for the case where c(t) and s(t) are periodic, d(t) is also periodic. A 
nontrivial illustration of this is the Jacobi elliptic functions: cn(t), sn(t), and 
dni(t). 


5.3 MATHEMATICAL STRUCTURE OF 6(t) 


In addition to being odd, what other restrictions must be placed on 0(t) 
to ensure that c(t) and s(t) are periodic? Now, the periodicity conditions for 
c(t) and s(t) are the requirements 


c(t+T)=c(t), s(t+T) = s(t), (5.23) 
where T' > 0 is the period. This means that 


e(t4T) — oft + T) + is(t + T) 
= c(t) +is(t) 


Ses (5.24) 
Using 
1=¢7"*  & = integer, (525) 
it follows that 
eI (t+T) =f eA(t) = oot hi ’ ei9(t) = cilO(t) 427k] (5.26) 


Taking the “minimum period” corresponds to selecting k = 1; doing this gives 


6(t + T) = O(t) + 2n, (5.27) 
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which is a first-order, linear inhomogeneous difference equation. The general 
solution to this equation is 


6(t) = A(t) + (=) t (5.28) 
where A(t) has the properties 
A(-t) = A(t), A(t+T) = A(t). (5.29) 


Since 0(t) is assumed to have a continuous, second derivative, this condition 
applies also to A(t). One consequence of this analysis is that A(t) has a sine- 
type Fourier representation, i.e., 


~ Qnkt 
A(t) = )~ ay sin (=) ; (5.30) 
k=1 
where the coefficients satisfy the upper bound 
M 
lan] < => ye (5.31) 


where the positive constant, M, depends on A(t). Note that 0(t) is not peri- 
odic. In fact, we have 


O¢+T) = AE+T)+ (=) (t+T) 


2 
= A(t) + (F) t +20 
= 6(t) + 2r. (5.32) 
From Equations (5.11), we have 
del) = -(¢ a) sin O(t 
dt 
; (5.33) 


dst ant = = (#) cos O(¢ 


and using these expressions in Equation (5.22) gives the result 
. do 
[(cos 0)? + (sin 0)?] = ae" (5.34) 


However, without loss of (too much) generality, we will restrict the function, 
6(t), such that 
d6(t) 


ar 


> 0. (5.35) 
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Thus, it follows from Equation (5.28) that 


dA(t) 2m 
(jee ee 
= 4042 (5.36) 
Since d(t) > 0, the derivative of A(t) satisfies the condition 
dA(t) Qr 


If we had chosen to use the absolute value in the determination of d(t), 
Le., 
d0(t) 


d(t) = Se (5.38) 


then no restriction would be placed on dA(t)/dt. In the following, we will 
assume that the conditions in Equations (5.36) and (5.37) hold. 


5.4 AN EXAMPLE: A(t) = asin (#2) ¢ 


The simplest mathematical structure for A(t) is one for which only a single 
harmonic is included, i.e., A(t) has the form 


A acai (F) t, (5.39) 


and 0(t) is given by the expression 


A(t) + ay sin (=) poh (=) t, (5.40) 


] (5.41) 


J (5.42) 


d(t) = (=) i + a1 cos (=) J (5.43) 


The restriction d(t) > 0, requires that 


Jay] <1. (5.44) 
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It turns out that explicit expressions can be found for the Fourier series of 
c(t) and s(t). These results are based on the following two formulas 


cos(z sin t) )+2 3 J%(z) cos(2kt), (5.45a) 


sin(z cost) = 2 > Jox41(z) sin(2k + 1)t, (5.45b) 
k=0 


and the trigonometric relations 


cos|i (t) + $2(t)] = cos $1 (t) cos d2(t) — sin 61 (t) sin ¢o(t), (5.46a) 
sin[é1(t) + d2(t)] = sin d1(£) cos d2(t) + cos ¢1 (t) sin da(t). (5.46b) 


These Fourier series are (for T = 27) 
c(t) = cos(t + a, sint) 


= —Ji(a1) + S- [Jor—1(@1) — Jor+1(a1)] cos(2kt) 
k=l 


+ 3 [Jon (a1) + Jon+2 (a1)] cos(2k + 1)t, (5.47) 


s(t) = sin(t + aq sin t) 
Co 


=) [Jor—1(a1) + Jor41 (a1) sin(2kt) 
k=1 


+ dX Jox (a1) + Jor42(a1)] sin(2k + LE. (5.48) 
=0 


In the above relations, J,(z) is the k-th order Bessel function. Figures 5.4.1 
to 5.4.2 are illustrative plots of A(t), c(t), s(t), and d(t). 


5.5 DIFFERENTIAL EQUATION FOR f(t) AND g(t) 


The functions f(t) and g(t) are solutions to the same linear, second-order, 
non-autonomous differential equation. To demonstrate this, we start with 


a an 2 (5.49) 
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© (t) =t+ [0.2] sin (2) 
25¢ 


20; 


15; 


10; 


d(t) =1 + 0.2 cos(t) 
1.2 


1.15 


1.0; 


0.9} 


5 10 15 20 25 


FIGURE 5.4.1: Plots of 6(¢) vs t and d(t) vs t, for a, = 0.2. 
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c(t) = cos(0.2 sin(t) + t) 
1.0 


0.5; 


s(t) = sin(0.2 sin(é) + f) 
1.0} 


0.5} 


-0.5} 


-1.0} 


FIGURE 5.4.2: Plots of c(t) and s(t). 
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This is just a relabeling of Equations (5.18). Taking the derivative of the first 
expression gives 


p= —tg—ag =(-a) (4) -aay=(F)r-eF 6.50) 


Since d(t) = 0’(t), we finally obtain 
d' 
f"’- (5) fi+@f =0. (5.51) 
Starting with the second expression in Equation (5.49), we, likewise, obtain 
d' 
g’ - (5) g +d’g=0. (5.52) 


Note that Equations (5.51) and (5.52) are the same differential equation, ex- 
cept for the label placed on the function. Therefore, this situation can be 
summarized as follows: 


cp tt —. fdlh ! 2 ie 
f(t):w" — (4) w' +@w=0, is 
w(0)=1, w'(0) = 
watt — (da ! Di hee 
g(t): w" — (4) w! + Pw =0, a 
w(0)=0, w'(0) = 
Making use of the fact that d(t) = 6’(t), we have 
e a i 6')?w =0 5.55 
wi — (> }w' + (8')°w =0. (5.55) 


5.6 DISCUSSION 


We now summarize and discuss what has been obtained in this chapter: 
i) The functional equation, f?(t) + g?(t) = 1, has a class of T-periodic 
solutions that can be expressed as 


O(t) =cos@(t), s(t) =sin A(t), (5.56) 


where 6(t) has the form 


6(t) = A(t) + (=) t (5.57) 
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with 
oS 20k 
= i ——e t. . 
A(t) d ay, sin ( 7a ) (5.58) 


ii) A knowledge of 6(¢) allows for the complete determination of the three 
periodic functions: c(t), s(t), and d(t) = 6’(t). 

iii) c(t) and d(t) are even functions of t, while s(t) is an odd function. 
These results are a consequence of 0(t) being an odd function of t. 

iv) Given 6(t), then f(t) and g(t) satisfy the following coupled, first-order 
differential equations 


df(t) 


a —O'(t)g(t), ae O'(t) f(t), (5.59) 
and the following second-order equation 
2 [O"E]) dal (FO _ 
(ae [relat OF} (Ga) = i 


Examination of Equation (5.1) shows that it is invariant under t > —t, and 
such a condition also holds for the differential Equation (5.60). Since this equa- 
tion is linear, the previous result implies that its solutions can be constructed 
to be either even or odd functions of t. 

Observe that if 6(¢) = t, then Equation (5.60) becomes 


(S) iN)» 


f(t) =cos(t), g(t) = sin(¢). (5.62) 
v) The Jacobi elliptic functions, (cn, sn,dn), depend not only on t, but 
also a parameter k, where 0 < k < 1. However, in the discussion to follow, its 
functional dependencies, with regard to the Jacobi elliptic functions, are not 
important. 
Assume that 


with solutions 


6’(t) = dn(t), (5.63) 
then it follows that : 
O(t) = ? dn(z)dz, (5.64) 
0 
and 
ec) — C(t) + iS(t) = cos 6(t) + isin 0(t). (5.65) 
Therefore, 
t 
C(t) = cos if an(2)de| ; (5.66a) 
0 


S(t) =sin | i an(2)de| (5.66b) 
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and taking the respective derivatives of C(t) and S(t) gives 


aco = —dn(t)S(t), (5.67a) 
ao = dn(t)C(t), (5.67b) 


and these are the coupled, first-order differential equations for the Jacobi 
elliptic functions, i-e., 


C(t) =cn(t), S(t) = sn(t). (5.68) 


Thus, we conclude that the Jacobi cosine and sine functions are the real and 
imaginary parts of the complex function, exp[‘0(t)], where @(t) is related 
to the Jacobi dine function through the relationship expressed in Equation 
(5.64). This means that cn(t) and sn(t) are solutions to both a nonlinear, au- 
tonomous, second-order differential equation having a cubic nonlinearity and 
also solutions to a linear, non-autonomous, second-order differential equation. 


5.7 NON-PERIODIC SOLUTIONS OF /?(t) + 9°(t) =1 


Equation (5.1) has non-periodic solutions. An example for —oo < t < +00, 
is 


f(t) = 3. = ravonal number, (5.69a) 
0, = irrational number; 
iGy= 0; t= raven! number, (5.69b) 
1, ¢ = irrational number. 
Other elementary examples are (—co < t < +00) 
f®)=1, g(t) =0, (5.70) 
and 
als aay ill : 
f(t) =-1, g(t)=0, fort >0. 


Note that the functions, f(t) and g(t), are not continuous at t = 0. 
A more interesting set of functions are those that are oscillatory, but not 
periodic. So, taking 
A(t) = 0, (5.72) 


we have 
e() — cos(t?) + isin(t?), (5.73) 
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\ AAMIMM 
~VVTIANT 


sin(t+ ¢) 
1.0} 


“ZAMAN, 
= VVAITATN 


-1.0; 


FIGURE 5.7.1: Two oscillatory, non-periodic functions, corresponding to 
6(t) =# and O(t) =t+ #?. 
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ef 
: 2 4 6 8 10 12 


FIGURE 5.7.2: Two oscillatory, non-periodic functions, corresponding to 
O(t) =t?/(14+#?). 
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and 
f(t) =cos(t?), g(t) = sin(t”). (5.74) 

The effective period at time ¢ can be calculated. Start with ¢? and write it as 
follows 

Piao (5.75) 

T(t) ” , 
where T(t) is the effective period. Therefore, T(t) is 
2 
T(t) = -_ (5.76) 


and we find that T(t) decreases with T. The results shown in Figures 5.7.1 
and 5.7.2 illustrate this phenomena, for 0(t) = t? and 0(t) =t+t?. 


NOTES AND REFERENCES 


Section 5.1: This chapter is based on work presented in a previous article: 
R. E. Mickens, Periodic solutions of the functional equation f(t)? + g(t)? = 1, 
Journal of Difference Equations and Applications 22 (2016), 67-74. 

Section 5.2: Starting from c? + s? = 1, we have the result 


on taking the derivative. If we wish to split this into two expressions, one each 
for dc/dt and ds/dt, then a possibility is 


dc ds 
de = —d(t)s, dt = d(t)s. 


The simplest choice is that d(t) is a constant. But, as we show in this section, 
the more interesting and more productive selection is having d be a function 
of t. It should be noted that in Hamiltonian dynamics d(t) is taken to be one. 

Section 5.3: Methods to solve the first-order difference Equation (5.27) 
are given in R. E. Mickens, Difference Equations: Theory, Applications and 
Advanced Topics, 3rd edition (CRC Press; Boca Raton, FL; 2015), Section 2.2. 

The results expressed in Equations (5.30) and (5.31) may be justified using 
techniques from the theory of Fourier analysis. See, for example: R. E. Mick- 
ens, Mathematical Methods for the Natural and Engineering Sciences (World 
Scientific, London, 2004), Section 2.5. 

Section 5.4: The mathematical relations, given in Equations (5.45) are 
from I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series and Prod- 
ucts (Academic Press, New York, 1965). See formulas 8.511.3 and 8.511.4. 
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Taylor & Francis Group 
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Chapter 6 


RESUME OF (SOME) PREVIOUS 
RESULTS ON GENERALIZED 
TRIGONOMETRIC FUNCTIONS 


6.1 INTRODUCTION 


This chapter provides a very brief summary of work on generalized trigono- 
metric functions and related topics. We do not seek completeness in terms of 
all the various formulations and results which have appeared in the scien- 
tific and mathematical research literature. Upon examining these works, it 
becomes clear that there is not full unity with regard to the definition of 
the generalized sine and cosine functions. This is also demonstrated in the 
present work of the author within this book. Consequently, the main task of 
the chapter is to introduce and summarize some of the publications that have 
influenced my views of generalized trigonometric functions. 

Section 6.2 approaches the generalized sine and cosine functions by defin- 
ing them as solutions to two coupled, nonlinear, differential equations, with 
specified initial conditions. Many physical systems are formulated this way, in 
particular, systems which can be modeled by Hamiltonians. 

Section 6.3 introduces the generalized sine and cosine by way of an integral 
representation. This way of proceeding is similar to how the standard sine and 
cosine functions are “constructed” in a calculus course. 

Section 6.4 is devoted to defining generalized sine and cosine functions in 
terms of the geometric properties of curves in a two-dimensional plane. We 
show that in addition to the usually expected sine- and cosine-type functions, 
there is an associated third periodic function, the “dine.” This origin of the 
dine function was explored in all of the previous five chapters. 

Section 6.5 is not directly connected to most of the published articles on 
the generalized sine and cosine functions, discussed in Sections 6.2, 6.3, and 
6.4. However, the use of symmetry considerations has played an important role 
in the work of the author, and in this section, we show how the invariances of 
the phase-space curves can be used to deduce very strict constraints on the 
properties of the generalized sine and cosine functions. 

Finally, Section 6.6 provides a general, but brief discussion of the sinp, and 
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COSpq functions, and their relationship to the standard trigonometric sine and 
cosine functions. 


6.2 DIFFERENTIAL EQUATION FORMULATION 


The standard trigonometric sine and cosine functions may be defined by 
means of the following system of coupled, first-order, linear differential equa- 


tions 
dx(t) 


¢ = —y(t), 
dy(t) = x(t), (6.1) 
x(0)=1, y(0)=9, 
where 
x(t) =cost, y(t) =sint. (6.2) 


Note also that both x(t) and y(t) satisfy the following set of second-order, 
linear differential equations 


d*x(t) 7 _ dx(0) _ 
ae a(t)=0, 2x(0)=1, rm 0, (6.3a) 
Tutt) +y(t)=0, y(0) =0, au) =, (6.3b) 


Also, if the first and second expressions in Equation (6.1) are multiplied, 
respectively, by x(t) and y(t), and integrated, then we obtain 


z(t)? + y(t)? = 1, (6.4) 


which corresponds to the Pythagorean identity. 
These results may be generalized (Shelupsky, 1959) by considering the 
following system of coupled, nonlinear differential equations 


aot) — _ f(y), 
du) — f(x), (6.5) 


where 
fz)=|eP 2, p> (6.6) 
Since 4 r 
f= =-f@FW), FW =fW)FO), (6.7) 


it follows that 
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If we now make the definitions 
x(t) =sinp(t), y(t) = cosp(t), 
then the following generalized Pythagorean identity is obtained 
|sin,(t)|? + |cos,(t)|? = 1. (6.9) 


Thus, this definition of the generalized sine and cosine function is based on 
the differential Equations (6.5), which in turn gives rise to the first-integral 
presented in Equation (6.9). 

The above arguments can be extended to cover the situation for which the 
corresponding sine function depends on two parameters, p and q. This gives 
rise to the identity 

Ising (t)!? + [cospy(t)| = 1, (6.10) 


which can be used to define cos,,(t), ie., 


COSpq (t) = (+) [1 — |sinpg(t)|"]/%. (6.11) 


6.3. DEFINITION AS INTEGRAL FORMS 


Generalized trigonometric functions may also be defined in terms of inte- 
grals. This method is based on an extension of the definition of the regular 
sine function, i.e., 


edz 
t= | —~ 0<|e/<1, 6.12 
fs ox (6.12) 
and 
x(t) =sint, (6.13) 
to the expression 
t dz 
; 0O<t<1, 
sing} = fo Oy: (6.14) 
(“) Jo Gasper ~ bse sO, 


where cos,(t) is related to sin,(t) by the relationship 
|sin,(t)|? + |cos,(t)|? = 1. (6.15) 


Two parameter generalizations have been constructed for the sine function. 
An example of such an extension is 


t dz 
OS i Gaenr0 StS, 


= eee. (6.16) 
Ae qa st, 
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with 
Isinpa(t)? + [co8pq(¢)[? = 1. (6.17) 


Generally, the (p,q) are taken as positive. However, particular requirements 
may force p> 1 andq> 1. 


6.4 GEOMETRIC APPROACH 


In Chapters 1 to 4 of this book, we have investigated the construction of 
periodic solutions to the functional equation 


lyl? + ||? =1, (6.18) 
where x and y are functions of an angle 9, ice., 
cx=2(0), y=y(0), O<A< 2r. (6.19) 


This general technique was also used by Mahdi et al. (2014) to study the 
properties of functions x(6) and y(@) satisfying the constraint 


jal" + |yl™=1, m>O. (6.20) 


As we have shown, the derived functions are always periodic with period 27, 
but some of the periodic solutions may be unbounded. 
Considering the relation in Equation (6.20), « and y are replaced by 


x(0) =r(0)cos0, y(0) = r(@)siné. (6.21) 


Therefore, substitution of these representations into Equation (6.20) and solv- 
ing for r(@) gives 


1 
— 22 
ava [|cos |” + |sin A)” ]1/™? 24) 
and we have 9 
cos 
= 008 (0) = —— 2 
7(8) = cosin(®) = Tap + jain dm oer 
sin 0 
8) = sin», (0) = —————__, 6.24 
ye) she) [|cos 6|™ + |sin O|™]1/™ ea 
1 
0) = din, (¢) = ————_—__-. 2 
nie) = atal?) [|cos 0] + |sin o|7]1/™ eee) 


Note first that there are three basic periodic functions. Second, cos,,(@) and 
sin;,(@) have a period 27, the related dine function, din,,(@), has a period of 
Tr. 
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6.5 SYMMETRY CONSIDERATIONS AND CONSE- 
QUENCES 


This section gives a brief summary of results which can be derived from 
the application of symmetry transformations to the formula 
yP+|eit=1, @21,¢2 1), 
x(0)=1, y(0) =0, 


where x = x(t) and y = y(t). This means that Equation (6.26) is taken to 
represent a dynamic system in the x-y phase-plane, and ¢ represents the time. 


(6.26) 


6.5.1 Symmetry Transformation and Consequences 


Consider the following transformation of the x and y variables: 


T1:287-@, yoy, (6.27a) 
To:t72, yr-y, (6.27b) 
T3 = T,T> = ToT; 22> Lh, Y>mey. (6.27c) 


T, corresponds to reflection in the y-axis; T> to reflection in the x-axis; and 
T3 to inversion through the origin. For 


A(x, y) = yl? + |al*, (6.28) 


we have 
T1H(z,y) = H(-2,y) = H(z,y), 
12H (z,y) = H(z,—-y) = H(z,y), (6.29) 
T3H(z,y) = H(—2, —-y) = H(2,y). 


Thus, it can be concluded that H(«,y) is invariant under the three transfor- 
mations (71, T2, 73). This fact allows the following conclusions to be made: 


(a) Equation (6.26) is a simple, closed, convex curve in the z-y phase-plane. 


(b) Considered as a dynamic system, x(t) and y(t) are periodic, i.e., there 
exists a T;,, such that 


u(t+Tpq) =2(t), y(t + Tra) = yb). (6.30) 
(Ipq will be calculated in Section 6.5.3.) 
(c) a(t) and y(t) satisfy the bounds 
-—l[S<2)s), ~legG)s L (6.31) 
(d) x(t) and y(t) are, respectively, even and odd functions of ¢, i-e., 


a(—t)=a(t), y(—t) = —y(). (6.32) 
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6.5.2 Hamiltonian Formulation 


Let H(a,y) be considered the Hamiltonian of a dynamic system; see Equa- 
tion (6.28). The equations of motion are 


dx _OH(x,y) dy __ OH(x,y) 


dt »=6 Oy. sot CS Ox? 8) 
and from these relations the following expressions result 
dx ot 
HT PlyP sony), 2(0) = 1, (6.34) 
dy 1 
where sgn(z) is the “sign” function, i-e., 
+1 2z>0, 
sgn(z)=40 z=0, (6.36) 
-l1 z<0. 


The first-order differential equation for the trajectories in the x-y phase- 
space, y = y(x), is given by the equation 


dy | En 
Ls oe 6.37 
dx -) piylP-* | Lsgn(y) een 
y 
oo a 
(O,-1) 


FIGURE 6.5.1: Plot of the phase-space trajectory in the fourth-quadrant 
beginning at (1,0) and ending at (0,—1). The time to do this is T,,/4. 
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6.5.3 Area of Enclosed Curve 


Let us now calculate the area interior to the closed curve represented in 
Equation (6.26). Note that because of the invariance under the symmetry 
transformations, (T,, 72,73), this area is four times the area included in the 
first quadrant of the x-y plane, i.e., 


1 
area = A = 1/ y(x)dx, (6.38) 
0 
where 
y(z) = +[L— [el]. (6.39) 
Therefore, 
1 
A= if (1 — |a|9]!/? de. (6.40) 
0 
Now, let 
1/ 1 i-4 
u=e'>ar=ul? => dz —Ju du. (6.41) 
q 
With this change of variable, Equation (6.40) becomes 
4 1 ag 1 
A=(- , ut (1—u)/Pdu. (6.42) 
G/ Jo 
Using the definition of the beta function 
: DP(a)T'(b) 
B(a,b) = ety) t= 6.43 
(a0) = fw a— wrt = ee (6.43) 


where I'(z) is the gamma function, we obtain for the area the result 


r(+)r(1+4 
A= (2) Pa)ross) (6.44) 
q) 7 (1 +24 1) 
Using the relationship 
I(z+1) = 2T(z), (6.45) 


we finally have 


(6.46) 


and the area enclosed by the curve |y|? + |z|2 = 1. 
Observe that the area, A = A(p,q), is symmetric with respect to p and q, 
Le., 
A(p, 4) = A(q, P). (6.47) 
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6.5.4 Period 


The time taken to travel from (1,0) to (0, —1) in phase-space, is one-fourth 
of the period T,,. This result is a consequence of the symmetry properties of 
the phase-space trajectories; see the discussion presented in Section 6.5.1. 


Along this portion of the curve 
y(x) =(-1)[1—lalt"?, 20. 


Therefore, 
dx p-1 i 
—, =ply|Psgn(y) = (—)p [1 — |2|*] 


dt 
and it follows that 
( 1 ) dx 
—dt = ~~ yp? 
p [1 — x4| pe 


and 


Tpq/4 0 
-{ ee (=) / _ ats, 
0 Pp 1 [l— 29]? 


Therefore 


4 1=p 
Log = | — | 1—a%|? dz, 
(5) ) 


which with the transformation of variable 


becomes 


1 — #% ie 
tn= (4) fa-wtuian= (4) r()PG) 
pa) Jo pq r (+2) 


p—-1 
Pp 
? 


(6.48) 


(6.49) 


(6.50) 


(6.51) 


(6.52) 


(6.53) 


(6.54) 


Note that, just as for the area formula, T,, = T (p,q) is symmetric in p and q, 


i.e., 
T (p,q) =T(4,P)- 


(6.55) 


From the above analysis, it follows that x(t) and y(t), which we can denote 


as 
r(t) = coSpq(t), y(t) = sinpg(t), 


have the following special values: 


x(0) =1 (0)=0 

x Toa =0 y Tea =-1 

x Fea =-l y Toa =0 
3T pq — 3T pq — 

u(t) = y(—)=1 


(6.56) 


(6.57) 
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6.6 SUMMARY 


The expressions in Equations (6.10), (6.18), and (6.26) may be used to 
define generalized trigonometric sine and cosine functions. They differ in what 
is considered the independent variable for the x and y functions. Of importance 
is the fact that the general relation 


fe +|2|?=1, 2(0)=1, y(0)=0, (6.58) 


p21, ql, 


is invariant under the transformation stated in Equation (6.29). Since the 
standard sine and cosine trigonometric functions correspond to the case where 
p=q= 2, it is not at all surprising that the generalized sine and cosine func- 
tions have similar properties, as reflected in the results stated in Sections 6.5.1 
and 6.5.4. 


NOTES AND REFERENCES 


Section 6.1: An excellent introduction to “sine” functions is the book 
listed below. It covers the generalized sine, as taken by the author to be the 
inverse function defined by an integral of the form 


[ dz 
) — SS ——————s 
0 Vl+tmz24+nz4 
The four cases considered are 
e circular sine: m = —1, n = 0; 
e hyperbolic sine: m = 1, n = 0; 
e lemniscate sine: m = 0, n = —1; 
e sine amplitude of Jacobi: m = —(1+k?), n=k?, forO<k <1. 


(1) A. I. Markushevich, The Remarkable Sine Function (Elsevier, New York, 
1966). 


Section 6.2: The topics of this section are discussed in the references: 


(1) D. Shelupsky, A generalization of the trigonometric functions. The 
American Mathematical Monthly 66 (1959), 879-884. 
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(2) D. Wei, Y. Liu, and M. B. Elgindi, Some generalized trigonometric 
sine functions and their applications. Applied Mathematical Sciences 6 
(2012), 6053-6068. 


(3) R. M. Rosenberg, The Ateb(h)-functions and their properties. Quarterly 
of Applied Mathematics 21 (1963), 37-47. 


Section 6.3: In addition to ref. [2], given under the items in Section 6.2 
above, see also the following: 


(1) J. Lang and D. Edmunds, Figenvalues, Embeddings and Generalised 
Trigonometric Functions (Springer, Heidelberg, 2011), Section 2.1. 


Section 6.4: A good summary of the “purely” geometric approach is given 
in the paper: 


(1) H. Mahdi, M. Elatrash, and S. Elmadhoun, On generalized trigonometric 
functions. Journal of Mathematical Sciences and Applications 2 (2014), 
33-38. 


Section 6.5: The geometry of two-dimensional (x — y) phase-space is 
discussed in the following books: 


(1) M. Humi and W. Miller, Second Course in Ordinary Differential Equa- 
tions for Scientists and Engineers (Springer-Verlag, New York, 1988). 
See Chapter 8. 


(2) R. E. Mickens, Mathematical Methods for the Natural and Engineering 
Sciences, 2nd Edition (World Scientific, London, 2017). See Sections 4.4 
and 4.5. 


Chapter 7 


GENERALIZED 
TRIGONOMETRIC FUNCTIONS: 


yl? + [elt =1 


7.1 INTRODUCTION 


This chapter presents and discusses a method for determining periodic 
solutions to the functional equation 


ly|P+|2l?=1, p>0, g>0. (7.1) 
We show that functions x(t) and y(t) can be selected such that 
e(t+T)=<a(t), yt+T)=y(t) (7.2) 


for any T > 0. Further, our procedure constructs x(t) and y(t) that mimic 
the standard trigonometric cosine and sine functions, i.e., x(t) and y(t) are, 
respectively, even and odd, and bounded. Further, x(t) and y(t) have zeros 
and maximum/minimum values at the t values (0,7/4,T/2,3T/4,T). 

Section 7.2 outlines the general methodology, while the explicit functions 
are given in Section 7.3. Finally, in Section 7.4, we present a gallery of plots 
of a(t) and y(t) for a selected set of (p,q) values. 


7.2 METHODOLOGY 


Our purpose is to construct explicit solutions to Equation (7.1) such that 
x(t) and y(t) have exactly the same essential properties as the standard cosine 
and sine functions. This means that the following hold: 


(i) 
x(0)=1, y(0)=0. (7.3) 


(ii) 
x(t) = a(t) y(—t) = -y(2). (7.4) 
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(iii) For given T > 0, the periodicity condition of Equation (7.2) holds. 


(iv) At the indicated specific values of t, we have 


2(0)=1, «(%)=0, 2(#)=-1, 

Ber (7.5) 

nA =O, ah) 1, 
y(0) a (z) > ¥(s) =9, (7.6) 

y()=-1, y7)=0 

Now, make the following transformation of dependent variables 

ue)? = |a@)I*, 7a) 
[v()]? = lwo? b) 


and note that Equation (7.1) can be written as 
[v(t)}? + [u(t))? = 1, (7.8) 


in the new variables. The work presented in Chapter 5 informs us that there 
are an unlimited number of periodic solutions to Equation (7.8) provided that 
u(t) and v(t) take the form 


u(t) =cosy(t), v(t) = sin v(t), (7.9) 
where 
p(-t) = —v(t), 
ee W(t) +2n cer 


and the general solution for 7)(t) is 


=yom sin (+ =~) + (F) f. (7.11) 


As a consequence of these results, z(t) and y(t) can be selected as follows 
a(t) = [cos¥(t)]?/4 sgn(cos w(t), (7.12a) 
y(t) = [sine (t)]?/? — sgn(sin(t)). (7.12b) 
Comments: 


(a) Observe that in Equation (7.1), the absolute value operation appears. 
This takes care of potential problems of having the sign functions ap- 
pearing in Equation (7.12). 
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(b) When actually computing either x(t) and y(t), from the expressions 
given in Equations (7.12), the following mathematical formats must be 
used 


[cos h(t) 2/4 + {[cosa(t)2}/" , (7.13a) 
[sin o(t))?/” > {[sin v(t))2} 7". (7.13b) 


Doing this prevents the appearance of complex values at an intermediate 
step of the calculations. 


7.3 SUMMARY 


For purposes of convenience, we give below a concise summary of the 
results derived in Section 7.2. Also, hopefully, without generating confusion, 
we rewrite and renumber several previously given equations. 

The functional equation 


Hi + [el = 1, cn 
p>0od, q>Qd0, 
has the following class of periodic solutions 
Lpq(t, T) = [cos w(t)]?/4 sgn(cos w(t), (7.15) 
Ypq(t,T) = [sinv(t)}?/” — sgn(sin Y(t), (7.16) 


where w(t) has a Fourier sine series as given by Equation (7.11). 

Note that in Equations (7.15) and (7.16), the dependence on (p,q,T) is 
indicated explicitly. 

It is important to note that xpq(t,T) and ypq(t,7) were explicitly con- 
structed to mimic the trigonometric cosine and sine functions. However, unlike 
most of the previous work on generalized trigonometric cosine and sine func- 
tions, our method of construction is not based on an a priori inverse integral 
relation. 


7.4 GALLERY OF PARTICULAR SOLUTIONS 


We end this chapter with a gallery of solutions for our class of periodic 
functions. On the following pages, we give plots of |y|? + |z|4 = 1, and zp,(t) 
and Ypq(t) versus t, for selected (p,q) values. 
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FIGURE 7.4.1a: Plot of |y|':? + ||? =1, ie., p= 1.2 and q=2. 
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FIGURE 7.4.1b: Plots of 21.2 2(t) and y1.2,2(t), versus t. 
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FIGURE 7.4.2a: Plot of |y|? + |z|4/3 = 1, ie., p = 2 and q = 4/3. 
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FIGURE 7.4.2b: Plot of #2,4/3(t) and y2,4/3(t), versus t. 
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FIGURE 7.4.3a: Plot of |y|? + |z|* = 1, i.e., p= 3 and q=4. 
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FIGURE 7.4.3b: Plot of x3,4(¢) and y3,4(t), versus t. 
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Chapter 8 


GENERALIZED 
TRIGONOMETRIC HYPERBOLIC 
FUNCTIONS: |y|? — |x|? =1 


8.1 INTRODUCTION 
We will use the functional equation 
ly? — |x|? =1 (8.1) 


to construct generalized hyperbolic functions, in correspondence with the stan- 
dard trigonometric hyperbolic functions, cosh(at) and sinh(at), 


ewt + e7at 


hyperbolic cosine(at) = cosh(at) = 5 ; (8.2) 
hyperbolic sine(at) = sinh(at) = a a (8.3) 

These functions satisfy the functional relation 
[cosh(at)]? — [sinh(at)]? = 1, (8.4) 


which is a special case of Equation (8.1) with p = 2 and q = 2. 
Note that the cosh and sinh functions have the following properties: 


(i) cosh(at) and sinh(at) are, respectively, even and odd, i-e., 


cosh(—at) = cosh(at), sinh(—at) = — sinh(at). (8.5) 


(ii) Since cosh and sinh are continuous functions, we have (as a trivial con- 
sequence) 


cosh(0) =1,  sinh(0) = 1. (8.6) 


(iii) The following asymptotic relations hold (for a > 0) 


t3 +o: eee * (2) ae (8.7) 
F (8.8) 
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8.2 SOLUTIONS 


Denote by ypq(t) and z,,(t) the generalized hyperbolic cosine and sine 
solutions of Equation (8.1). We wish these functions to have the following 
properties: 


(a) 
Ypq(0) =1, &pq(0) = 0. (8.9) 
(b) Ypq(t) = 1 and ap,(t) = 0 are to be, respectively, even and odd functions 
of t, i-e., 
Ypq(—t) = Ypq(t), Xpq(—t) = —Zpq(t). (8.10) 


Carrying out a similar set of substitutions and related procedures as was done 
in Chapter 7, the following expressions are gotten for ypq(t) and &pq(t): 


Ypq(t) = [cosh(at)]?*/”, (8.11) 


2pq(t) = {[sinh(at)]?}"/" [sgn(t)], (8.12) 


where a is taken to be positive. Note that the expression, enclosed in {... }, 
for Zpq(t) is written in a manner such that x,,(t) is always real. 

Inspection of Equations (8.11) and (8.12) show that they satisfy the a pri- 
ori desired properties expressed in Equations (8.9) and (8.10). These functions 
also have the following asymptotic behaviors 


Ypq(t) ~ ( 


2 8.13 
Lpq(t) ~ ( 


too: | 


t > (—)oo: ee (3)""e a es (8.14) 


8.3 GALLERY OF SPECIAL SOLUTIONS 


We now present a gallery of yp, and %pq for the same values of (p, q) as pre- 
sented in Section 7.4 of the previous chapter. Observe that from a “pictorial” 
viewpoint, the graphs are similar in shape. All plots are for a = 1. 
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FIGURE 8.3.1a: Plot of |y|+? — |a|? =1, for p = 1.2 and q = 2. 
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FIGURE 8.3.1b: Plots of yi.2.2(t) and 2.2,2(t), versus t. 
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FIGURE 8.3.2a: Plot of |y|? — |a|4/ = 1, for p = 2 and gq = 4/3. 
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FIGURE 8.3.2b: Plots of y24/3(t) and 22,4/3(t), versus t. 
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FIGURE 8.3.3a: Plot of |y|? — |z|* = 1, for p = 3 and q = 4. 
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FIGURE 8.3.3b: Plots of y3,4(t) and x3,4(t), versus t. 


Chapter 9 


APPLICATIONS AND 
ADVANCED TOPICS 


9.1 INTRODUCTION 


This chapter examines a number of issues related to conservative dynam- 
ical systems for which the Hamiltonian takes the form 


A(x,y) = ly? + |a|? = 1 
PAM tg 24 (9.1) 
x(0)=1, y(0) =9, 
where + = x(t) and y(y(t). The independent variable t is the time. The cor- 
responding equations of motion are 


dx _OH(x,y) dy OH(x,y) 


Se 9.2 
dt Oy ° dt Ox we) 
and 
dx pA 
& =lyl*sgn(y), (0) = 1, (9.3a) 
d 
<i = —gle|**sgn(x), y(0) =0. (9.3b) 
Observe that for p = q = 2, Equations (9.1) and (9.3) take the form 
A(z,y)=y+2° =1, 
ie — dy __ (9.4) 
ad Yy, dt a, 
while for p = q = 1, we have 
H(x,y) =|yl + lal =1, Ae 
& = sgn(y), $f = —sgn(z). 


Note that Equations (9.3) are invariant under the dependent variables 
transformation 


(9.6) 
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We denote such systems as having odd parity, and call them odd-parity sys- 
tems. For the general case of a second-order differential equation, which we 


write as 
F(a, dx/dt, d’x/dt?) = 0, (9.7) 


it is of odd parity if and only if 
F(—a, —dx/dt, —d?x/dt?) = (+)F (a, dx /dt, d?x/dt?). (9.8) 


The following are examples of odd-parity systems defined in terms of differ- 
ential equations 


: : ia 
harmonic oscillator : —- + « = 0, (9.9) 
dt? 
do 3 
Duffing oscillator : qatar + bya" = 0, (9.10) 
d2 
van der Pol oscillator : a +2 =e(1- ye, (9.11) 


where (a,b,¢€) are parameters. In Section 9.3, we show that for odd-parity 
periodic systems, the Fourier expansions have only odd harmonics. 

For many dynamical systems, modeling a broad range of phenomena in 
the natural and engineering sciences, Equation (9.7) takes the special form 


x dx 


and the system is of odd parity if and only if 


f (- -<) =-—f («. =) : (9.13) 


See the expressions in Equations (9.9), (9.10), and (9.11). 

We now present arguments to indicate why the harmonic oscillator, Equa- 
tion (9.9), and the Duffing oscillator, Equation (9.10), appear so often as 
mathematical models of dynamical systems. To begin, assume that x is the 
extension of a one-dimensional system away from a position of equilibrium. 
Further, assume that based on Newton’s force law, we have 


dex 
=I). (9.14) 


Expanding f(a) in a Taylor series gives 


d2 
7 = fot fix + fox? + fax? + O(c*) (9.15) 


where 


k; 
Le (=) ei k = (0,1,2,...). (9.16) 
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Now fo is zero, since our reference point is an equilibrium point. Also, for 
fo = 0, the system is of odd parity. Therefore, neglecting terms of O(x*) and 
higher, we obtain 
d’z 
dt? 
which is the Duffing equation. Note that the harmonic oscillator equation is 
the special case for when fs = 0. 
If y is taken to be y = da/dt, then in the z-y phase-plane 


+ fiet fax® =0, (9.17) 


dz dy(x) 
de = y(z) de? (9.18) 
and Equation (9.17) can be written as 
d 
Vs + fit t+ fax? =0. (9.19) 


Integrating this expression gives 
= il 
A(z,y) = (5) y+ (4) a + (2) a* = constant, (9.20) 


which is the energy function or Hamiltonian of the system. Therefore, the 
equations of motion are 


da OH (x,y) 


dt Oy =% (9 21) 
B= —UpeP = — fie — fae’, 


which is equivalent to the second-order differential Equation (9.17). 
Observe that if f; = 0, then H(x,y) becomes 


H(2,y) = (5) y+ (2) az = constant (9.22) 


and with a rescaling of the x and y variables, this becomes Equation (9.1) 
with p= 2 andq=4. 


For the remainder of this chapter, we will investigate dynamical systems 
for which the Hamiltonian has only two terms, both of which are power laws 
of their respective variables. Thus, we directly demonstrate the importance of 
generalized, periodic, trigonometric functions for the modeling, analysis, and 
understanding of this class of dynamical systems. 

In Section 9.2, we give arguments to show that odd-parity periodic sys- 
tems have Fourier expansions for which only the odd harmonics or frequen- 
cies appear. Section 9.3 introduces the concept of a “truly nonlinear oscillator 
(TNLO)” and we examine and solve special cases of this type of oscillator. In 
Section 9.4, the so-called Ateb periodic functions are introduced and briefly 
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discussed. Discretization of differential equations plays an important role in 
the numerical integration of such equations. Section 9.5 presents exact dis- 
cretizations of the differential Equation (9.17) for which the corresponding 
solutions are the three Jacobi elliptic functions. The next two Sections, 9.6 
and 9.7, are concerned with the method of harmonic balance and its gener- 
alization. This technique is of great value when accurate analytical approxi- 
mations to periodic solutions are required. Another methodology created to 
calculate approximations to the periodic solutions of the equations of motion, 
see Equations (9.3), is based on iteration procedures. In Section 9.8, we outline 
these types of procedures and illustrate their use by applying these methods 
to several examples. 

Finally, the main reason for the necessity of this chapter is the fact that in 
general, exact analytical solutions to Equation (9.1), in terms of finite combi- 
nations of “elementary functions,” are not possible. Therefore, we must either 
obtain numerical representations of the solutions and/or construct analytical 
approximations. The latter issue is the focus of this chapter. 


9.2 ODD-PARITY SYSTEMS AND THEIR FOURIER 
REPRESENTATIONS 


In this section, arguments are presented to show that for odd-parity sys- 
tems the Fourier series for periodic solutions only contain odd (angular) fre- 
quencies. 

To begin, consider the following dimensionless Duffing equation 


Tz +e +r? =0 
di : (9.23) 


da(0 
a(0)=1, S2=0. 


For 0 < ¢ < 1, the corresponding perturbation solution is 
€ 


x(6,€) = cosO+ (= 


) (— cos @ + cos 36) 


2 
+ (=) (23 cos 8 — 24 cos 30 + cos 50) 


+ O(e?), (9.24a) 


and 


O(e,t) =w(e)t = 1+ = -—~+O0(e)|, (9.24b) 


3e = 21? 
8 256 


where the angular frequency, w, is related to the period, T, by the relation 
2a 


7 (9.25) 


i 
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Observe that, at least to terms of O(e?), only odd angular frequencies occur. 
Likewise, for the van der Pol oscillator 


ax dx 


which is of odd parity, the perturbation solution for its unique limit-cycle 
solution is 


x(0,€) =2cos0 + (=) (3 sin 6 — sin 30) 


2 
+ (=) (—13 cos 6 + 18 cos 30 — 5cos 56) 
+ O(e*), (9.27a) 
where 
Pe 
O(e,t) =w(e)t = ae + o)| t. (9.27b) 


Again, only odd angular frequencies appear in the solution. 
We now demonstrate that this is a general feature of odd-parity periodic 
systems, i.e., if x(t) is such a solution, then x(t) has a Fourier representation 


x(t) = are cos(2k — 1)wt + By, sin(2k — 1)wt]. (9.28) 


The demonstration is based on the following three requirements: 


(i) The equation of motion 


ax dx 


is of odd parity and has a periodic solution. 


(ii) In the two-dimensional phase-space, (x,y = dx/dt), the periodic solu- 
tions occur about the fixed point, (0,0). 


(iii) The periodic solutions of Equation (9.29) are essentially unique, i-e., if 
x = $(t) is a non-trivial periodic solution, then for to 4 0, z = o(t — to), 
is also a periodic solution. 


From the perspective of the phase-space, the moving point 


(a(d.uto) = (0, 2) (9.30) 


132 GENERALIZED TRIGONOMETRIC & HYPERBOLIC FUNCTIONS 


traces out a closed curve. This means that 
dz(t)\ _ do(t — to) 
(20, 2) = (oe to), SE), (9.31) 


traces out the same closed curve, except that it is shifted in phase. 
Now assume Equation (9.29) has the complex Fourier representation 


x(t) = Slane" +a,e"], (9.32) 


where the a, are complex valued coefficients. Note that if x(t) is a periodic 
solution, then z(t) is also a periodic solution, where z(t) is defined to be 


shee (« + ) (9.33) 


This result is a consequence of the fact that both a(t) and —2(t) are solu- 
tions, and thus x(t — to) and —2(t — to) are also periodic solutions. From the 
uniqueness requirement 


z(t) = a(t) (9.34) 
and 


: (« + 5) ee) (9.35) 


If Equation (9.32) is substituted into this last expression and a comparison 
is made of the coefficients of the two exponential terms, then the following 
result is found 

(—1)¥a, = —ax, (9.36) 


and this implies that all the even labeled coefficients are zero, 
@om = 0, m= (1,2,3...). (9.37) 


Therefore, x(t) has a Fourier representation only in terms of odd harmonics 
or angular frequencies; see Equation (9.28). 

This result will play an important role in the formulation of appropriate 
methods to determine approximate analytical solutions for odd-parity sys- 
tems. 

A final comment. Consider the function N (a), defined as 


a 2 yan ; 


HI 


(9.38) 
=O. V2) eg Te = "0, 12) cee 


A direct calculation shows that N(a) is an odd function, i-e., 


N(x) = —N(x), (9.39) 
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and, as a consequence, the following differential equations are of odd parity: 


d2 

ar +21/3 =, (9.40) 
On, 3 dx\ 8 

we le 59) Ses 2 peas 

ete e(1 — |2]) (=) ; (9.40b) 
x gil 


9.3 TRULY NONLINEAR OSCILLATORS 


The dynamic systems modeled by Hamiltonians having the structure given 
in Equation (9.1) represent “truly nonlinear (TNL)” oscillators. 

To understand what these oscillators are, we begin with the following def- 
inition. 
Definition 9.3.1 A function f(x) isa TNL function, at x = 0, if and only if 
f(z) has no linear approximation in any neighborhood of « = 0. 


Explicit examples of TNL functions are 
file)=2¥8, foe) =a,  fa(x) = Ia, 
3 
fala) =a+a'3, f5(x) = pop. 


Definition 9.3.2 If f(x) isa TNL function and if solutions to the differential 
equation 


(9.41) 


dx 
Te + Fle) =, (9.42) 
are periodic, then Equation (9.42) is a TNL oscillator. 
Note that if H(x,y), in Equation (9.1), is taken as the Hamiltonian of 
a system, then its “solutions,” « = x(t) and y = y(t), are periodic, and 
the corresponding equations of motion, Equations (9.3), are a coupled set of 


nonlinear, first-order, ordinary differential equations. 
Using the transformation 


(x,y) = (=a =y), (9.43) 


in Equations (9.3), it is easy to see that these equations are of odd parity. 
For the remainder of this section, we will investigate exact solutions to 
several cases of Equation (9.1) for which p = 2, i.e., 


H(a,y)=y? + |2|4=1, 2(0)=1, y(0)=0. (9.44) 


The selected g-values correspond to dynamical systems of importance to var- 
ious phenomena in the natural and engineering sciences. 
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9.3.1 Antisymmetric, Constant Force Oscillator 


The relevant equation of motion for this case is 


2 
a + sgn(x) =0 (9.45) 
where, as a reminder, 
+1, «>0, 
sgn(z)= 40, «x=0, (9.46) 
-1l, «<0. 


Consequently, Equation (9.45) is equivalent to the following two linear differ- 
ential equations 


41 =0, z>0, (9.47a) 
1 =0, «<0, (9.47b) 
which have the respective solutions 
(== ($) P+ Att By (9.48a) 
te) = (5) + Aat + Bo (9.48b) 


where (Aj, Ag, Bi, Bo) are integration constants. 
To construct a solution, start with the initial condition 


dx(0) 


=> —— A - 9.49 
20)=0, “9 _ aso (9.49) 
Since A > 0, we must use 2‘+)(t) and impose the initial conditions given in 


Equation (9.49); doing this gives 


(H(0) = By =0 
x 1 5) 
dab (0 (9.50) 
= , = AL =A, 
and, therefore, 
a(t) = — (5) t(t-2A), O0<t<2A. (9.51) 


Inspection of Equation (9.51) shows that «'+)(t) has the following features: 


(0) = 
(t) > 
a(+)(2A) = 0. 


+) 0, 
+) 0, 0<t<2A, (9.52) 
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These results clearly imply that the period, 7, of the oscillation, is 


T =4A. (9.53) 
The conditions on x(t) at t = 2A and t = 4A are 
a'-)(2A)=0, a) (4A) =0, (9.54) 
and from Equation (9.48b), we find 
Az =-3A, Bo, =4A?, (9.55) 
and «‘—)(t) is 
a(t) = (5) #2? (3A)t+4A?, 2A<t<4A. (9.56) 


In summary, the periodic solution to Equation (9.45), with the initial condi- 
tions in Equation (9.49), is 


— ($) t(t — 2A), Us 6S 2A; 
One ies : (9.57) 
(5)? —(BA)E+4A?, 2A<t< 4A, 
with the periodicity condition 
a(t +4A) = a(t). (9.58) 


From this piecewise continuous representation for x(t), its Fourier series 
can be easily calculated and is given by the expression 


a= (=) ss ope “ae (9.59) 


k=0 


Observe first that only odd angular frequencies occur in the expansion, i.e., 
27 
wr = (2k +1) (3) = (2k + 1)wo, (9.60) 


where the period is T = 4A. Second, x(t) is an odd function of t. Third, the 
Fourier coefficients, b,, where 


16.A2 1 
b= ( 7 ) (2k + 1)3’ oP) 


decrease asymptotically as k~°. 
Finally, this system corresponds to the following H(z, y) 


1 
(ev) = (5) 9? + lel (9.62) 
which on comparison with Equation (9.1) and rescaling the y-variable give 


pH 2s. qa 1: (9.63) 
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9.3.2 Particle in a Box 


Many physical systems can be approximated as taking place in a one- 
dimensional box. To begin, assume that the box is located between x = 0 and 
x = L. Thus, the size of the box is L. Assume that a (classical) particle is 
located in the box and that at t = 0, it is at x = 0, with velocity vo, i-e., it is 
moving to the right. After a time T*, where 


Tt= (9.64) 


VO , 
the particle hits the wall at « = L, and reverses direction. Then, after a time 
T*, it hits the wall at x = 0, and reverses direction. Assuming all the collisions 
are elastic, the overall motion is periodic with period T’, where 


T =2T* = —. (9.65) 
In terms of x(t), this feature is given by the periodicity condition 

x(t+T) = a(t) (9.66) 
with a related expression for the velocity u(t) = da(t)/dt, i-e., 

v(t + T) = v(t). (9.67) 


The time dependencies of x(t) and v(t) are 


t 0<t<f 
a(t) = te poe (9.68) 


w(T—-t), £<t<T, 
and 
vw, O<t<4, 
t) = 9.69 
“(0 ea T<e<v. ee) 


A rescaling of the z(t) and v(t) variables to units such that v9 = 1 and L =7, 
simplifies the Fourier series representations for x(t) and v(t). With these units, 
we have the (rescaled) period, T = 27, and the expressions 


T A\  cos(2k — le 
4 sin(2k — 1)t 
= »-(93 1) Qk 21) (9.70b) 


Observe the following: 
(i) The Fourier series contains only odd harmonics. 


ii) Equation (9.70a) gives the following well-known expression for 7? 
& 


Co 


7 1 ow 1 
T 8), pap 8 Gera (9.71) 
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9.3.3 Restricted Duffing Equation 
The restricted Duffing equation 


aa 


3 
is a special case of the general Duffing equation 
Cx dx 2 3 
Gp tty toate +dx° =0, (9.73) 
where 
a=0, b=0, c=0, d=1. (9.74) 
For the initial conditions 
dx(0 
x(0) = A, a ) ag (9.75) 


the exact solution can be expressed in terms of the Jacobi cosine elliptic 
function, cn(At; 1/2), ie., 


a(t) = A en(At;1/V2). (9.76) 
The Fourier series for “cn” is rather complex. Let (k, k’) satisfy the relation 
(k)* + (A")? = 1, (9.77) 
and define g(k) as 
TE (k’) 
k) =e — : 
a(t) =e {- |] h. (9.78) 
where F'(k) is the complete elliptical integral of the first kind, i.e., 
m/2 d 
F(k) = ‘ a (9.79) 
0 1—k? sin? 


With these definitions, the general Jacobi elliptic function has the Fourier 
series representation 


7 = m+3 
cn(u, k) = Ee » (| cos(2m + 1)u, (9.80) 
where 
v = v(k,u) = Spa U. (9.81) 
2F(k) 


For the restricted Duffing equation, we have 


k 


x Sil, cee 
=e F (=) 1.854074..., nas 
q 


(+) = 0.043213... 


Using the results given in Equations (9.77) to (9.82), the Fourier series can be 
constructed. Again, note that only odd harmonics appear. 
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9.4 ATEB PERIODIC FUNCTIONS 


Consider the following coupled, first-order, nonlinear differential equations 


with initial conditions 

x(0)=1, y(0)=0. (9.84) 
In the «-y phase-plane, there is a fixed point at (%,y) = (0,0). Also, x(¢) 
satisfies the following second-order differential equation 


d?x(t) 2 
t)* =0. ; 
qe +(S) ae 0 (9.85) 
If this equation is to be of odd parity, then a must either have the form 
— 2n+1 
- ee (9.86) 
W=051,2,3,2643- m= 051,253 346, 


or the term x(t)® replaced by 
x* — |x|“sgn(x). (9.87) 


The trajectories of the curves in the x-y phase-space are solutions to the 
first-order differential equation 


dy at+l1 oe 2n+1 

oe Gru cy) Bere ae 
or 

d 1 a 

sf LY IPN ye pce (9.88b) 

dx 2 y 


A first integral can be found by rewriting Equation (9.88) to the forms 


dy + (24) x%dx = 0, q= tt. 
ydy + (S*) ante (9.89) 
ydy + (SF) |a|*sgn(x)dx =0, a= real > 0, 
and integrating to obtain 
Ge a. ($) |z|°*+ = constant, a =real > 0 
2 2 = , — < 
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Imposing the initial conditions, 7(0) = 1 and y(0) = 0, gives 


Tia goth 1 — 2n+1 

2 °° 2 — 2 OS 2m+1? (9 91) 
: 4 

wey (3) ier = 4, a= real > 0. 


Comparing these expressions with Equation (9.1) gives p= 2 andg=a-+l. 

In an important paper by Rosenberg, published in 1963, he associated the 
name “Ateb functions” with x(t) and y(t), and examined the most important 
of their mathematical properties. In our notation, we will call z(t) and y(t), 
respectively, the Ateb-cosine and Ateb-sine functions, i.e., 


At th=alt 
ola, t) = a(t), iiss 
Ats(a,t) = y(t). 
Let II(a@) be defined as 
1 1 
II(a) = B | ——, = 9.93 
(@)=8(—.5). (9.93) 
where B(a,b) is the beta function, ie., 
T'(a)P(6) 
Bia, b) = ———=- 9.94 
(0.0) = Te (9.94) 
then the following properties can be easily established: 
(i) Atc(t) and Ats(t) are, respectively, even and odd, ice., 
Atc(—t) = Atc(t), Ats(—t) = —Ats(t). (9.95) 
(ii) Both Atc(t) and Ats(t) are periodic, with period 2II(q), i.e., 
Atc(t + 210 = Atc(t 
oft + 211(a)) = Ate(t), a 
Ats(t + 2II(a)) = Atc(t). 


(iii) Atc(t) and Ats(t) satisfy the following coupled, first-order pair of differ- 
ential equations 


Catett) = - (=) Ais (9.97a) 
< Ats(t) = (Ate(t))®. (9.97b) 


(iv) Atc(t) = x(t) is a solution to the following nonlinear, second-order dif- 
ferential equation 


dx 2 a 2n+1 
et(saz)e es oO Om +1’ Cee) 


d? 2 
ait ( -) lel*sgn(2) = 0, a=real>0. — (9.98b) 
= 
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It should be indicated that the periodic Ateb functions have found a signif- 
icant use in the analysis of certain nonlinear, conservative oscillatory systems. 


9.5 EXACT DISCRETIZATION OF THE JACOBI 
ELLIPTIC DIFFERENTIAL EQUATIONS 


The Duffing differential equation provides a nontrivial mathematical model 
for a broad range of phenomena in the natural and engineering sciences. It also 
corresponds to the simplest generalization of the so-called simple harmonic 
oscillator equation 


d?y 
— =0. 9.99 
arty (9.99) 
For our purposes, we consider the Duffing equation in the form 
d*y 
meow. by? = l 
72 + ay + by” = 0, (9.100) 


where (a,b) are real parameters. The following initial conditions and parame- 
ter signs give solutions that are, respectively, related to the Jacobi cosine and 
sine elliptic functions: 


y(0) = A, ato) =0; a>0, b>0; y(t) = Acn(t) (9.101a) 
y(0) =O, at) HA ase b<0> Ysa,  “1018) 


The main goal of this section is to construct discretizations of the Duffing 
differential equation such that at the discrete-time lattice points, the solutions 
of the difference equations are exactly equal to the solutions of the differential 
equations. This effort can play an important role in the numerical integration 
of the equations of motion for dynamical systems modeled by the Duffing 
equation. 


9.5.1 Rescaled Duffing Equation 


The Duffing differential equation can be rescaled using the substitutions 
y(t) = Ax(t), t=Tt, (9.102) 
where (A,T) are the scaling parameters, and x and # are the new dependent 
and independent variables. Substitution of this result into Equation (9.100) 
gives 
Px 2 2 42,3 
In the next sections, we show how to select A and T. 
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9.5.2 Exact Difference Equation for CN 


We now determine the exact difference equation for the Jacobi cosine func- 
tion cn(t,k). In the calculations to follow, the parameter k will not be indi- 
cated. 

To begin, we make use of an addition theorem for cn/(t). 


(cn u)(cnv) — (snu)(dnu)(snv)(dnv) . 


cn(u+v) = 1—Flsnu)"env? (9.104) 
Now define 7 7 
oe h=At, wn =cn(ty), (9.105) 
N =0,1,2,..., 
and use Equation (9.104) to obtain 
ear ere os (9.106) 
To obtain this latter result, we made use of the relation 
[sn(tw)]? = 1—[en(tn)]? =1—-2%,. (9.107) 


Adding (—22y) to both sides of Equation (9.106), and replacing (cn h)xy by 


(cnh)azn = an — [1 —-(cnh)an, (9.108) 
we obtain 
UN41 — 2tn+2tn-1 ip { - = a ee pos + aa \ 
(sn h)? (sn h)? 2 
4 2k? pt) 2, =0. (9.109) 


This second-order, nonlinear, difference equation has 
tn =cn(tn) = cn(ty,k) (9.110) 


as its solution. 
The differential equation corresponding to this difference equation is de- 
termined by taking the following limits 


h-+0, N->co, £=hN = fixed. (9.111) 


Using also the results 


ee 9.112) 


enh =1- © + 0(h'), 
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Equation (9.109) becomes in these limits the differential equation 


dx 2 2.3 

Gp + (1 — 2h? )a + 28708 = 0. (9.113) 
Comparing this with Equation (9.103) and equating the similar coefficients of 
x and x°, gives the relations 


lear, 2 Sar A’, (9.114) 
If these equations are now solved for T? and k?, we obtain 


1 bA? 
2 2 
T =e k = 30a + 6A)’ (9.115) 


where from the second expression, we find that k? satisfies the bounds 
1 
0<k< a (9.116) 


The time scale T can be interpreted as the “physical time scale” for a 
system modeled in the form of the Duffing equation. 
Finally, it should be indicated that 


x(tn) =n, (9.117) 


where x(t) and xy are, respectively, the solutions to Equation (9.103) and 
(9.109), with T? and k? given in Equation (9.115). 


9.5.3 Exact Difference Equation for SN 


The methodology of Section 9.5.2 can also be applied to the derivation of 
an exact discretization for the Jacobi sine function. We do not provide the 
details, but give a summary of the procedure. 

The sn addition theorem 


(sn.u)(cnv)(dnv) + (snv)(cnu)(dn wu) 


sn(ut+v) = 1 — k2(snu)2(snv)2 


(9.118) 


along with the relation 
(dnu)? =1—k?(snu)?, (9.119) 


can be used to derive the following second-order, nonlinear difference equation 
for xy = sn(ty, k) 


tnN41 —2an +2N-1 ee 
(enh? (sn h)2 7 


— 2k? (este) x, =0, (9.120) 
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and the limiting differential equation is 
d2 
oF 4 (14k?) — 2k223 = 0. (9.121) 
dt? 
Likewise, for this case, where a > 0 and b < 0, we have 


2 2 2 |b|.A? 


~ Qa — |b A2’ ~ Qa — |b) A2” one) 


Since we must have 0 < k? < 1, it follows that the amplitude is restricted in 


magnitude, i.e., 
|A] < Vi (9.123) 


9.6 HARMONIC BALANCE: DIRECT METHOD 
9.6.1 Methodology 


Harmonic balance procedures are a set of techniques which allow the cal- 
culation of analytical approximations to the periodic solutions of differential 
equations. These methods are based on representing a given periodic solution 
by means of a truncated Fourier series, i.e., a trigonometric polynomial. We 
now outline this procedure. 

For our purposes, assume that the differential equation under consideration 
can be written in the form 


dx d?x 


Further assume that it is of odd parity, i.e., 


dx d’x dx d?x 
F (-«. 2. -5) a, («. an =) (9.125) 


and it has a periodic solution which can be expressed as a Fourier series, i.e., 


x(t) = S {Ax cos[(2k — 1) M4] + By sin[(2k — 1)Q4]}}. (9.126) 
k=1 


Let z(t) be the N-th-order harmonic balance approximation to «(t). It cor- 
responds to the representation 


N 
y(t) = S> {Ap cos(2k — 1)Qnt + By sin(2k — 1)Qn(t)}, (9.127) 
k=1 
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where (AX, Br Quy) are approximations to (Az, Br, ) for k = 1,2,..., N. 
In general, for conservative, one-space-dimensional systems, Equation 
(9.124) takes the form 


d?x 
qe f(x), (9.128) 
and the initial conditions can be selected to be 
d 
(0) = A, nO) 5G. (9.129) 


with the consequence that only cosine terms appear in the Fourier expansion, 
Le., 


a(t) = S~ Ax cos(2k — 1)9t, (9.130) 
k=1 
and es 
y(t) = S> Ay’ cos(2k — 1)Ont. (9.131) 
k=1 


In the work to follow, only conservative dynamical systems will be discussed. 
From Equation (9.1), this situation corresponds to p = 2. 

Note that x(t) has (NV-+1) unknowns: the N-coefficients, (AY, AY,..., AN) 
and Qy. These quantities can be determined by performing the following 
steps: 

(i) Substitute Equation (9.131) into Equation (9.128) and expand into 
harmonics, such that it has the form 


x Hy, cos|(2k — 1)Ont] + HOH ~ 0. (9.132) 
k=1 
The H;, will be functions of the N-coefficients and Qn, i-e., 
Hy, = Hy (AY, Ad’,..., AN, Qn). (9.133) 
(ii) Now equate the N-functions to zero, i.e., 
A,=0, k=1,2,...,N. (9.134) 
(iii) Solve these N-equations for Ad’, A’,... AN, and Quy, in terms of AN. 


(iv) Use the initial conditions, see Equation (9.129), to determine AN, i.e., 


N 
ay(0)=A=AN+5— Ap (AD). (9.135) 


We now illustrate the use of the harmonic balance method by applying it to 
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two systems for which their Hamiltonians have the structure given in Equation 
(9.1). For all examples, we have p = 2. We will also use the following compact 
notation for time derivatives 


(9.136) 


Finally, it should be pointed out that the same initial conditions are im- 
posed on each zy (t), ie, en (0) = A and &y(0) = 0. 


9.6.2 #+2°=0 
The equation of motion 
#+2°=0, 2(0)=1, 2#(0)=0, (9.137) 


has the Hamiltonian function 
tae = (3) 9? a (3) t= i (9.138) 
a )=1, 


where y(t) = dx(t)/dt = «(t). Note that with a rescaling of the variables, 
H(a,y) corresponds to Equation (9.1) with p = 2 and gq = 4. However, for our 
calculations, we will keep the results in Equation (9.138). 

The N = 1 approximation, x1(t), is 


x1(t) = AcosMit, (9.139) 


where, from now on, we drop the “bar” on coefficients and the angular fre- 
quency Q. 
Substituting Equation (9.139) into Equation (9.137) gives (0 = Qyt) 


(—AQ7) cos 6 + (Acos 6)? ~ 0 (9.140) 
and 
A 92 + (3) a cos0 ++ HOH ~ 0. (9.141) 
Therefore, 
= (3) A’, (9.142) 


and the first-order harmonic balance solution is 
3\ 1/2 
— At|. 
() 


21(0)=1 and #,(0) =0, (9.144) 


x1(t) = Acos (9.143) 


With the initial conditions 
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we have 
3\ 1/2 
x1(t) = cos (3) t| . (9.145) 
The next approximation, x2(t), has the form 
X2(t) = Ay cos@ + Ag cos30, 8 = Not. (9.146) 


If this is substituted into Equation (9.137) and the necessary trigonometric 
expansions are done, then the following result is obtained 


H, cos@ + Hz cos30+ HOH ~ 0, (9.147) 
where 
H, = A, & = (3) A? — (3) Ay As — (3) ay (9.148) 
2 l\ 42 3) 42 3) 43 

Ay = —9A205 + Zz Aj + 3 Aj Ao + 7 As. (9.149) 

If z is defined as 
a (9.150) 

—— 3 
Ay’ 


and if Hj, is set to zero, we obtain 


3 1/2 
% = (3) Ail 24 927)*, (9.151) 


Now setting H2 to zero and using 22 above, we get the following cubic equation 
for z 
512° + 2727 + 21z-1=0. (9.152) 


While this equation has three roots, the one we seek must be real and “small,” 
ie., z is the ratio of the amplitudes of the third harmonic, cos(36@), to that of 
the first harmonic, cos @. We expect this ratio to be small if this method is to 
be valid, ice., 

lz] <1. 9.153) 


If we denote this root by 21, then 
zy = 0.044818.... 9.154) 


Therefore, x2(t) is 


X2(t) = Ay[cos 6 + 21 cos 36]. 9.155) 


Using the initial condition, x2(0) = 1, gives 
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or 7 
A, = ——.,, (9.157) 
1+ ZA 
and fa ; 
3 (1+ 21 + 223)1/2 
Q2 = | - ——_——_——_=+— | = 0.8489. 9.158 
az @ | +2 O28) 
The periods T = 27/Q, for the first and second approximations, are 
2 2 
tT, =a 70554, TS ST ANIG, (9.159) 
Qy Q2 


and these are to be compared to the “exact” result 
Texact = 7.4163. (9.160) 
The corresponding percentage errors, defined as 


Texact aa 


E = percentage error = | | - 100, (9.161) 


exact 


are 
E, =2.2% and Ep =0.20%. (9.162) 


In summary, the harmonic balance method provides a quick, efficient tech- 
nique for calculating analytical approximations for the periodic solutions of 
#4+.2° = 0. In essence, the procedure converts the solving of a strongly nonlin- 
ear differential equation to one where nonlinear algebraic equations must be 
solved. Further application of this technique gives even more accurate results 
for the period. In addition, the ratio of neighboring amplitude coefficients 
remains small, a result needed to insure the validity of the method. 


9.6.3 ¢+27/3=0 


The “cube-root” oscillator 
#+a3=0, 2«(0)=1, (0) =0, (9.163) 


has the Hamiltonian 
1 3 3 
H =(-)y,? ~)e*2 == 9.164 
(4) (5) +(3)2 7 (9.164) 


where y(t) = «(t). 
The exact angular frequency can be calculated and is found to be 


Qexact = 1.070451. (9.165) 
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There are several ways to proceed to determine first- and second-order har- 
monic balance approximations to the periodic solutions of Equation (9.163). 
We start with a cubed version of this equation, i.e., 


(#)? +2 =0. (9.166) 


For 
x1(t) = Ay cos(Qit) = A; cos, (9.167) 


we have on substitution into Equation (9.166) 


[—(1)?A; cos 6] + Ay cos4 ~ 0, (9.168) 
and ; 
Ay f - (3) Aa)" cos0 + HOH ~ 0. (9.169) 
This gives for Q1 the value 
Ta a 
Q1(A1) = (5) (=) : (9.170) 
1 
Since 21(0) = 1, it follows that 
4\ 1/6 
A; =1 > 2;(t) =cos (2) | ; (9.171) 
and 
4 (1) = 1.049115. (9.172) 


Therefore, the percentage error in the angular frequency, 1);, is 


- 100 = 2.0%. (9.173) 


If we write x2(@) as 


(9.174) 


x2(t) = Ay[cos@ + zco0s36], 8 =Qat 
2(0)=1, 4#(0) =0, 


substitute into Equation (9.166), and calculate the coefficients of the cos @ and 
cos 38 terms, then the following expressions are obtained 


(Q2)° A? (3) - (2) z+ (=) 2| =, (9.175a) 
(Q2)° A? (+) + (2) z+ (3) ‘| = 2: (9.175b) 
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Dividing these equations gives 
(1701)z* — (27)z? + (51)z +1=0, (9.176) 


and the smallest (in magnitude) root is 


z, = —0.019178. 9.177) 
Also, from #2(0) = 1, 
Ai = = 9.178) 
and (leaving out a lot of the details), we calculate 
Q2 = 1.063410. 9.179) 
The percentage error 
Ey = ag - 100 = 0.7%. (9.180) 


In summary, by rewriting the equation of motion a rather accurate value 
has been determined for the angular frequency (and by extension, the period). 
Note that the ratio of the amplitudes of the first two harmonics is small. 


9.7 HARMONIC BALANCE: RATIONAL APPROXI- 
MATION 


9.7.1 Methodology 


The standard harmonic balance method is based on constructing truncated 
approximations to the Fourier series of periodic x(t) and, as a consequence, 
contains only a finite number of harmonics. An alternative technique is to 
start with a rational form for x(t), such as 

A; cos @ 
tr(t) = —— ———,, _ 0 = 1, 9.181 
r(t) 1+ B, cos 20 ; ( ) 
where (A;, B,, 91) are to be determined. The general initial conditions to be 
satisfied are taken to be 


xrR(0) = A, £R(0) =0. (9.182) 


Note that this approximation for x(t) contains all the odd harmonics. In fact, 
it can be demonstrated that 


p(t) = 5° axcos(3k+1)0, @=Mt, (9.183) 
k=0 
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where 
k 
(-1)FA,-|. (t= x/PSB? | Bp tse Be? 
Gp ee | | ee | | | ee... O1ee 
V1— 3B? By By 
Observe that this expression and that in Equation (9.182) only makes sense if 
|Bi| <1. (9.185) 
If we assume 
|Bi|<1, (9.186) 


then the coefficients ax, decrease exponentially with k, i.e., let “a” be defined 


as ‘, 
e* = rae (9.187) 
2 
then 


-1)FAye"**, 0< By <1 
m4 eee Mpeg (9.188) 


Aje~%, 0< (-—B,) <1. 


For purposes of the harmonic balance, ar(t) must be calculated. A 
straightforward, but lengthy calculation gives 


(1+ B, cos 20)?#p = -«e#.{ a + Bi)- (=) Bi | cos 0 


B 
+ 3B, a = 1 cos 30 + nou}, (9.189) 


We now illustrate the use of this harmonic balance procedure by applying 
it to two examples. 


9.7.2 ¢+2°=0 


To begin, if we rationalize Equation (9.181) and the cube of the resulting 
expression, we obtain 


(1 + B, cos 20)?x?, = (A; cos 6)? 
A$ A? 
= (2) cos 6 + (+) cos 30. (9.190) 


This result and that from Equation (9.189) allow us to write 


é+a° =0, (9.191) 
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as 
{-(o2) a + By) — (5) Bi | a al cos 0 


+ ¢ —(07A;)(3B,) Pig gee cos 30 
ia 


+ HOH ~ 0. (9.192) 


Setting the coefficients of cos @ and cos 30 to zero gives the two expressions 


0? ja i oe ($) a = (3) A2, (9.193a) 
07(3B1) a - 1 = a (9.193b) 


Dividing these two equations gives a relationship involving only By, i.e., 


(2) B? —10B, -1=0. (9.194) 
Solving for B, gives, for the root having the smallest magnitude, the value 

B, = —0.090064. (9.195) 
Now, solving Equation (9.193b) for 02, i.e., 


2 (3) At 


"= Tae) - ay saa 
Using the above calculated value of B, gives for the Q; the value 
Q) = (0.930982) A). (9.197) 
If we take the initial conditions to be 
x(0)=1, «(0) =y(0) =0, (9.198) 
then from Equation (9.181) we obtain 
Aj =1+8h, 9.199) 
and 
Q, = 0.847134, 9.200) 


with the period 
2 
eS = Ar: 9.201) 
Qy 
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The exact period for this example, with the initial conditions stated in Equa- 
tion (9.198), is 
Texact = 7.4163, (9.202) 


and, therefore, the percentage error in our calculated period is 


Texact —T 
Texact 


| - 100 = 0.01%, (9.203) 


which, for most practical purposes, is very small. 
To this first-order of approximation, the rational harmonic balance proce- 
dure gives, for # + x? = 0, the following result 


(0.909936) cos[(0.847134)E] 
t) = ee _. 9.204 
@R() = 9 990064) cos|(1.694268)F] ee!) 
9.7.3 4+ 27sgn(x) =0 
The quadratic nonlinear oscillator can be written in two equivalent ways 
&+|r\x =0, (9.205a) 
# + |a|?sgn(0) = 0. (9.205b) 


We will construct an approximation to the periodic solution using the rational 
harmonic balance technique. Writing this approximation as 


A; cos@ 


SS OF 9.206 
1 + Bi cos 20’ a ( ) 


tp(t) 
our task is to determine A;, B,, and Q1, subject to the initial conditions 
tr(0)=1, &r(0) = yr(0) = 0. (9.207) 


Note that a first integral or Hamiltonian for this system is 


1 1 1 
H =(=)y?4+(=)|eP == 2 

a= (3)? +($) l= 5, (9.208) 

which, on comparing with Equation (9.1) gives p = 2 and q = 3. 

Assuming 
0< |Bi| <1, (9.209) 
then Allee 
1|cos 

t)| = —————— 9.210 
len (t)| 1+ B, cos 20 ( ) 


where the |cos @| has the Fourier expansion 


jot = (4) (2) +()owre(L) coors J. oan 
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Therefore, 


q 3B 
(1+ B, cos 20)? |rr|r2R = (=) i+ cos 6 


cs (=) E & we cosan + ont) (9.212) 


5 


Substituting Equations (9.190) and (9.212) into Equation (9.205a) and equat- 
ing to zero the coefficients of the resulting expressions in cos @ and cos 380, gives 
the following two expressions to be solved for B, and 91, in terms of Aj: 


A? B 11 
cael eae = OFA, |(1+ Bi) -—( — ) Bil, (9.213) 
37 i) 2 
8A? 17B, 2 9 2 
(=) i+ 7 = O{A1 1 By -— 3B). (9.214) 
Dividing these two equations, with simplification of the resulting expression 
gives 
563 169 129 
— | Bi+(—)Bi-(—)]B-1= Bal 
Cale ae Ge a ea 


for which the smallest-magnitude real solution is 
By, = —0.052609. 9.216) 
We also have, from xr(0) = 1, 
A=1+8,, 9.217) 
and can use these to calculate 
Qy = 0.95272, 9.218) 


and (0.9474) cos[(0.9527)t] 


t) = —— 9.219 
tR(Y) = 719.9526) cos[(1.9054)F ) 
The exact period for this problem is 
Qa 1 1\]° 
Texac a Sh T(- De . 9.220 
Hee). ee 


9.8 ITERATION METHODS 


Nonlinear differential equations modeling periodic systems may also be 
solved by constructing iteration schemes to obtain approximations to their 
solutions. We now provide an outline as to how this can be done. 
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First, we formulate a direct iteration procedure. The general idea is that 
the original nonlinear differential equation is transformed into an infinite sys- 
tem of linear inhomogeneous equations such that at a particular step of the 
calculation, knowledge of the solutions of the previous members of the se- 
quence determines the inhomogeneous term. The basic steps to accomplish 
these goals are now listed: 

(a) Assume the differential equation 


F(#,2)=0, 2(0)=A, #(0) =0, (9.221) 


is of odd parity and assume that it can be rewritten to a form (which may 
not be unique) 
&+ f(%,x) =0. (9.222) 


(b) Next add 0?z to both sides of Equation (9.222) and obtain 
#+ 0x = 072 — f(#,2) = G(4, 2). (9.223) 


Note that, at this point, 0? is not known. 
(c) The iteration scheme is given by 


Brsr + Ojter1 = G(Fx, ek), (9.224) 
RSA Dies 
xo(t) = Acos(Qot), (9.225) 
where the x,41(t) satisfy the initial conditions 
K41(0) = A, &x41(0) = 0. (9.226) 


(d) At each step of the iteration, Q, is determined by the requirement 
that secular terms do not appear in the solution for 7,41(t), ie., ep41(t) is 
periodic and bounded. 

Note that at step-(&k +1), Q, appears in Equation (9.224) and, further, at 
the next step of the iteration (k +2), 0,41 appears, but is unknown and must 
be calculated. The worked example will help to clarify this point. 


Comments: 

(1) The inhomogeneous term in the linear differential equation for 7,+1(t), 
depends only on the solutions for k less than k + 1, and consequently, are 
known. 

(2) The linear differential equation for x,41(t) allows the determination of 
QO, by requiring that secular terms be absent. Therefore, the angular frequency, 
“Q,” appearing in %, and 2, in the inhomogeneous term for Equation (9.224) 
is Ox. 

(3) Since, Equation (9.222) is of odd parity, the calculated x, will only 
contain odd multiples of the angular frequency. 
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9.8.1 Direct Iteration Scheme: # + x° = 0 
A possible direct iteration scheme for 
#+a° =0, (9.227) 
is 
Fee + Ree = O22, — 23. (9.228) 


Since, xo(t) = Acos(Qot), where Qo, for the time being, is not known, Equa- 
tion (9.228), for k = 0, is 


#1 + 0221 = O2a9 — x8 = N2(Acos 6) — (Acos6)* 


3 A’ 
— ng _ (3) a Acos 6 — (=) cos 30, (9.229) 


where 6 = Qot. Secular terms will not occur in the solution for x(t) if the 
coefficient of the cos@ term is zero, i.e., 


92 — (3) A? =0, (9.230) 


Q% = [oA (9.231) 


Therefore, Equation (9.229) becomes 


which applies 


As 
#, +0071 = - (=) cos 36. (9.232) 


The general solution of this last equation is 
a, = af) 4 2), (9.233) 
where 


(9.234) 


and of) and a? ) denote, respectively, the homogeneous and particular solu- 


tions, with C an integration constant. Since 71(0) = A, it follows that 


23 
C (=) (9.235) 
and 2, is 
23 1 
ai(t)=A (3) cos 6 + (=) cos 20 ; (9.236a) 


6=Not = [240 (9.236b) 
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Thus, the first iteration approximation to the periodic solution of #+ 2x? = 0, 
x(0) = A, and «(0) = 0 is 


23 3 1 3 
= jad ; At — ; —At]}|. 2 
ay(t)=A (=) cos (2 + (=) cos (23 (9.237) 
To extend the calculation to k = 1, we use 
Bo + 0225 = O%0) — 23, (9.238) 


where 2 (t) takes the form of Equation (9.236a), but now @ is equal to it, 


iLe., 
r=A (3) cos(Qyt) + @ cos(30h)| 


AA (=) eee (=) cos(3) (9.239) 


If this expression for x1 is substituted into Equation (9.238), if the resulting 
coefficient of cos@ is set to zero, and if the complete solution to x2 is then 
calculated, the following results are obtained, after a long series of algebraic 
manipulations: 


v(t) = A{(0.955) cos + (4.29) - 107? cos 36 
+ (1.73) - 107% cos 54 + (3.60) - 10~° cos 70 


+ (3.13) - 1077 cos 96}, (9.240a) 
6 = 04(A)t = (0.849326) At. (9.240b) 


Observe that from 


(A) = (0.866025) A, 
Q,(A) = (0.8493276).A, (9.241) 
exact (A) — (0.847213) A, 


we calculate the corresponding percentage errors to be 


exact (A) = Qo (A) a. 

Pe 100 = 2.2%, (9.242a) 
exact (A) = Qo (A) 
ana = ET ee ae 242 
| Qexact(A) 02% (9.242b) 


Our conclusion is that the second-order iteration gives an order of magnitude 
decrease in the percentage error for the angular frequency. 
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9.8.2 Extended Iteration: *+ 2° = 0 


The direct iteration scheme, given by Equation (9.224), can be generalized 
to the following extended iteration formulation 


2 = 
¥ep1 + OZCe41 = G(Epe-1,2r-1) 


+ Gy (€p-1, Tk-1)(@k — Le-1) 


+ Ge (@e—1, Te-1) (Ee — Le-1), (9.243) 
where k = 0,1,2,..., and 
G,= oe G;= ae (9.244) 
Ox zt 


and the x4, satisfy the initial conditions 
K41(0) = A, &x41(0) = 0. (9.245) 


Inspection of Equation (9.243) indicates that this iteration method requires a 
knowledge of x_1(t) and xo(t), and these are taken to be 


x_(t) = xo(t) = Acos(Qot). (9.246) 


However, a computationally better scheme to use is the following modification 
of Equation (9.243). 


fnp1 + Urns = G(0, to) + Go(#o,20)(e — to) 
+ Ge(Zo, 20) (Ex — 2). (9.247) 


Therefore, for # + x? = 0, the extended iteration scheme, as represented 
by Equation (9.247), is 


Bro + %@e41 = G(x) + Gz (xo) (Te — Xo) 
= (N22 — x3) + (N2 — 3x2) (aE — 20), (9.248) 


where, and this is important, 
xo(t) = AcosO, 0 = Oxt. (9.249) 
For k = 1, we have 
fo + 2x2 = (N29 — £3) + (0? — 3a2)(x1 — 20), (9.250) 


where 
xo(t) = cos, 


az1(t) = A [ (24) cosd + (4) cos 36] , (9.251) 
6 = Mit. 
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A direct, but lengthy calculation gives 


((A) = (0.846990) A, (9.252) 
13, 244 595 
L(t) = a4 (S ==) cos 6 — 5 =) cos 36 
23 
——_ 356 >. 12 
+(sas) cos 5 ! (9.253) 


Two comments need to be made. First, the fractional error for (A) in Equa- 
tion (9.252) is very small, i.e., 


Qexac or Q 
an - 100 = 0.03%. (9.254) 


exact 


Second, the ratio of the neighboring coefficients is of order 10~?, ice., 


|S] = ean a(S) lr @, 
lee 


(9.255) 


ae = (3.9) -107?, 


and these results suggest that the extended iteration scheme may have higher 
harmonic coefficients, which rapidly decrease to zero. Clearly, the extended 
scheme is more accurate than its direct version. 


9.9 DISCUSSION 


It should be noted that for all the applications considered in this chapter, 
the equation of motions have the form 


aa 


qt alz|“sgn(x) = 0, (9.256) 


where the parameter a is positive, and a > 0. Written as a system of coupled, 
first-order, differential equations, we have 


dx dy i 
ape ae —a|x|“sgn(z). (9.257) 
We also used the following initial conditions 
_ dx(0) 
dt 


x(0)=1, (0) =0, (9.258) 
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and this corresponds to the Hamiltonian function 


He) = +( 


; ) jzjot? = (9.259) 


atl atl 


An interesting and very important feature of the solutions to Equation 
(9.256) is that all are periodic, and the exact value of the period is 


4/yT+ al (5) 
ee ee 


a) (9.260) 

a -ayr [53] 
Further, this differential equation is of odd parity. This implies that with the 
initial conditions stated in Equation (9.258), the Fourier series for x(t) has 


the form 
co 


a(t) = S> ax cos[(2k + 1)4], (9.261) 
k=0 
where 2 = 27/T. 

A comparison of Equations (9.1) and (9.259) shows that p = 2 and q = 
a+1. (Note that renormalizing the variables x and y, Equation (9.259) can 
always be transformed into the expression of Equation (9.1).) The reason why 
p = 2 is because the systems investigated in this chapter model phenomena 
appearing in the physical universe and, for such systems, the Newton force 
law equation is the starting point. Since the Newton force law takes the form 


(9.262) 


Force = (mass) acceleration, 
F=ma, 


then, for a point particle in one-space dimension, the acceleration can be 
expressed as 


a=—5 (9.263) 
Defining the momentum as 
dx 
momentum = ma =P (9.264) 
and using units of mass such that m = 1, then 
da(t 
p(t) = 28 <n, (9.268) 


Consequently, in the «-y phase-space, 


d?x = dy dy 
aay ig re (9.266) 
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If the system is conservative, then the force depends only on 2, i.e., 


Force = F(x), (9.267) 
and it follows that d 
F(x) =y—. (9.268) 
dx 


Integrating this expression gives 


ve + U(x) = constant (9.269) 
U(a) = — f F(az)dz. 
If U(x) is selected to be 
alx|ott 
= 9.270 
U@) =m (9.270) 


then with the initial conditions given in Equation (9.258), we arrive at the 
result of Equation (9.259), and comparing with Equation (9.1), we find 


p=2, q=atl. (9.271) 


A detailed examination of the Hamiltonian function, Equation (9.259), 
shows that it corresponds to a closed curve in the z-y phase-space if 


a>-tl. (9.272) 
Thus, a system modeled by the differential equation 


dx 1 


ta) (9.273) 


has bounded phase-space curves and all solutions are periodic; observe that 
the Hamiltonian function 


fae =£ + (3) (3 = 3 (9.274) 


However, a = —1 requires a very detailed and careful analysis to prevent 
difficulties from arising. 


Another way of characterizing Equation (9.256) is to observe that it is 
a “truly nonlinear (TNL)” oscillator differential equation; see Section 9.3. 
Equations of this type have the feature that no “small parameter” can be 
introduced such that any of the standard perturbation techniques may be ap- 
plied to calculate analytical approximations to their periodic solutions. There 
are basically two ways to get around this limitation. First, we can numeri- 
cally integrate the equation of motion, Equation (9.257). Second, schemes can 
be created which allow the determination of analytical approximations. Two 
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particular useful procedures are the method of harmonic balance and itera- 
tion techniques. Sections 9.6 and 9.7, and 9.8 respectively, deal with these 
two methodologies. With regard to harmonic balance methods, generally, an 
increase in the number of Fourier harmonic terms increases the precision of 
the calculated x(t) and the values of their periods. However, iteration meth- 
ods face some potential difficulties. The main problem being the selection of 
a scheme to be iterated. This is currently an active area of research. 

Finally, there are the so-called Ateb functions corresponding to the ex- 
act solutions of Equation (9.256). While there has been some work done on 
constructing analytical approximations to the Ateb functions, it is not en- 
tirely clear (to me) that this provides generally a better approach than the 
harmonic balance methods. Additional investigations are required to give a 
definitive resolution of these issues. 


NOTES AND REFERENCES 


Section 9.1: Excellent introductions to nonlinear oscillators are given in 
the following three books: 


(1) R. E. Mickens, Introduction to Nonlinear Oscillations (Cambridge Uni- 
versity Press, New York, 1981). 


(2) I. Kovacic and M. J. Brennan, The Duffing Equation: Nonlinear Oscil- 
lators and their Behaviors (Wiley, New York, 2011). 


(3) L. Cveticanin, Strong Nonlinear Oscillators: Analytical Solutions 
(Springer, Berlin, 2017). 


Section 9.2: This section is based on the paper 


(1) R. E. Mickens, Fourier representations for periodic solutions of odd- 
parity systems, Journal of Sound and Vibration 258 (2002), 398-401. 


Section 9.3: A full discussion of the concept of “truly nonlinear (TNL)” 
oscillators is the book 


(1) R. E. Mickens, Truly Nonlinear Oscillations: Harmonic Balance, Pa- 
rameter Expansions, Iteration, and Averaging Methods (World Scientific, 
Singapore, 2010). 


Section 9.4: The classical paper on Ateb periodic functions is 


(1) R. M. Rosenberg, The Ateb(h)-functions and their properties, Quarterly 
of Applied Mathematics 21 (1963), 37-47. 
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See also the book by Cveticanin (Section 9.1, above), Appendices A and B, 
and the following papers: 


(2) I. V. Andrianov, V. I. Olevski, and Y. B. Olveska, Analytical approx- 
imation of periodic Ateb-functions via elementary functions, Proceed- 
ings of the 5th International Conference on Nonlinear Dynamics (ND- 
KhPI2016, September 27-30, 2016, Kharkov, Ukraine), 260-267. 


(3) I. Kovacic, On the response of purely nonlinear oscillators: An Ateb-type 
solution for motion and an Ateb-type external excitation, International 
Journal of Non-Linear Mechanics 92 (2017), 15-24. 


Section 9.5: This section is based on my work as presented in the paper 


(1) R. E. Mickens, A note on exact finite difference schemes for the differen- 
tial equations satisfied by the Jacobi cosine and sine functions, Journal 
of Difference Equations and Applications 19 (2013), 1042-1047. 


See also the following: 


(2) R. B. Potts, Best difference equation approximation to Duffing’s equa- 
tion, Journal of Austria Mathematics Society 23B (1982), 349-356. 


(3) R. E. Mickens, Nonstandard finite-difference schemes for differential 
equations, Journal of Difference Equations and Applications 8 (2002), 
823-847. 


Sections 9.6 and 9.7: The discussions in these sections are based closely 
on the following papers: 


(1) R. E. Mickens, Comments on the method of harmonic balance, Journal 
of Sound and Vibration 94 (1984), 456-460. 


(2) R. E. Mickens, A generalization of the method of harmonic balance, 
Journal of Sound and Vibration 111 (1986), 515-519. 


(3) R. E. Mickens and D. Semwogerere, Fourier analysis of a rational har- 
monic balance approximation for periodic solutions, Journal of Sound 
and Vibration 195 (1996), 528-530. 


Section 9.8: Good discussions of the general methodology of iteration 
techniques applied to nonlinear oscillators can be found in the following arti- 
cles: 


(1) R. E. Mickens, Iteration procedure for determining approximate solu- 
tions to nonlinear oscillator equations, Journal of Sound and Vibration 
116 (1987), 185-190. 


(2) R. E. Mickens, A generalized iteration procedure for calculating approx- 
imations to periodic solutions of “truly nonlinear oscillators,” Journal 
of Sound and Vibration 287 (2005), 1045-1051. 


APPLICATIONS AND ADVANCED TOPICS 163 


Section 9.9: A general analysis of “truly nonlinear (TNL)” oscillators and 
various techniques to obtain approximations to their solutions is 


(1) R. E. Mickens, Truly Nonlinear Oscillations: Harmonic Balance, Pa- 
rameter Expansions, Iteration, and Averaging Methods (World Scientific, 
London, 2010). 


A rather large set of works has been devoted to the case of fractional power 
nonlinearities. Details and summaries are given in the following two articles: 


. E. Mickens, Oscillations in an x*/° potential, Journal of Sound an 
2) R. E. Mickens, Oscill ee Ld l of Sound and 
Vibration 246 (2001), 375-378. 


(3) H. Gottlieb, Frequencies of oscillators with fractional-power nonlineari- 
ties, Journal of Sound and Vibration 261 (2003), 557-566. 


Examples of such oscillators are 


The oscillator 


has been investigated by Mickens and several other researchers; see, for ex- 
ample, the following: 


(4) R. E. Mickens, Harmonic balance and iteration calculations of periodic 
solutions to 7 +y~! = 0, Journal of Sound and Vibration 306 (2007), 
398-401. 


(5) Johanna Denise Garcia Saldafia, A Qualitative and Quantitative Study of 
Some Planar Differential Equations — A thesis submitted in fulfillment of 
the requirements for the degree of Doctor of Mathematics, Departament 
de Matematiques, Universitat Autonoma de Barcelona (February 2014). 


Ateb-type periodic functions are being used to help in the representation of the 
periodic solutions of strongly nonlinear oscillators. Since the Ateb functions 
cannot be expressed in terms of a finite combination of elementary functions, 
approximations are required. An attempt in this direction is the paper 


(6) I. V. Andrianov, V. I. Olevskyi, and Y. B. Olsevska, Analytical approx- 
imation of periodic Ateb-functions via elementary functions, Proceed- 
ings of the 5th International Conference on Nonlinear Dynamics (ND- 
KhPI2016, September 27-30, 2016, Kharkov, Ukraine), pp. 260-267. 


Taylor & Francis 


Taylor & Francis Group 


http://taylorandfrancis.com 


Chapter 10 
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10.1 GOALS 


The major goals of this book were to present and discuss some prelimi- 
nary investigations on new families of periodic functions and, in some cases, 
also examine the properties of the associated hyperbolic functions. This work 
centered on an examination of the formulas 


IyP+|al?=1, 2(0)=1, y(0)=0, (10.1) 
ly? —|al?=1, 2(0)=0, y(0)=1, (10.2) 

where, in general, 
p>0od, q>Q0. (10.3) 


From our perspective, what we have given is an extension and new interpre- 
tation of previous results obtained on the so-called “generalized trigonometric 
functions.” 

The essential features of our work derive from considering the expressions 
in Equations (10.1) and (10.2) from three different aspects: 


A: as geometrical structures in either the (a — y) or (r — 6) planes, 


B: as a Hamiltonian function for a one-space-dimensional dynamical 
system and 


C: as functional equations. 


10.2 RESULTS 


The following is a summary of the major results obtained from the various 
methodologies used in this book: 

1) Equation (10.1) has periodic solutions. For the geometrical case, where 
the independent variable is 0, 0 < @ < 27, these solutions are periodic, with 
period 27. Further, three classes of periodic solutions arise and these functions 
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are direct generalizations of the standard trigonometric cosine, sine, and dine 
functions. 

2) If the expression in Equation (10.1) is taken to be the Hamiltonian of 
a dynamical system, for which t, the time, is the independent variable, then 
x(t) and y(t) are periodic functions of t. 

3) If now Equation (10.1) is considered to be a functional equation, then 
it has many families of solutions, some of which are periodic. 

4) For Equation (10.2), there are two major families of solutions. First, 
with 0, 0 < @ < 27, as the independent variable, periodic solutions exist, but 
are not defined over specific intervals of x. Also these periodic functions may 
become unbounded for particular values of @. 

5) Considered as a dynamical system, Equation (10.2) can have continuous 
solutions which become unbounded as t > (+)oo. 

6) Many of the most important applications in the natural and engineering 
sciences correspond to 


p=2, q>0, (10.4) 


and this result is a consequence of the Newton definition of force. 
7) The expressions in Equations (10.1) and (10.2) are invariant under the 
transformations 


Ti: 27 -2,y> y, 
To: 07 2,y>-Y, (10.5) 
T3 =T7,T> = ToT) 12> LLY > -Y, 


and, as a consequence, both the periodic and non-periodic solutions have many 
common features imposed on their respective solutions, either x(@) and y(6), 
or x(t) and y(t). 


10.3 SOME UNRESOLVED TOPICS AND ISSUES 


We end the book with a brief presentation of various currently unresolved 
topics and related issues, all of which are of importance in research efforts di- 
rected toward the creation, analysis, and understanding of generalized trigono- 
metric functions. 

(I) Equations (10.1) and (10.2) were originally formulated for real, non- 
negative values of p and q. One consequence of doing this is that absolute- 
values must appear. However, in actual physical universe-related matters, 
where these equations may appear, they can only be known to a finite num- 
ber of numerical places. For such cases, it makes (experimental) sense to have 
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rational representations for p and q. One possibility is to use 
eee 
P= Svat? 
My = 2 Beh (10.6) 
my SO eT Dek a's 


and 
= 2n2 
4 = Bmo41? 


ng = 1,2,3,..., (10.7) 
m2 = 0,1,2,.... 


Note that 


Cs — 2(ni-—mi1)-1 
Dt oa == cere ’ (10.8) 
2n2 _ 1 — 2(ng—me2)—-1 : 

2me2+1 PS Qmetl ? 


goals 


and 


1 
a ea) (10.9) 
(q-1)>0 nz >m2+5 
(II) A major extension of the geometrical approach would be to consider 
general simple closed, convex curves, in the x-y plane, for which the origin lies 
within the curves. We expect the existence of three types of periodic functions 
similar to the trigonometric cosine, sine, and dine functions. However, for an 
arbitrary curve, the details of the analysis would be expected to be much more 
complex than what we have done in our work. 
(III) The geometrical periodic functions associated with 


ey ST, (10.10) 


and 


NA 2 
(=) +y=1, a>l, (10.11) 
a 


differ greatly in the functions they generate. For Equation (10.10), the asso- 
ciated differential equation is 


dx 
du? 


+u=0, (10.12) 


while for the second equation, we have 


d?x 3 
—+Ar+ Br’ =0, (10.13) 
du? 
where A and B are constants, related to a. Inspection of these differential 
equations shows one to be linear, the second to be nonlinear. 
Consider a simple, closed, convex curve, surrounding the origin, taking the 
form 


f(x,y, a) =0, (10.14) 
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where a is a parameter. The results from Equations (10.12) and(10.13) indicate 
that at specific values of the parameter a, a “bifurcation” in the structure of an 
associated differential equation may take place. Is this a general phenomenon? 

(IV) Do “interesting” and useful solutions for the modeling of dynamical 
systems exist for Equation (10.2) which are not (asymptotically) power-law 
or exponential in mathematical structure? 

(V) A general issue: Is there a finite set of families of periodic functions, 
each with a possible infinite number of members? 

(VI) For the case p = 2, q > 1, it would be of value for the techniques such 
as the method of harmonic balance to determine the exact Fourier series of 
the following expressions 

cos 0+ A, cos 30 
1+ B, cos 20+ Bo cos 40? 


cos 0 _ (10.15) 
1+C}, cos 20+C2 sin 20? 


|n|°sgn(x), 


where (A, By, Bz,C, C2) are parameters such that 
ByeBy <i, OP4Oy <4, (10.16) 


and 
a>0. (10.17) 


(VII) For Equation (10.1), with p = 2, can new methodologies be for- 
mulated for the construction of approximations to the exact solutions of the 
equations of motion? (This means that Equation (10.1) is to be taken as a 
Hamiltonian function.) This issue involves the creation of extensions to the 
standard procedures for harmonic balance and iteration schemes. 

(VIII) Can generalized Ateb-type periodic functions be constructed for 


ly? + [alt = 1, 


a(0)=1, ¥(0) =0, 
p2l, q2l? 


NOTES AND REFERENCES 


Section 10.3: Issue (I) is introduced and briefly discussed in the following 
book: 


(1) R. E. Mickens, Truly Nonlinear Oscillations (World Scientific, London, 
2010). 


Note that 
ag 2n — 2n+1 


: PO Friel 


? 
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where (m,n) are non-negative integers having the properties 


(—ax)Pe = xP, (—a)Po = —(2)Po. 


Issue (III) may be related to an essential feature in the formulation of the 
equations of motion for the Hamiltonian formulation, i.e., for example, if 


Hew) = (5) +09) = 5. 


then the equations of motion are 


de OH dy OH _ 


a dy” ad ov 
However, 
dH (x, y) _ dz dy _ 
re a ie 
and the following three possibilities may occur: 
dx dy 
A:—= —=- 
a at a 
dx dy 
B:—=—= =— 
dt dt’ ~ m7 
dx dy 
= t ara t 
Ci = o(2,yty, a= O24 te, 


where ¢(z,y,t) is an arbitrary function of its arguments. Case A is the stan- 
dard answer for the equations of motion. Case B implies x(t) = y(t) = 0, and 
consequently, is not of interest. However, Case C is a “nonstandard” represen- 
tation and reduces to Case A when (x, y,t) = 1. 

Issue (VII) appears in the book by Mickens, see (1) above. 

Issue (VIII) seeks to extend the various investigations of the equation 


y+ |elt@=1, g2>1, 
to the general case 
WP +|2/@=1, p21, q2l. 
See the following references: 


(2) R. M. Rosenberg, The Ateb(h) functions and their properties, Quarterly 
Journal of Applied Mathematics 21 (1963), 37-47. 


(3) A. Beléndez, J. Francés, T. Beléndez, S. Bleda, C. Pascual, and E. Ar- 
ribas, Nonlinear oscillator with power-form elastic form elastic-term: 
Fourier series expansion of the exact solution, Communications in Non- 
linear Science and Numerical Simulation 22 (2015), 134-148. 
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(4) I. Kovacic, On the response of purely nonlinear oscillators: An Ateb-type 
solution for motion and an Ateb-type external excitation, International 
Journal of Non-Linear Mechanics 92 (2017), 15-24. 


APPENDIX 


A TRIGONOMETRIC RELATIONS 
A.1 Basic Properties 


(sin 0)? + (cos@)? = 1 
—l<sinOd<+1, -1l<cos9<+41 
sin(—0) = —sin@, cos(—@) = cos@ 


sin(? + 27) =sin@, cos(0 + 27) = cos0 


sin(0) = 0, cos(0) = 1 
sin (3) =. cos (3) =0 
sin(7) = 0, cos(7) = —1, 
sin (32) =-1l, cos 3m) =0 


ci 4 e-i0 
cos 9 = ————— 


and 
e’”’ =cosé+isind 
e a (e%*)* 
= cosé —isin#@. 


A.3 Functions of Sums of Angles 


sin(@, + 02) = sin 6; cos 2 + cos 6; sin 02 


cos(6; + 62) = cos 6; cos 62 + sin 6; sin 02 
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A.4 Powers of Trigonometric Functions 


sin? @ = (5) (1 — cos 20) 


OQ 
e} 
n 

Ny 
DS 
| 


(1 + cos 20) 


n 
H. 
i) 
w 
DS 
lI 


(3 sin @ — sin 30) 


o) 
e} 
an 

w 
DS 
II 


(3 cos @ + cos 36) 


(3 — 4 cos 20 + cos 40) 


le} 

O° 

DN 

ma 

D 

lI 
COo}r CO] BI BIR NIeR 


Se eee Se 


(3 + 4 cos 20 + cos 40) 


) (10 sin @ — 5 sin 30 + sin 56) 


ia 
JR Sle 


= 
aD 
Sa 


(10 cos 6 + 5 cos36 + cos 56) 


se 

Dp 

aw 

Dp 

| lI 
Log, ORR RON, OI ORNL 


A.5 Other Trigonometric Relations 


+6 6,4 
sinh sesin fy = 2sin (25 2) cos ( aa 


0,+0 0, —90 
e086 + e085 = 20s ( 2) cos ( . :) 


2 


6,4+06 0, —0 
cos) — c0565 = —2sin — 2) sin ( : 5 zi 


) 


sin ral cos Ag = (5) sin(6, a 62) + sin(61 a 62)| 


cos cal sin Ag = (5) sin(6, T 62) a sin(61 = 62)| 


cos 9; cos 62 = 


sin ral sin Ag = 


i 
Nl Mle MIF WIR 
Se Ne 


cos(9; + 62) + cos(@1 — @2)] 


cos(0; + 42) — cos(41 — 42)] 
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A.6 Derivatives and Integrals 


d 
qg (088 = — sind, ap nd = cos? 
2 d2 
gz C8? = — cos, qo2 in? = — sind 


Both cos@ and sin @ are solutions to the following second-order, linear differ- 
ential equation 


d?u(0) 
7 +u(0) =0, 
u(0) = 1, m0) =0, gives u(@) = cos8, 
u(0) = 0, au) =1, gives u(@) = sind. 


Also 
[cos out =sind+C}, 


[su0ae = (—)cos0+ Co, 


where C; and C2 are arbitrary constants. 


A.7 Taylor Series 


22. g2n 
d= —)"—— 
cos DI ) (Oni 
02 04 g2n 
=] paren _yn a 
am Fo 
oP g2n+1 
O= — 
sin DI ) (2n+1)! 
63 6° g2ntl 
=@-—-—+—4+.. —)jr 
oe Ga 
A.8 Other Derived Trigonometric Functions 
sin 6 
tangent : tan(0) = 
cos 8 
cos 6 
cotangent: cotan(@) =— 
sin 6 
secant : sec(6) = 
cos 
cosecant : cosec(0) = 
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B TRIGONOMETRIC HYPERBOLIC RELATIONS 
B.1 Exponential Definition of Hyperbolic Functions 


Hyperbolic sine : sinh(«) = > 
Hyperbolic cosine : cosh(x) = ——— 


Therefore, it follows that 


e* = cosh(«) + sinh(x) 


e~” =cosh(a) — sinh(z). 
The sinh(x) and cosh(z) satisfy the relation 
cosh?(2) — sinh?(a) = 1. 
B.2 Basic Properties 
cosh(—2) = cosh(x), sinh(—a) = — sinh(z) 


cosh(0) =1, sinh(0) = 0 


Lim cosh(a) = +00 
w— oo 


bin sinh(x) = (+)oo 


B.3 Derivatives and Integrals 


d d 
aa cosh(a) = sinh(x), ae sinh(x) = cosh(z) 


ae d? 

i cosh(a) = cosh(x), a? sinh(«) = sinh(z) 

From these two latter relations, it follows that cosh(z) and sinh() are solu- 
tions to the second-order, linear differential equation 


w(x 
wt Ys w(x) =0, 
such that iii 
w(0) =1, wk ) =0, gives w(x) = cosh(z), 
d 
w(0) = 0, wo) =1, gives w(x) = sinh(z). 
x 
Also 


J snh(a)ae = cosh(x#) + Ch, 
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J cosh(e)ae = sinh(x) + Co, 
where C; and C2 are arbitrary constants. 


B.4 Related Hyperbolic Functions 


inh 
Hyperbolic tangent :tanh(«) = sinh(z) 
cosh(x) 

h 
Hyperbolic cotangent : coth(x) = ae 
INN ax 

1 
H boli ia) Se 
adie ag a 1c) 

1 
Hyperbolic cosecant : csch(x) = ane 


B.5 Sum and Subtraction Formulas 


sinh(a + 


y) () sinh(y) 
cosh(x + y) = cosh(x) cosh(y) + sinh(z) sinh(y) 
sinh(« — y) = sinh(a) cosh(y) — cosh(«) sinh(y) 
cosh(a — y) = cosh(x) cosh(y) — sinh(z) sinh(y) 


sinh(z) + sinh(y) = 2sinh € = “) cosh (: 7 “) 
cosh(x) + cosh(y) = 2 cosh (= 5 “) cosh (- 5 “) 
sinh(«) — sinh(y) = 2 cosh (= 5 “) sinh € - “) 


cosh(a) — cosh(y) = 2 sinh (4) sinh (: . ) 


B.6 Relations between Hyperbolic and Trigonometric Functions 

Using the exponential functions definitions of the cosine and sine functions, 
then the following relations exist between the trigonometric and hyperbolic 
functions: 


lI 
3. 
lo) 
ie} 

2 ot 
— 
SS 
8 
Near 
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B.7 Taylor Series 


co gent 
sinh(x => On +i)! 
sere 
rn" ge ye? % rn g2ntl 
aaa Tata (Qn +1) 
ised 2n 
cosh(x) = ys nl 
n=0 
a at ger 


C SPECIAL FUNCTIONS 
C.1 Properties of Even and Odd Functions 
Consider functions f(x) and g(h), defined on a symmetric interval (—a, a), 
a > 0, where a may be unbounded. 
e A function f(x) is an even function on the interval (—a, a), if and only 
if 
f(-a) = f(a). 
e A function f() is an odd function on the interval (—a, a), if and only 
if 
f(-a) = —f(@). 


e Given an arbitrary function g(x) defined on the interval (—a,a), then it 


can be written as 
g(x) = g(x) + g(a), 


where ena 
_g g 
g(a) aa D) ) 
aia Gla) ole) 
g(a) ; 
and 


gP)(-a) = g(a), gO (-2) = -9 (2). 

The functions g'*)(2) and g‘—)(x) are the even and odd parts, respectively, 
of the function g(x). 

The following propositions are true: 

a) If f(x) and g(x) are both even functions, then h(x) = f(x)g(a) is an 
even function. 

b) If f(x) and g(x) are both odd functions, then h(x) = f(x)g(a) is an 
even function. 

c) If f(x) is an even function and g(x) is an odd function, then h(x) = 
f(x)g(x) is an odd function. 
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d) Let f(a) be an even function, integrable on the interval, (—a, a); then 
: f(x)dx = af f(x)da. 
e) If f(a) is an odd function, integrable over the interval (—a, a), then 
: f(x)dx = 0. 


—a 


f) Let f(x) be an even (odd) function defined over the interval (—a, a), 
then its derivative is an odd (even) function on the same interval, i.e., 


f(z) even > a) dd, 
f(x) odd > ae) even. 


C.2 Absolute Value Function 
Let a be a nonzero real number. The absolute value of a is defined to be 


a, ifa>0, 
lal= 40, ifa=0O, 
-a, ifa<0. 
An alternative definition is 
la] = Va?. 
The following rules apply: 
i) ja] > 0 
ii) |—a| =|a 
iii) |ab| = |al|b| 
iv) |¢| = i b£0 
v) |a”| = Jal” 


vi) |a+ }| < |a| + [0], triangle inequality. 
The above rules generalize to the situation where a and 6b are functions. 
A graph of |z|, the absolute value function, is given in Fig. C.2.1. 


178 APPENDIX 


[x| 


x 


Fig. C.2.1 Plot of the absolute value function, ||. 


C.3 Sign and Theta Functions 


The theta and sign functions are defined, respectively, as follows 
1, ifa>0, 
0, ifa <0, 


; 1, if x > 0, 
ae) = -1, ife<0 


A(x) = 


Since both @(x) and sign(x) are discontinuous at x = 0, they are not defined 


at this point. However, depending on the circumstances, the following values 
may be selected 


6(0) = > sign(0) = 0. 


The following relation holds for the theta and sign functions 
sign(x) = O(a) — 0(—2). 
Plots of these functions appear in Fig. C.3.1. Also, note that 


: d|x 
sign(x) = slip 
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sign(x) 


(a) 


Q(x) 


1/2 


(b) 
Fig. C.3.1 (a) Plot of sign(x) vs x. (b) Plot of 0(x) vs x. 
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C.4 Delta Function 

A function, (x), is a “delta function” if it formally satisfies the three 
conditions: 

a) [blade = 0, 

b) d(z) =0, «#0, 

0) f°, 6(@) f(a)de = f(0). 
The use of a “delta function” is based on several important properties: 

d) The “delta function” is even, i.e., 


e) f%, 6(ax) f(x) = (4) FO. 
f) [%_d(w — a) f(w)de = f(a). 


This result follows directly from c). 
g) Let g(x) be a function having a finite number, N, of simple zeros at 
(a1,X2,...,0Nn). We then have 


N 


je ye 


26 Vgi(ax) 


h) x(x) = —d(z). 
i) i 5”) (x) f(x)dx = (=e roy) (0), where 


_ @ f(x) 


dx” 


OM (¢) = 


D FOURIER SERIES 
Let f(x) be defined on —co < x < co and have period 2L, i-e., 


f(a + 2L) = f(x). 


If the following integrals exist, 


2L 
(x) cos (=) dx, (x) sin (=) dx, 
0 


for k = 0,1,2,..., then the formal Fourier series of f(z) on 0 < a < 2D is 
given by the expression 


2L 


0 


where 
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Definition 1: A function f(x) is said to be piecewise continuous on a finite 
interval, a < x < 5b, if this interval can be partitioned into a finite number 
of subintervals such that f(x) is continuous in the interior of each of the 
subintervals and f(x) has finite limits as x approaches either end point of 
each subinterval from its interior. 


Definition 2: A function f(x) is said to be piecewise smooth on a finite 
interval, a < x < b, if both f(x) and f’(«) are piecewise continuous on a < 
x<b. 


Theorem 1: Let f(x) be piecewise smooth on the interval, 0 < x < 2D, then 
its Fourier series is given by 


= k k 
~~ + pe a. cos (=) + by sin (=)| . 


The Fourier series converge at every point x, in the interval 0 < x < 2D to 


the value 
f(a*) + f(a) 
9 y] 
where f(a) is the right-hand limit of the function f at x and f(x~) is the 
left-hand limit of the function f at x. 
If f is continuous at x, then the Fourier series of f at x converges to f(x). 


E HAMILTONIAN DYNAMICS 

Consider a dynamic system having N-degrees-of-freedom (NDOF). As- 
sume that there exists a function H, called the Hamiltonian, which depends 
on 2N + 1 variables, i.e., 


H(x,p,t) = constant, (E.1) 
where 
21(t) pi(t) 
x(t) = oe , p(t) = me (E.2) 
tn (t) p(t) 


The x(t) and p(t) will be called the generalized coordinates and momenta, 
and they both depend on time, t. It is from the Hamiltonian that all the laws 
governing the dynamics of the system must have their genesis. In other words, 
dx(t)/dt and dp(t)/dt must be derivable from a knowledge of H(x,p,t) and 
some other physical and mathematical requirements. 

Experience with the physical universe leads to the following conclusions: 


In general, the outcome of an experiment does not depend on where it is done 
(location) or when it is done (time). 


Since the Hamiltonian is assumed to determine the dynamics, it follows from 
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the second observation that H(«,p,t) must be invariant under an arbitrary 
time translation, to, i-e., 
t—+t+to. (E.3) 


It is sufficient to consider an infinitesimal time translation, t > t + dt. Thus, 
the time translation invariance can be expressed as the condition 


H(x,p,t + dt) = H(2,p,t), E.4) 
ie ANGed 
H(2,p,t) + SEP 64 4 O64] =H(a,n.1), (BS) 
and this implies that 
OH(a,p,t) _ 
aoe =. E.6) 


The conclusion is that the Hamiltonian does not depend on time 
H = H(x,p) = constant. (E.7) 
The first observation of space translation invariance, i.e., 
roa+on, (E.8) 


gives the condition 
A(x Ws 6x, p) = H(z, p), (E.9) 


where dx is the same constant for each coordinate, i.e., 


Li ata, (=1,2,...N). (E.10) 


Comment: One has to be careful here. Certainly, if the «’s are components 
of Cartesian coordinates, then it can always be arranged to have all of the 
separate (x, y, z) “components” be translated by the same magnitude. 

For infinitesimal dz = a, we have 


OH 

H(z,p)+a)> ee + O(a?) = H(a,p), (E.11) 

i=1 : 

and it follows that 
OH (a,p) _ 
» ag Se (E.12) 
Returning to 

H(x,p) = constant, (E.13) 


and taking the total time derivative of H(x,p), gives 


=0. (E.14) 


aH _ x (OH dos | OF dn 
Ox; dt = Op; dt | — 


dt 


i=l 
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This latter equation allows us to determine a connection between the Hamil- 
tonian function, H(x,p), and the evolution equations for the coordinates and 
the momentum. The direct and simplest possibility is to make the choice 

{i = 2Htep) (j= 1,2,..., N) 


dt Op; 
Spit) _ _Obep) (i = 1,2,...,N). 


Li 


(E.15) 


Therefore, given H(a,p), we have a calculational procedure to determine the 
equations of motion. They are 2N, coupled, ordinary differential equations for 
the N-coordinates and N-momenta. 

Equations (E.15) are the Hamiltonian equations of motion. 

Finally, comparing Equations (E.12) and (E.15), we conclude that the 
result in Equation (E.12) is just the conservation of momentum, i.e., 


N N . 
2 aa . wala <- (=) do pilt) =0, (E.16) 


and this implies 


N 
Spi (t) = constant. 
i=1 


NOTES AND REFERENCES 

Sections A and B: Knowledge of the trigonometric periodic and hyperbolic 
functions is ubiquitous and any of a number of the standard textbooks and 
handbooks will provide the full details of the properties of these functions and 
the associated derivations. 

Section C: The delta function was “created” by P. A. M. Dirac and dis- 
cussed (non-mathematically) in this book: 


(i) P. A. M. Dirac, The Principals of Quantum Mechanics, 2nd ed. (The 
Clarendon Press, Oxford 1930). 


A formal examination of such “generalized functions” was done by 
M. J. Lighthill and appears in his book: 


(ii) M. J. Lighthill, Introduction to Fourier Series and Generalized Functions 
(Cambridge University Press, London, 1958). 


For a general, heuristic discussion, see: 


(iii) R. E. Mickens, Mathematical Methods for the Natural and Engineering 
Sciences, 2nd ed. (World Scientific, London, 2017). 


Section D: Two basic, concise, and good introductions to Fourier series 
are the books: 


(iv) S. L. Ross, Differential Equations (Blaisdell, Waltham, MA; 1964). 
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(v) J. D. Logan, Applied Mathematics, A Contemporary Approach, (Wiley- 
Interscience, New York, 1987). 


Section EF: One of the best introductions to the Hamiltonian approach to 
classical mechanics is the book: 


(vi) H. Goldstein, Classical Mechanics (Addison-Wesley; Cambridge, MA; 
1951). 
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